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Introduction 


The  purpose  of  this  monograph  is  to  present  the  theories  of  basic 
tensor  analysis  and  of  the  differential  geometry  of  surfaces  for  the 
purpose  of  formulating  problems  of  coordinate  generation  in  regions 
bounded  by  arbitrary  curves  or  surfaces.  Since  the  writing  of  the  first 
memoir  on  the  subject  of  tensor  analysis  by  Ricci  and  Levi-Civita  [1] 
in  1901  some  very  significant  developments  in  the  theory  of  tensor 
analysis  have  taken  place,  though,  the  major  applications  of  the  subject 
have  only  been  confined  to  the  general  theory  of  relativity  and  to  the 
continuum  mechanics.  In  this  monograph  an  attempt  has  been  made  to 
utilize  the  theories  of  classical  tensor  analysis  and  differential 
geometry  of  surfaces  in  developing  new  methods  for  the  generation  of 
coordinates  in  arbitrary  regions.  Only  those  results  of  tensor  theoretic 
and  differential  geometric  significance  have  been  explained  which  are 
needed  in  the  development  of  the  subject  in  a  fruitful  manner.  However, 
it  turns  out  that  for  a  better  understanding  and  a  sound  conceptual 
orientation  some  basic  ideas,  by  the  way  of  definitions  and  notations, 
have  also  to  be  introduced.  Though  this  elementary  exposition  forms  a 
small  part  of  the  total  effort,  and  is  explained  much  better  in  the 
references  given  below,  nevertheless,  its  inclusion  imparts  a  sort  of 
continuity  to  the  whole  presentation. 

Almost  all  the  material  explained  in  Parts  I  and  II  of  this  monograph 
is  available  in  the  standard  texts,  such  as,  Levi-Civita  [2],  Weatherburn 
[3],  McConnell  [4],  Eisenhart  [5],  [6],  Tolman  [7],  Graustein  [8], 


i 


Synge  and  Schild  [9],  Brand  [10],  Spain  [11],  Truesdell  and  Toupin 
[12],  Struik  [13],  Sokolnikoff  [14],  Willmore  [15],  O'Neill  [16], 
and  Kreyszig  [17],  [18],  on  the  classical  topics  in  tensors  and  dif¬ 
ferential  geometry.  Some  other  texts  and  monographs  which  can  be  used 
with  advantage  are  Aris  [19],  Borisenko  and  Tarapov  [20],  Stoker  [21], 
Spivak  [22],  do  Carmo  [23],  Fltigge  [24],  Howard  [25],  and  Eiseman  [26]. 

Part  III  of  this  monograph  is  the  culmination  of  the  ideas 
developed  in  Parts  I  and  II.  Specifically  §  §  2  and  8  of  Part  I  and 
§  § 2  and  3  of  Part  II  provide  the  necessary  material  for  the  development 
of  new  methods  of  coordinate  generation.  It  is  the  hope  of  the  author 
that  the  material  of  Part  III  will  form  a  framework  for  further  research 
in  the  area  of  mesh  generation  for  physical  problems,  based  on  partial 
differential  equations. 
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Part  I 


Fundamental  Concepts  and  Basic  Tensor  Forms 

§1.  Preliminaries. 

In  this  section  we  summarize  some  elementary  operations  on  vectors 
and  tensors  with  the  assumption  that  an  Euclidean  space  is  available  in 
which  a  set  of  rectangular  Cartesian  coordinates  has  been  introduced. 
Further,  to  maintain  a  sort  of  continuity  of  exposition  with  the  rest 
of  the  sections,  we  also  clarify  the  nomenclature  of  some  commonly  used 
terms.  For  further  details  the  reader  is  referred  to  References  [4], 

[11]  and  [14]. 

In  this  report  the  vectors  and  tensors  will  be  denoted  by  using  the 
symbol  -  under  and  above  a  letter,  respectively.  Thus,  the  vectors  are 
denoted  as  u,  <)>,  etc.,  and  the  tensors  as  T,  x,  etc. 

A  rectangular  Cartesian  system  of  coordinates  in  a  three-dimensional 
Euclidean  space  will  usually  be  denoted  by  (i  =  1,2,3),  or,  occasionally 
as  x,y,z.  The  orientation  of  axes  will  always  be  assumed  to  be  right- 
handed.  . 


Figure  1. 
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The  basis  of  a  rectangular  Cartesian  coordinate  system  will  be  denoted 
by  a  system  of  constant  unit  vectors  e^  (i  =  1,2,3).  The  components 
of  a  vector  v  with  respect  to  a  rectangular  Cartesian  system  will  be 
denoted  by  v  (i  *  1,2,3). 

§1.1.  Summation  Convention  on  Cartesian  Components. 

A  repeated  index  on  quantities  either  appearing  as  a  single  entity 
or  as  products  will  imply  summation.  Thus 


aibi  =  +  a2b2  +  a3b3 


aiTij  "  alTlj  +  a2T2j  +  a3T3j 


Tii  T11  +  T22  +  T33 


while  no  summation  is  implied  in 


T  +  T 

ij  ji 


§1.2.  Vector  Multiplications  Using  Cartesian  Components, 
(i)  Scalar  or  dot  product: 


a  •  b  =  ab  cos0 


where  0  is  the  angle  between  a  and  b,  and  a,  b  are  the  magnitudes 
of  the  vectors  a  and  b  respectively.  Obviously 


a  =  |a|  “  /a  •  a  =  /a^a^ 
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(li)  Cross  product: 


(8) 


a  •  (b  «  c)  =•  b  •  (c  x  a)  =  c  •  (a  x  b) 

(iv)  Vector  triple  product: 

(a  x  b)  x  c  =  b(a  •  c)  -  a(b  •  c)  (9) 

a  x  (b  x  c)  =  b(a  •  c)  -  c(a  •  b)  (10) 

(v)  Lagrange  identity: 

(a  x  b)  •  (c  x  d)  =  (a  •  c)  (b  •  d)  -  (a  •  d)  (b  •  c)  (11) 

51.3.  Placement  of  Indices  (Covariant  and  Contravariant) . 

For  simplicity  consider  a  two-dimensional  rectilinear  but  skew 
coordinate  system  in  a  plane  as  shown  in  Fig.  3a. 


Figures  3a,  b,  c. 

Let  a  vector  w  emanate  from  the  point  0.  We  now  decide  not  to  use  any 
subscripted  variables  on  the  components  of  w  since  we  want  to  develop  a 
consistent  method  of  index  notation.  Obviously  there  are  two  ways  to 
write  the  vector  w  in  a  linear  form. 
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(i)  Parallel  projection: 


Let  A  and  y  be  a  basis  for  the  coordinate  axes  Ox  and  Oy 
respectively  as  shown  in  Fig.  3b.  From  the  tip  of  w  draw  lines 
parallel  to  OX  and  OY  to  have  a  parallelogram  OAPB.  Thus 

w  =  Ap  +  pq  (12) 

where  p  and  q  are  the  components  of  w  with  respect  to  the  basis 
(A,p).  To  find  the  lengths  OA  and  OB,  we  introduce  unit  vectors 

u  =  A/ | A  j  ,  v  =  p/ J  M |  . 

Then 

w  =  u|A |p  +  v | p | q  (13) 

so  that  I A | p  and  |p|q  are  the  respective  parallel  projections  of  w 
on  the  coordinate  axes. 

(ii)  Orthogonal  projection: 

Another  method  of  writing  w  in  a  linear  form  is  to  draw 
perpendicular  lines  PD  and  PC  on  the  coordinate  axes  as  shown  in 
Fig.  3c.  We  now  draw  lines  OX'  and  OY'  parallel  respectively  to 
PD  and  PC.  Obviously  the  axes  OX'  and  OY'  are  perpendicular  to  OY 
and  OX  respectively.  Let  il*  and  \  be  the  basis  for  this  new  coordinate 
system.  Then 


w  =  iji  R  +  X  S 

But  since 


(14) 
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so  that  writing 


we  obtain 


A*X  =  0,p*ij>  =  0 


r  =  w-  A,s  =  w*p 


w  = 


r  + 


A  •  1 1>  p  •  X 


(15) 


Because  of  the  two  possible  linear  representations  of  the 
same  vector  w,  viz.  (13)  and  (15),  it  is  important  to  introduce  a 
new  system  of  labeling.  It  is  a  standard  convention  to  write 


A  =  ^  ,  p  *  a2 


— : _ =  ai  _ =. — 

A  •  ijj  ~  ’  p  •  X 


so  that 


a^  •  a*  =  1  ,  a^  •  a2  =  1  ,  a^  •  a2  =  a2  •  a*  =  0 


l  2 

p  =  w1  ,  q  =  w^ 


r  =  w^  ,  s  =  w^ 


Thus  (13)  and  (15)  can  be  written  as 


w  “  a^w1  +  a2w2 


(16) 
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(17) 


w  =  a^  +  §2W2 

The  quantities  w^  are  called  the  covariant  components,  and  are 
called  the  contravariant  components  of  the  same  vector  w.  Similarly 
the  vectors  a^  and  a1  are  respectively  the  covariant  and  contravariant 
base  vectors.  It  is  easy  to  conclude  that  if  the  axes  OX  and  OY  are 
orthogonal,  thus  forming  a  rectangular  Cartesian  coordinate  system, 
then  there  is  no  distinction  between  the  covariant  and  contravariant 
components. 

§1.4.  Dyads. 

An  indefinite  product  of  two  vectors  a  and  b  written  as 

l  =  ab  (18) 

is  called  a  dyad.  Some  authors  put  the  symbol  ©  between  a  and  b.  It 
is  instructive  to  view  dyads  as  operators  since  their  utility  lies  in 
the  area  of  operations  with  vectors  or  other  dyads.  In  Cartesian  coor¬ 
dinates  we  can  also  write 

$  =  3ibj ' i' j  (19> 

The  dyad  ba  is  the  transpose  of  $  written  as  .  Thus 

=  ba  (20) 

so  that 

0*  =  a ,b .e  e .  (21) 

l  j'j'i 


9 


* 


1 


SI. 4.1.  Operations  with  Dyads. 


The  scalar  product  of  a  dyad  with  a  vector  u  is  a  vector. 

v  =  $  *  u 

In  general  the  vector  w  obtained  by  pre-multiplication 

w  =  u  •  $ 


is  different  from  v.  However,  it  is  easy  to  verify  that 

v  =  |  •  u  =  u  • 

For  two  arbitrary  vectors  v  and  w,  we  have  the  result 

*•  ~  . 

w  •  (<p  •  v)  =  v  •  (<p  •  w) 

The  unit  dyad  is  defined  as 

I  =  e  .e  . 

~i~i 

since  its  dot  product  with  a  vector 

is  again  v,  viz . , 

v  =  I  •  v 

The  scalar  product  of  two  dyads 

$  =  ab  ,  <J/  =  cd 


That  is, 
(22) 

(23) 

(24) 

(25) 

(26) 


(27) 
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is 


<t>  •  ip  =  (ab)  •  (cd)  =  (b  •  c)ad  (28) 

The  double  scalar  or  inner  product  is  defined  as 
<J>  :  =  (ab)  :  (cd) 

=  (a  •  c)(b  •  d)  (29) 

§1.5.  Curvilinear  Coordinates. 

General  curvilinear  coordinates  introduced  either  in  Euclidean  or 
non-Euclidean  spaces  (cf.  §2)  will  always  be  denoted  by  x^.  As  stated 
earlier,  the  rectangular  Cartesian  coordinates  will  be  denoted  by  x^. 

The  general  coordinates  also  form  a  right-handed  system. 

In  general  coordinates  x1,  a  repeated  lower  and  upper  index  on 
quantities  either  appearing  as  a  single  entity  or  as  products  will 
always  imply  summation.  Thus 


A^a .  =  Ala,  +  A?a„  +  A^a_ 
J  i  J  1  J  2  j  3 


but  no  summation  is  implied  in  the  expression 


A1  +  A'?  on  either  i  or  1 . 
J  1 


All  quantities,  with  the  exception  of  x^  and  x  ,  with  subscripts 
are  termed  covariant  components,  while  all  with  superscripts  are  termed 
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contravariant  components.  It  is  customary  in  all  the  standard  works 
cited  before  to  call  and  v*  as  covariant  and  contravariant  vectors 
respectively.  Similar  is  the  case  with  tensors  too.  However,  it  is 
clear  that  they  are  the  respective  components  of  the  same  entity  v  or 
T.  Occasionally  it  is  helpful  to  write  the  entity  form  such  as 

v  =  v1^  =  (30) 


T  =  Tijia.a.  =  T.  aka^  (31) 

-i-j  k£-  ~ 

where  as  described  before,  a*  are  the  contravariant  base  vectors,  or  the 
reciprocal  basis  to  a^.  The  two  bases  are  related  as 


i 

a 


a . 


~3 


6 


i 

j 


where  6 


i 

j 


is  the  mixed  Kronecker  delta  defined  as 


“  1  if  i  =  j 

=  0  if  i  t  j  . 


(32) 


§1.5.1.  Various  Representations  in  Terms  of  a^  and  a^. 

All  quantities  which  follow  certain  transformation  of  coordinate 
rules  are  called  tensors.  Thus  scalars  and  vectors  are  also  tensors  of 
orders^  zero  and  one  respectively.  However,  it  is  customary  to  name  the 
quantities  of  orders  greater  than  or  equal  to  two  as  tensors.  The  order 
of  a  tensor  is  determined  by  the  total  number  of  free  indices  used  in 
+ 

Also  called  ranks. 
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the  description  of  its  components.  Thus  T  ,  is  a  third  order  tensor; 

jK’ 

covariant  of  order  two  and  contravariant  of  order  one.  The  total 
number  of  components  of  a  tensor  T  are  given  by  N™,  where  N  is  the 
space  dimension  and  m  the  order  of  the  tensor.  Thus  T^  has  27  components 
in  a  three-dimensional  space. 

The  dyadic  representation  of  the  unit  tensor  I  (also  called  the 
Idem  tensor)  in  general  coordinates  is 


I  =  a. a 


(33) 


In  eq.  (30),  v  and  v^  are  the  contravariant  and  covariant  components  of 
a  vector  v,  while  in  eq.  (31)  T1-1  are  the  contravariant  components  of  a 
second  order  tensor  T.  The  covariant  components  of  T  are  given  by 


*  ■  VV 


(34) 


A  tensor  is  said  to  be  symmetric  if 


T. .  =  T, 
iJ  ji 


and 


(35) 


In  entity  form  symmetry  implies 

T  =  T*  . 


A  tensor  is  said  to  be  antisymmetric  if 


f 


I 
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and 


(36) 


rJi 


or. 


iji  _  _(j3i 


Note  that  the  transpose  of  the  representation  (34)  is 


T*  =  T. .a^a1 
13  ~  - 


and  of  (31)  is 


T*  =  Tlja.a. 

-j-i 


(37) 


1.6.  Differential  Operations  in  Curvilinear  Coordinates. 

The  continuity  and  differentiability  of  vector  and  tensor  functions 
in  general  coordinates  follow  the  same  rules  as  those  by  functions  of 
real  variables  in  multivariate  calculus.  Thus,  let  (^(x1)  be  a  scalar 
function  of  general  coordinates.  Then  its  first  differential  is  given 
by 


d<j> 


9 

9X1 


dx 


(38) 


where  dx^  are  the  differentials  of  the  coordinates  x^.  Also  since  dx* 
are  the  contravariant  components  of  the  displacement  vector  dr  (cf.  §2.), 
we  have 


dr  =  a^dx 


(39) 
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Scalar  multiplication  with  a''  on  both  sides  of  (39)  and  a  use  of  eq 
gives 


dx1  =  a1  *  dr 


Thus 


d<j>  =  (■  — a"*")  •  dr 
\  l  ' 

3x 


Equation  (41)  defines  the  operator  V  or  grad  as 


dcji  =  (grad  <j>)  •  dr 


where 


.  n  i  3 
grad  =  V  -  a  — - 

3X1 


The  divergence  (div  or  V  •)  of  a  vector  function  is  given  by 


i  3y 

div  v  =  7  •  v  =  a  •  - - 

Sx1 


and  the  curl  or  Vx  as 

i  3v 

curl  v=Vxv=a  x  — y 

3x 


1.6.1.  Gradient  of  Vectors  and  Divergence  of  Tensors. 
Let  v  be  a  vector  function  of  x^,  then 


3v  ■( 
dv  =  — dx 
3xX 


.  (32) 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


(46a) 
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(46b) 


dx 


JY_ 

Sx1 


Using  (40)  in  (46a)  and  (46b) ,  we  get 


3v  i  3v  i 

dv  *  — =r  a  •  dr  =  ( — =t  a  )  •  dr 

3xl  3x 


,  .  i-,  SY  .  ,  i  3v  v 

(dr  •  a  ) — *7  3  dr  •  (a  - r) 

Sx1  3xX 


Thus  there  is  a  duality  in  the  representation  of  grad  v.  It  can  be 
represented  as 


3y  i  i  3v 


In  this  report,  we  take  the  first  representation  to  represent  grad 


3v  i 

grad  v  =  — =-  a 
3x 


then  its  conjugate  is 


>*  _  i 


3v 


(grad  v)  =  a  — V 
3x 


As  is  obvious  from  eqs.  (47),  we  have  the  identity 


dv  =  (grad  v)  •  dr  =  dr  •  (grad  v)v 


In  the  same  manner,  we  define  the  gradient  of  a  tensor  as 


grad  T  =  — 7  a 


3x 


i  - 


(47a) 

(47b) 

either 

v,  i.e. , 

(48) 

(49) 

(50) 

(51) 
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The  divergence  of  a  tensor  is  then  the  trace  of  (51) ,  that  is 


-  3T  i 

div  T  =  — -  •  a  (52) 

3x 

Below  we  list  some  important  vector  and  tensor  formulae  involving 
vectors  and  tensors  under  the  operations  of  grad,  divergence  and  curl, 
[27]. 


(i)  For  two  vectors  u  and  v,  the  divergence  of  the  dyad  uv  is 


div (uv)  =  (grad  u)  •  v  +  (div  v)u 

(53) 

(ii) 

If  f  is  a  scalar,  then 

div (If)  =  grad  f 

(54) 

(iii) 

The  Laplacian  of  a  vector  u  is 

V2i 

u  =*  div(grad  u)  =  grad(div  u)  -  curl(curl  u) 

(55) 

(iv) 

(grad  u)  •  u  =  grad(1j|uj2)  +  (curl  u)  x  u 

(56) 

(v) 

For  two  vectors  u  and  v 

(curl  u)  x  v  =  [grad  u  -  (grad  u)*]  •  v 

(57) 

(vi) 

div (grad  u)*  =  grad (div  u) 

(58) 

(vii) 

grad(u  •  v)  =  (grad  u)  •  v  +  (grad  v)  •  u 

(59a) 

(viii) 

curl(u  x  v)  =  div(uv  -  vu) 

(59b) 

i 


(60) 


(ix)  div[grad  u  -  (grad  u)  ]  =  -curl(curl  u) 


(x) 


I  :  grad  u  =  div  u 


(61) 


(xi)  For  a  tensor  T  and  a  vector  u, 

div (T  •  u)  =  (div  T*)  •  u  +  T*  :  (grad  u) 


(62) 


(xii)  If  T  is  a  symmetric  tensor  and  r  is  the  position  vector,  then 
div(r  x  f)  =  r  x  (div  T)  (63) 

(xiii)  Let  ft  be  a  skewsymmetric  tensor  in  a  three-dimensional  space, 
then 


ftu  =  0  ,  ft22  =  0  ,  ft33  =  0 


^12  =  — ^21  ’  ^13  =  ~ ^31  ’  fi23  ~  “^32 


(64) 


With  ft  we  can  associate  a  vector  w,  such  that 
(xiv)  div  ft  =  -curl  w 


(xv) 


div  ft  =  curl  w 


For  an  arbitrary  vector  v 

(xvi)  ft  •  v  =  w  x  y 


(xvii) 


ft  :  (grad  v)  =  w  •  (curl  v) 


(65) 

(66) 

(67) 

(68) 


For  an  arbitrary  tensor  T 

(xvii')  w  x  T  =  ft  •  T 


(69) 


52.  Euclidean  and  Riemannian  Spaces  . 


t 


Refer  to  58  of  Part  I  for  a  complete  discussion  on  spaces. 
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Spaces  of  various  kinds,  abstract  as  well  as  perceivable,  are 
needed  to  analyze  mathematically  the  basic  nature  of  practically  all 


problems  in  engineering  science.  The  most  widely  studied  is  the 
N 

Euclidean  space  E  of  dimension  N.  We  shall  usually  be  interested  in 
E3  space,  though  most  of  the  results  are  immediately  extensible  to  any 
value  of  N.  The  most  important  property  of  an  Euclidean  space  is  that 
in  this  space  rectangular  Cartesian  coordinates  can  always  be  introduced 
on  a  global  scale.  As  an  aid  to  form  some  intuitive  ideas  about  spaces, 
it  is  worthwhile  to  realize  that  the  two-dimensional  space  formed  by  the 
surface  of  a  sphere  is  non-Euclidean  since  rectangular  Cartesian  coordi¬ 
nates  cannot  be  introduced  in  it  on  a  global  scale. 

In  rectangular  Cartesian  coordinates  the  position  vector  r  of  a  point 
from  the  origin  is  obviously  given  by 

-  =  ?ixi  =  ?ixi  +  ?2x2  +  ?3X3  *  (70) 

This  type  of  global  expression  for  the  position  vector  r  is  not  available 
in  terms  of  general  coordinates  x1  either  in  the  Euclidean  or  non-Euclidean 
spaces.  Nevertheless,  the  infinitesimal  vector  dr,  which  is  the  directed 
segment  between  two  infinitesimally  close  points,  is  fundamental  to  all 
geometric  considerations.  According  to  Lanczos  [28],  the  line  element  ds 
(s  is  the  arc  length)  which  is  the  magnitude  of  dr,  viz., 

ds  =  | dr | 

is  the  fountainhead  of  metrical  geometry.  To  fix  ideas,  let  in  E3  a 


Cartesian  coordinate  system  has  already  been  introduced.  Then  the 


infinitesimal  vector  dr  is  given  by 


dr  =  e.  dx 
~k 


where  for  the  purpose  of  further  comparison  we  have  written  x,  *  xK 

K 


The  magnitude  ds  is  then  given  by 


(ds)2  =  dr  •  dr  =  e  •  e  dx^dx^ 

*“  *"  K.  ~  X. 


x  a  k  A  fc 

6.  dx  dx 
k£ 


(71) 


where  6  is  the  Kronecker  delta, 


6k£  =  1  if  k  “  * 


*  0  if  k  #  £  . 


In  E3  we  now  introduce  a  curvilinear  coordinate  system  x1.  The 

i 


infinitesimal  vector  dr  is  then  a  function  of  x  ,  so  that 

3r  i  i 

dr  =  — —  dx  =  a.dx 

ax1  '1 


(72) 


and 


(ds)2  =  dr  •  dr  =  (a^  •  a^dx^dx^ 


Writing 


?i  *  =  gi) 


(73) 


we  obtain 
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(ds)2  =  g^dxidxJ 


The  coefficients  g..  are  functions  of  x  and  are  called  the  fundamental 
ij 

metric  coefficients  of  the  chosen  space.  The  chosen  space  is  still 
Euclidean  and  fundamentally  its  metric  coefficients  are  the  constants 
6  but  because  of  the  introduction  of  a  curvilinear  coordinate  system 
the  g  are  expressed  as  functions  of  x*.  The  introduction  of  a  curvi- 

"5  "f* 

linear  coordinate  system  in  EJ  does  not  change  the  nature  of  space  . 

It  is  obvious  from  (73)  that  g_  is  symmetric,  i.e.. 


8iJ  =  8Ji 


and  using  the  condition  that  ds  is  an  invariant,  we  can  equate  (71)  and 
(74).  This  equation  immediately  yields  the  expressions  for  the  g  in 
terms  of  the  derivatives  of  the  Cartesian  coordinates  with  respect  to 
the  curvilinear  coordinates. 

The  name  "Euclidean"  for  a  space  is  due  to  the  fact  that  in  this 
space  the  five  axioms  and  five  postulates  of  Euclid  are  assumed  to  be 
true.  Some  important  consequences  on  the  basis  of  these  axioms  and 
postulates  are  summarized  below. 

(i)  The  Pythagorean  theorem  for  right  triangles  can  be  proved  both 
in  the  infinitesimal  and  global  regions. 

(ii)  Possibility  of  introducing  rectangular  Cartesian  coordinates 
both  in  the  infinitesimal  and  global  regions. 

(iii)  Global  parallelism.  That  is,  a  vector  in  E3  or  E°  can  be 


tAccording  to  Lanczos  [28]  p.  236,  this  understanding  was  the  motive  force 
in  overcoming  the  difficulties  which  Einstein  faced  in  the  year  1914  while 
working  on  the  theory  of  general  relativity. 


\ 


displaced  parallel  *•<  itself  on  any  space  curve  without  a 
change  in  magniti  Thus  in  an  Euclidean  space,  a  parallel 
field  of  vectors  is  constant  in  magnitude  and  direction. 


Now  though  the  metric  equation  (74)  has  been  obtained  by  introducing 
curvilinear  coordinates  in  an  Euclidean  space,  we  have  only  used  the 
results  (i)  and  (ii) ,  implicitly  through  Eq.  (72),  on  an  infinitesimal 
basis.  Further,  no  where  in  the  derivation  of  (74)  the  result  (iii)  on 
global  parallelism  has  been  used.  Following  Riemann,  we  now  take  (74) 
as  the  one  and  only  axiom  of  a  geometry  in  which  the  functions  g  are 
arbitrary  but  continuous  and  at  least  twice  continuously  differentiable 
functions  of  the  coordinates  x\  Because  of  the  general  nature  of  such 
g_'s  this  geometry  will  be  non-Euclidean.  However,  the  possibility  of 
introducing  a  rectangular  Cartesian  system  locally  in  this  general  space 
in  an  infinitesimal  region  still  exists.  These  assertions  have  been 
proved  in  §8.2.  Spaces  in  which  the  Euclidean  background  has  been 
deleted  and  the  formula  for  the  metric  is  as  given  in  (74)  are  called 
Riemannian. 

The  purpose  of  the  preceding  two  paragraphs  has  been  to  bring  out 
the  subtle  differences  between  the  Euclidean  and  non-Euclidean  spaces. 

It  so  happens  that  a  majority  of  analytical  constructions  (such  as  the 
g^j  considered  before  and  obtained  essentially  from  Euclidean  considera¬ 
tions)  can  immediately  be  interpreted  in  the  sense  of  a  Riemannian  space. 
This  technique  eliminates  some  of  the  abstractness  surrounding  the  tensor 
theory  and  allows  us  to  obtain  all  the  essential  formulae  of  Riemannian 
geometry  while  essentially  remaining  in  the  Euclidean  space.  This  is  the 
scheme  for  further  development  of  the  subject  in  this  report.  It  must  be 


realized  that  most  of  the  Riemannian  constructions  are  analytic  by 
nature  since  the  human  mind  is  not  capable  of  imagining  a  curved 
surface  of  dimensions  greater  than  two. 


§3.  Fundamental  Tensor  Structures  and  Transformation  Laws. 

The  fundamental  metric  tensor  g  was  earlier  defined  through  the 
use  of  the  base  vectors  as 


g  .  .  =  a .  •  a . 
ij  'i  '3 


(76) 


If  the  contravariant  base  vectors  a  are  multiplied  scalarly  then  we 
define  the  new  second  order  components 


ij  l  i 
r  J  =  a  •  aJ 


(77) 


To  find  the  relation  between  g_  and  g1^  ,  we  first  write  an  arbitrary 
vector  A  both  in  the  covariant  and  contravariant  components 


A  =  A. a1  ,  A  =  A^a. 
i'  '  ~J 


so  that 


A .  a 1  =  AJ  a  , 

i~  -i 


(78) 


Taking  the  dot  product  of  both  sides  of  (78)  by  a^  and  using  (32),  we  get 


Ak  "  8jkA 


(79) 


Similarly,  taking  the  dot  product  with  a  ,  we  get 


» k  ik. 

A  =  g  A, 


(80) 
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Solving  the  set  of  equations  (79)  for  A  ,  we  have 


,  A.G 

,k  i 


ik 


where 


g  =  det(g^) 


ik 

=  g^G  (by  Kramer's  expansion) 


(81) 


(82a) 

(82b) 


ik 

and  G  is  the  cofactor  of  g  in  the  determinant  g.  Comparing  (80)  and 
(81),  the  required  relation  is 


g 


ik 


Gik  cofactor  of  g^^  in  g 
8  g 


(83) 


ik 

The  tensor  components  g  are  called  the  conjugates  of  the  metric  compo¬ 
nents  Equation  (83)  can  also  be  expressed  as 


g 


ij 


(grsgU  -  8rt8is^g 


(84) 


where  the  groups  (i,r,2)  and  (j,s,t)  separately  have  their  indices  in 
the  cyclic  permutations  of  1,2,3,  in  this  order.  Obviously  g1J  are  also 
symmetric  in  i,j. 

Having  defined  g1J  ,  we  find  from  Eqs.  (79)  and  (80)  that  if  the 
contravar iant  components  of  a  vector  are  known  then  the  covariant  compo¬ 


nents  can  be  obtained,  and  vice  versa.  In  (79)  the  index  has  been  lowered 
and  in  (80)  the  index  has  been  raised.  These  operations  are  called  low- 


ering  or  raising  an  index  respectively  through  g^,  and  g 


ij 
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§3.1.  Relations  Between  the  Base  Vectors. 

From  eq.  (32)  we  have  the  fundamental  relation  between  the  co¬ 
variant  and  contravariant  base  vectors,  which  is 


?i 


(85) 


Equation  (85)  shows  that  one  vector  from  the  reciprocal  basis  a'*"  is 
orthogonal  to  two  vectors  from  the  basis  a^.  Thus,  for  example,  x  a^ 
must  be  parallel  to  a1,  so  that 


a1  =  p(a2  x  a^)  (86a) 

where  p  is  a  scalar  function  of  the  coordinates.  Further,  since  a1  •  a^  =  1, 
we  have 


1 

P 


?1 


(a2  x  a^)  . 


Using  the  vector  formula  given  in  (8),  we  also  have 


"p  =  -1  *  X  ~3^  =  ~2  ’  ^'3  x  '1^  =  ~3  ’  (~l  x  ~2^  '  (86b) 


All  the  possible  forms  such  as  in  (86a)  can  therefore  be  written  as 


a  =  p(a^  x  ak)  (87) 


where  i,  j,  k  are  in  the  cyclic  permutations  of  1,  2,  3,  in  this  order. 
Similar  considerations  show  that 


/  ®  n , 
=  «l(a  x  a  ) 


(88) 


ii 
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where  £,  m,  n  are  in  the  cyclic  permutations  of  1,  2,  3,  in  this  order. 
From  (87)  and  (88)  using  (11) ,  we  get 


6 


i 

l 


pqCAJ  -  «“6J) 


n.m. 


(89) 


so  that 


pq  =  1  . 

Further,  on  the  basis  of  the  result  (84)  it  is  easy  to  show,  using  (87) 
and  (11) ,  that 


P 


q  =  • 


Having  obtained  the  values  of  p  and  q,  we  can  rewrite  (87)  and  (88)  in 
the  following  useful  forms 


a. 

-J 


x  5k 


(90a) 


k 

x  a  = 


4 


eijka. 


(90b) 


Note  that  from  (90a, b)  we  also  have 


a1  =  —  eljk(a.  x  a  )  (90c) 

2/i  ~J  ~k 

?i  *  e_R(a-i  x  ak)  (90d) 

i  i  k 

where  e  is  also  a  permutation  symbol  written  in  contravariant  form 
so  as  to  be  consistent  with  the  summation  convention. 

We  now  use  the  rule  of  lowering  and  raising  an  index  to  base  vectors 


It  is  obvious  that 


If  eq.  (91)  is  rewritten  using  different  indices  as 


a . 


g  a 


then  the  dot  product  of  (91)  and  (93)  gives 


k£ 

8ik8j  £8  "  8ij 


xk 

~  8ik6j 


(93) 


From  this  we  obtain  the  important  result 

gj,g  =  (94) 

In  (33),  the  idem  tensor  I  referred  to  general  coordinates  was 
defined.  We  can  also  write  (33)  as 

I  =  6*a.aj  (93) 

J~i~ 

Using  (91)  and  (92),  the  other  two  representations  are 

I  =  gjjfV  (96) 

=  (97) 

The  use  of  base  vectors  also  allows  us  to  write  vectors  and  tensors  in 
the  entity  forms.  The  choice  of  a  particular  form  of  components  is  usually 
dictated  by  the  user  according  to  his  needs.  For  writing  tensors,  the 
following  forms  can  be  used. 


(i)  All  components  in  contravariant  form. 

(ii)  All  components  in  covariant  form. 

(iii)  All  components  in  mixed  covariant-contravariant  form. 

In  eqs.  (30)  and  (31)  we  expressed  a  vector  and  a  second  order 
tensor  in  component  forms.  Thus,  for  a  vector 


v  =  v  a^ 


(98a) 


v.a 

3' 


(98b) 


For  a  tensor  of  the  second  order. 


T  =  T1Ja  a 

~i~j 


(99a) 


V1?3 


(99b) 


=  ^j§iaj 


(99c) 


=  T-Vaj  (99d) 


Using  (92)  for  a1  in  (99b)  and  then  equating  with  (99c)  we  get 

“  8lnTnj  (100) 


Similarly,  using  (91)  for  a^  in  (99a)  and  equating  with  (99d) ,  we  get 

i;J  ■  8tain,J  don 

The  dot  placed  before  an  index  indicates  which  index  has  been  raised  or 
lowered . 
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r 


Tensors  of  various  orders  can  be  written  by  using  as  many  base 
vectors  as  the  order  of  the  tensor.  When  there  is  no  confusion  in 
recognizing  which  symbol  has  been  raised  or  lowered,  we  may  suspend 
the  use  of  dots.  Thus 


rlV 


P 

i 

q 


are  the  mixed  components  of  a  tensor  of  (p  +  q)th  order;  covariant  of 
order  q  and  contravar iant  of  order  p. 

§3.2.  Transformation  Laws  for  Vectors  and  Tensors. 

We  now  consider  the  transformation  laws  for  vectors  and  tensors  under 

i  — i 

a  change  of  the  coordinate  system  x  to  x  .  Implicitly  there  is  a  func¬ 
tional  relation  between  the  two  coordinate  systems,  viz.,  the  coordinates 
are  functions  of  the  coordinates  x^.  Thus 

x1  =  /(xj)  (102) 


We  assume  that  the  mapping  or  transformation  (102)  is  nonsingular  so  that 
the  functions  4> 1  are  continuously  differentiable  and  their  functional 
determinant  (i.e.,  Jacobian) 


J  =  det  (— ^t)  (103) 

9xJ 

is  no  where  vanishing.  This  implies  that  the  functions  (102)  can  be 
inverted  to  have 

x1  =  </(xj)  (104) 
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(i)  Vectors: 

The  simplest  but  fundamental  is  the  vector  dr.  When  the  coordinate 
system  is  x*,  then  as  noted  earlier 


dr  =  — -  dx 

ax1 


=  a  .dx 
~  1 


(105a) 


(105b) 


On  change  of  coordinates  from  x1  to  x1,  the  same  vector  dr  can  be  written 


Obviously 


3r  -i 

dr  =  - r  dx 

ax1 


=  a^dx 


a.dx''  =  a.dx1 
-J  'i 


(106a) 


(106b) 


(107) 


By  the  chain  rule  of  differentiation,  we  have 


a  J  j— i 

dx  =  - r  dx 


(108a) 


i  ax1  j 

dx  =  - dx 


Using  (108a, b)  in  (107),  we  obtain 


ax1  - 
a  .  =  — r  a  , 

3xJ  '1 


(108b) 


(109a) 


—  3xJ 

ai '  ^  sj 


(109b) 


A  study  of  (108a)  and  (109a)  or  (108b)  and  (109b)  is  revealing.  It 

suggests  that  the  set  of  quantities  dx^  and  set  of  vectors  a^  follow 

different  transformation  laws  on  changing  from  x^  to  x^  or  vice  versa. 

Quantities  which  transform  in  the  manner  of  (108)  are  called  the  contra- 

variant  components  while  those  which  transform  in  the  manner  of  (109)  are 

called  the  covariant  components.  Thus  dr  is  a  vector  whose  contravariant 

components  are  dx^,  while  the  vectors  a^  are  called  the  covariant  base 

vectors.  (This  was  the  reason  for  denoting  the  coordinates  as  x*  and 

base  vectors  as  a^) .  Another  simple  vector  which  has  covariant  components 

in  a  natural  way  is  the  gradient  of  a  scalar.  If  f  is  a  scalar,  then  its 

i  3  f 

first  partial  derivatives  with  respect  to  x  are  — -  .  On  changing  the 

-i  9X 

coordinates  to  x  ,  we  have 


3f  =  3x2  JL 

3x^  3x^  Sx"*- 


(110) 


which  is  exactly  of  the  form  (109a).  Thus  the  first  partial  derivatives 
of  a  scalar  form  the  covariant  components  of  the  vector  grad  f. 

Based  on  the  above  deductions  we  now  state  the  transformation  laws 
for  any  vector  A. 

i  — i 

The  contravariant  components  A  of  the  vector  A  change  to  A  on  a 
change  of  coordinates  from  x*  to  x^  according  to  the  laws 


a1  =  is!  ^ 


(111a) 
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(111b) 


Similarly,  the  covariant  components  A^  transform  according  to  the  laws 


A 


i 


A. 

3X1  J 


(112a) 


A.  =  ~  A  (112b) 

J  3xJ  i 

It  must  be  noted  that  the  transformation  laws  for  vectors  are  linear. 

That  is,  the  vector  components  in  the  new  system  are  linear  functions 
of  the  vector  components  in  the  old  system. 

(ii)  Tensors : 

Consider  a  second  order  tensor  T.  Because  of  the  tensor  invariance, 
we  again  have 

T  =  Tijaa  (113a) 

-i~J 

=  (113b) 


Using  the  transformation  law  of  base  vectors  (109a)  in  (113a),  we  easily 
get 


=  Sjl!  x1:),,' 

3X1  3xj 


(114) 


Similarly 


Equations  (114)  and  (115)  are  the  transformation  laws  for  the  contra- 


variant  components  of  a  tensor 

For  the  covariant  components  we  again  have  from  the  invariance 
condition 


T  =  T. ,aXa 
ij 


J  =  T  iV 

L\a~  5 


Taking  the  dot  product  with  a^,  we  get 


T.  .aV  •  a 
ij~  '  ~P 


T  ^ 

1,  6 

k£~  p 


But 


a 

'P 


a 

'm 


3x 


m 


3? 


so  that 


T. 

im 


- 

3x 

3X7 


i 

a 


=  T  a 
kp~ 


Taking  another  dot  product  with  a^,  we  obtain 


T 

np 


3xS  3xm 
Sx*1  3?* 


T 

sm 


Similarly 


T 

sm 


3^ 


3x 


T 

np 


(116a) 


(116b) 


Equations  (116)  are  the  transformation  laws  for  the  covariant  components 
of  a  tensor. 
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Following  the  same  procedure,  we  list  the  transformation  laws  for 
the  mixed  components.  Starting  from 


T^a  aj 
•J~i~ 


- Z 

T  oai  a 


we  get 


3x^  9x^ 
a?5  axs 


(117a) 


a?3  axs 

3xm  3x^ 


(H7b) 


Similarly,  it  can  be  proved  that  the  transformation  law  for  the  components 


follow  the  same  rules  as  given  in  (117). 


Generally,  we  then  have 


i  a  i 

—  1  2  p  ,3x  3x  3x 

m.m  •••m  i.  i  m 

12  q  3x  1  3x  P  3x  1 


„  Jq  i. i„* • -i 
3x\t  12  p 

& 


(118) 


Metric  Tensor: 

Because  of  the  special  status  of  the  metric  components  ,  and  , 

we  consider  them  in  detail.  Looking  from  the  point  of  view  of  the 
definitions 


ij  i  j 

g  =  §  *  r  ,  =>  •  a 


we  immediately  conclude  that  they  are  symmetric  in  i  and  j.  If  base 


^  ***'■»-  - 
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vectors  are  not  brought  into  picture  and  the  g^'s  are  assumed  to  be 
functions  defining  a  metric  in  a  Riemannian  space  then  we  can  use  the 
formula 


(ds)2  =  g^dx^dx^ 


to  write 


(ds) 


2  =  +  8 j i ) dx Xdx*^  +  yCg^  -  gj;L)dx  dxJ 


(119) 


By  direct  expansion,  we  can  show  that  the  last  term  in  (119)  is  zero, 
proving  the  symmetry  of  g .  Similarly  g*^  is  also  symmetric. 

The  components  g^  are  covariant,  while  g^  are  the  contravariant 
components  of  the  metric  tensor.  Thus,  the  transformation  laws  for 
them  are 


-  _  3xf.  3x® 

5np  "  a?1  a?5  8sm 


(120) 


3x^  3x? 
3xS  3xm 


sm 

g 


(121) 


All  the  preceding  transformation  laws  are  linear.  We  therefore  list 
the  following  important  conclusions  regarding  the  nature  of  tensors. 

(I)  A  tensor  equation  or  expression  has  the  property  that  it  can 


be  obtained  in  any  legitimate  reference  system,  i.e.,  J  ^  0.  If  it  is 
correct  in  any  one  reference  system  then  it  must  remain  correct  in  any 
other  legitimate  reference  system.  The  above  property  is  due  to  the 


linearity  of  the  transformation  laws,  since  any  component  from  the  old 
system  is  a  linear  function  of  the  components  from  the  new  system. 


(II)  If  all  the  tensor  components  vanish  in  any  reference  system 
then  they  remain  zero  in  any  other  system. 

(III)  Because  of  the  linearity  of  transformation  laws,  a  symmetric 
tensor  remains  symmetric  on  coordinate  transformation.  Thus  symmetry 

is  an  absolute  property. 

§3.3.  Algebraic  Properties  of  Tensors. 

•ff 

1.  The  components  of  two  tensors  of  the  same  order  and  structure 
can  be  added  and  subtracted  according  to  the  usual  arithmetical  rules. 

If  A„  and  B_  are  the  covariant  components  of  the  tensors  A  and  B,  then 
on  addition  or  subtraction  we  generate  a  new  tensor  C  whose  covariant 
components  are 


C 


ij 


A.  .  ±  B. . 
iJ  iJ 


(122) 


2.  The  outer  product  of  two  tensors  of  any  order  or  structure  is 

obtained  by  arithmetical  multiplication  to  produce  new  tensors.  For 

example,  let  A^^  be  a  fourth  order  tensor,  and  let  Bm  be  a  third  order 

i  pq 

tensor.  Then  the  outer  product  is  the  new  tensor 


epq 


.  ii  km 

V  B 


pq 


(123) 


which  is  a  seventh  order  tensor,  contravariant  of  order  four  and  covariant 
of  order  three. 


* 

Variance. 
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3.  The  inner  product  of  two  tensors  is  obtained  by  equating  one 
index  of  the  first  with  one  index  of  opposite  variance  of  the  second, 
and,  then  summing  over  this  index.  Thus  in  (123)  if  we  set  l  -  m  and 
sum  over  m,  we  obtain 


Cljk  =  AijkBm 
pq  m  pq 


The  resulting  tensor  is  only  of  the  fifth  order. 


(124) 


§4.  Differentiation  of  Vectors  and  Tensors. 

One  begins  to  feel  the  power  of  the  method  of  tensor  analysis  after 
the  differential  aspects  of  tensors  have  been  completed.  In  this  connection 
we  first  consider  the  partial  derivatives  of  the  base  vectors. 

From  the  definition  of  base  vectors  a we  first  note  the  following 
result. 


3a.  s  3r  a  3r 

k  k'  i'  i'  k/ 

3x  3x  3x  3X  3x 


(125) 


for  any  values  of  i  and  k. 

We  now  select  any  three  indices,  say  i,j,k,  and  consider  the  following 


three  equations 


3g.,  a 

— =  — —  (o  •  a  ) 

,  i  .  i ~i  V 

3x  3x  J 


(127) 


°6ik  3  ,  , 

1?  ‘  I?(;‘  ‘ ^ 


(128) 


Adding  (127)  and  (128),  subtracting  (126)  from  it  and  using  (125),  we 


3a. 

3?  •^■U3’kI 


(129) 


where 


1  3®ik  ^®ik 
[ij,k]  =  + 

3xJ  3xX 


(130) 


The  quantities  [ij,k]  defined  in  (130)  are  called  the  Christoff el  symbols  of 
the  first  kind.  These  quantities  are  symmetric  in  i  and  j  but  they  are 
not  tensors. 


Equation  (129)  implies 


3~i  k 

— r  =  [ij,k]ak 

3xJ 


(131) 


Taking  the  dot  product  both  sides  by  a  ,  we  then  obtain 


3'i  i  ri 
— r  •  a  =  r. , 

3xJ  '  lj 


(132) 


where 


are  called  the  Christoffel  symbols  of  the  second  kind.  These  quantities 
are  also  symmetric  in  the  lower  two  indices  i  and  j,  but  are  not  tensors. 
Equation  (132)  implies 


3a . 


3x 


j 


riA 


(134) 


Equation  (134)  expresses  the  first  partial  derivatives  of  the  covariant 
base  vectors  in  terms  of  the  derivatives  of  g 

ij 

Having  established  the  preceding  definitions  and  results,  we  now 
consider  the  partial  derivatives  of  an  arbitrary  vector  A. 

Let  A  be  an  arbitrary  vector,  and  we  express  it  in  terms  of  its 
contravariant  components  A^ .  Thus 

A  =  A1ai  (135) 

Differentiating  with  respect  to  x  ,  we  get 


3A 


a.  +  A 


3x 


k  ~i 


3xk 


On  using  (134)  and  adjusting  the  dummy  indices,  we  get 


_  ,3A 
1,  ~  ' 


3x 


3x 


k+  rjVJ)5l 


We  use  a  special  notation  for  the  terms  in  parentheses, 


a1  =  -A  + 

,k  3xk  jk 


(136) 


which  is  called  the  covariant  derivative  of  the  contravariant  components 
A1 .  Thus 


k  =  A,k-i 


(137) 


To  find  the  covariant  derivative  of  the  covariant  components  A^  of 
the  vector  A,  we  start  differentia  Ing 


A  =  A. a 

i~ 


so  that 


3A  3Ai  i  8a1 

~  ?  +  Ai  ~ 

3x  3x  1  3x 


(138) 


To  obtain  — -  ,  we  differentiate  the  relation 
3xk 


i  x  i 

a  •  a .  =  6  . 
~j  J 


yielding 


3a  i  3?i 

~T  ’  +  -  ’  k  =  0 

3xk  J  3xk 


(139) 


Using  (134)  in  (139),  we  get 


—  •  a 

k  -j 

i r  J 


1 


which  implies 


-rjV 


Using  (140)  in  (138) ,  the  derivative  becomes 


(140) 


3A 


3x 


k 


(141) 


where 


A.  U 
i,k 


3A . 

— t  ~  H.a. 

3xk  lk  J 


(142) 


is  called  the  covariant  derivative  of  the  covariant  components  A^. 

The  comma  notation  will  always  imply  the  operation  of  covariant 
differentiation.  The  name  "covariant"  for  this  type  of  differentiation 
is  due  to  the  fact  that  the  differentiated  component  gains  one  covariant 
index  with  each  application.  For  example,  in  (142),  the  covariance  of 
che  components  is  now  of  the  second  order  due  to  the  covariant  differentia 
tion.  Similarly,  A1  is  a  mixed  tensor,  contravar iant  of  order  one,  and 

9 

covariant  of  order  one  due  to  differentiation. 

Following  the  method  described  in  §3.2,  we  can  introduce  a  transforma 
tion  from  to  x^  in  (136)  and  (142)  to  have 


(144) 


where 


3X"1  3xn 


A. 


i>£  3X1  3x*  j*n 


A1  -a^l+r1^ 


3x 


rn 


(145) 


3A. 

_ Lfl 

3x*  **  * 


(146) 


A  bar  on  the  quantities  in  (145)  and  (146)  denotes  values  in  the  new  coor¬ 
dinate  system  x1.  The  transformation  equations  (143)  and  (144)  prove  that 
A1  and  A.  are  tensors  of  the  stated  structures. 

,  K.  i,  k 

An  important  point  to  be  noted  is  that  on  comparison  of  (136)  with 
(145)  and  of  (142)  with  (146),  there  is  found  no  change  in  the  forms  of 
the  covariant  derivatives  in  changing  the  coordinates  from  x*  to  x^.  Thus 
there  is  no  preference  of  one  coordinate  system  over  any  other  as  far  as 
the  covariant  differentiation  is  concerned.  Also  in  the  case  of  rectangular 
Cartesian  coordinates,  since  the  Christoffel  symbols  are  zero,  the  covariant 
differentiation  reduces  to  partial  differentiation.  Because  of  these  pro¬ 
perties,  the  covariant  differentiation  is  also  called  "absolute  differentia¬ 
tion." 

In  obtaining  the  partial  derivatives  of  a  tensor  as  an  entity,  we 
again  encounter  the  covariant  derivatives  of  the  components  in  which  the 
tensor  has  been  expressed.  Thus,  for  a  second  order  tensor  written  in 
contravariant  components,  we  have 


j  _  j  1J 
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r 


Ricci's  Theorem:  The  covariant  derivatives  of  the  metric  tensor  g  , 

,  or  6^  are  identically  zero. 

This  theorem  can  be  proved  by  replacing  T  by  g  in  (151)  and  (152) 
and  using  the  expression  for  the  Christoffel  symbols  given  in  (133). 


=  0  ,  gij>k  =  0  ,  sjjk  =  0 


(155) 


and  the  metric  coefficients  behave  like  constants  under  covariant  differ¬ 
entiation.  Because  of  this  property,  e.g., 


(e^T  )  = 

V8  mn  ,k  8  mn,k 


(156a) 


From  (155) ,  we  have 


38ii  = 

3xk  rik8«.j  +  rjk8iH 


(156b) 


=  -r1  g£j  _  rj  gi£ 

.  k  £k8  *£k8 

3x 


(156c) 


§5.  Christoffel  Symbols:  Their  Properties  and  Transformation  Laws. 

The  definitions  of  the  Christoffel  symbols  of  the  first  and  second  kinds 
have  already  been  given  in  (130)  and  (133)  respectively.  It  must  be  restated 
that  these  symbols  are  not  the  components  of  any  tensor.  The  transformation 
laws  considered  in  this  section  will  prove  this  assertion. 

In  taking  the  divergences  of  vectors  and  tensors,  a  contracted 
Christoffel  symbol  of  the  form  appears.  To  find  its  value,  we  use 
equation  (82b) , 


8  =  g*mG 


im 


where  G  is  the  cofactor  of  g„  in  the  determinant  g.  Thus 

m 


3k 

3g  _  gJlm  B£m 

3x^  3x^ 


(157) 


and  since 


_£m  Jim 
G  =  gg 


(158) 


we  get 


3g 

3g  Jim  BJlm 

j  “  88  - T 

3xJ  3xJ 


(159) 


Now  in  (133)  setting  Jl  =  i,  and  summing  over  i  while  using  the  property 
that  g_  is  symmetric,  we  get 


_i  1  ik  3gik 

r .  .  =  -r  g  — r 

^  2  3xj 


On  using  (159),  we  have 


r1  .-LJt 

«  28  3xJ 


(160a) 


=  ~(£n/g)  (160b) 

3xJ 


15.1.  Transformation  Laws  for  Christoffel  Symbols. 

Let  x1"  and  x*  be  two  general  coordinate  systems.  We  assume  that  x^ 
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i  r-r  r~ 


and  are  functionally  related  and  that  the  Jacobian  of  the  transforma¬ 
tion  is  not  zero.  Recall  from  (120)  and  (121)  that  the  metric  components 
transform  as 


3xk  3x^ 
§ij  '  ^  3X1  3? 


(161) 


-ij  k£  3X1 

8=8  k  — I 

3x*  3x* 


(162) 


From  (130)  and  (133),  the  Christoffel  symbols  of  the  first  and  second 
kinds  respectively  for  the  coordinates  x^  are 


rT^Ti  l/3gik  .  3gjk  3gi1. 
[ij  ,k]  =  -=■( — ; -r  +  —±t - =]f) 

3xJ  3x  3x 


(163) 


— k£  r~ 


Ty  =  g  [ij  ,k] 


(164) 


If  we  now  use  (161)  in  (163)  and  perform  the  indicated  differentiations, 
we  get 


Sx1  3xJ  3xk 


[«.m,n]  =  [ij  ,k]  — -  — ~  +  g 


3X1  32xJ 


—l  — m  n  ij  — n  — £.-m 

3x  3x  3x  3x  3x  3x 


(163) 


Inner  multiplication  by  g*1*3  (given  in  (162)),  gives 


fP  -  Ts  3x~^  32x~^ 

=  ij  3xS  3x*  3^  3xj  3x^3^ 


(166) 


Equations  (165)  and  (166)  are  the  transformation  laws  for  the  Christoffel 
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symbols.  Because  of  the  appearance  of  the  second  derivatives  of  the 

coordinates  on  the  right  of  eqs.  (165)  and  (166),  the  symbols  do  not 

transform  like  the  components  of  any  tensor.  This  proves  that  the 

Christoff el  symbols  are  not  tensors. 

A  formula  expressing  the  second  derivatives  of  coordinates  can  be 

obtained  from  (166).  On  taking  the  inner  multiplication  of  (166)  by 
3x* 

— —  ,  we  get 
3xp 


32xr  _  yp  3x^  _  r  Sx*  3x~^ 

a— JL-m  £m  p  1  ij  —i  .-m 

3x  3x  3x  J  3x  3x 


(167) 


§5.1.1.  Formulae:'  Cartesian  to  Curvilinear  and  Vice  Versa. 

All  the  preceding  formulae  are  applicable  for  any  space  and  for  any 
two  general  coordinate  systems.  In  engineering  applications,  we  usually 
transform  from  a  rectangular  Cartesian  to  a  curvilinear  and  vice  versa. 

We  consider  two  cases. 

i  — ^ 

(i)  x  are  Cartesian  and  x  curvilinear. 

(ii)  x1  are  Cartesian  and  x1  curvilinear. 


Case  (i) : 

If  x1  are  Cartesian,  we  denote  x^  as  decided  earlier,  by  x^.  For 
this  case 


k«, 


’kZ 


=  6 


kZ 


so  that 


(168) 


For  brevity  of  notation,  we  denote  x'*'  =  £*,  and  also  remove  the  over 
bar  from  the  quantities  in  £* .  Thus 


3x,  3x 

g..  =  - : - -  (sum  on  k)  (169) 

3?1  3CJ 


,ij  =  ist  i£ 

3xk  3xk 


(170) 


To  find  the  partial  derivatives  of  the  curvilinear  coordinates  with 

respect  to  the  Cartesian,  take  the  inner  multiplication  of  (170)  with 
3x 


r 


3Cj 


Thus 


3x 


r 


. .  3x 

-  81J  -J 
3 


(171) 


Recall,  from  (94)  that 


gij8 


(172) 


From  (167),  the  second  partial  derivatives  are  given  by 


3  xr 


=  ^ 


dx 


U 


(173) 


Inner  multiplication  of  (173)  with  -r-*~  yields  the  formula  for  the 

OX 
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Christoffel  symbols  in  terms  of  the  second  derivatives  of  the  Cartesian 
with  respect  to  the  curvilinear  coordinates 


82> 


ij  3xs  3^3^ 


(174a) 


=  8  — 7  — i — 7  (sum  on  s  and  c)  (174b) 

H  30CJ 

As  an  application  of  the  preceding  results  in  a  two-dimensional  plane 
in  which  the  Cartesian  coordinates  are  =  x,  y±  =  y  and  the  introduced 
curvilinear  coordinates  are  S1  =  £,  £2  =  n,  we  have  the  following  formulae 
In  all  the  formulae  given  below,  a  variable  subscript  denotes  a  partial 
derivative. 


gll  =  x|  +  y|  *  812  =  Vn  +  Vn  »  822  =  xn  +  yn 


8 11  ^  +  Cy  ’  g12  -  +  K  n  ,  g—  =  +  iy* 

Ay  xxyy  x  y 


,22  =  n2  ,  2 


(175) 


11  =  122  12  =  ^12  22  _  811 


g**  =~  ,  g 


where 


g  "  gng22  “  (gi2)2  =  (Vn  '  Vc>2 


5  (C  n  -  n  (  )"2 
x  y  x  y 


~xn/’/g  .  nx  — y^/7g  ,  =  x^//g 


ri\  =  Eg' 


9g 


11  .  .Sgll  „  3g12, 


11  -22  +g12(^-2^f)l/28 


(176a) 

(176b) 

(177) 


(178a) 


(179b) 


Also 


=  (xry  ~  )!'/g 

(  nn  'C  no  b 


r1 

22 


r  n  3g12  9g22. 

IS99C2  1  ^  ™  r^r  )  ""  8 


3g 


22 


’22  3n 


ac 


12  3n 


]/2g 


=  (y  x  -xy  )/v^ 

n  nn  n7on  6 


r2  = 

11 


f  C2  3g12  agll>  98 

l8llV  3£  ~  3n  )  812  3? 


11 


]/2g 


'  <V«  -  ytx£S)/4 


r 1  =  r1 

*12  ‘21 


(a  38H  9g22'  , 

g22  3n  “  g12  3£  )/2g 


(ynX£n  "  Vcn)//g 


r2  =  r2 
12  21 


,  3g22  3gll 

^811  3£  "  g12  3n 


)/2g 


=  (Vcn  "  y£xCn)//g 


(180a) 

(180b) 

(181a) 

(181b) 


(182a) 

(182b) 


(183a) 

(183b) 


,  2  _ 

1 1 2  ~ 


_L  is 

2g  ac 


r22 


±  M 

2g  3  M 


(184a) 
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12 


(184b) 


Case  (ii): 

If  we  treat  x*  as  the  Cartesian  coordinates,  then  we  denote  them  as 
x^.  Also  writing  x1  =  E, 1  in  (167)  we  get 


32Cr  _  _rr  3C1  3Cj 

3x„3x  ij  3x.  3x 

Z  m  J  Z  m 

On  using  (171),  we  can  write  (185)  as 


(185) 


dhr 

3x  3x 
Z  m 


.r  1, 

lj  Hq 


(186a) 


or 


3x„  3x 

_rr  giP-Jq  — 1— “ 
lj  3CP  3Cq 


(186b) 


Equation  (186b)  expresses  the  second  partial  derivatives  of  general  coor¬ 
dinates  in  terms  of  the  first  partial  derivatives  of  the  Cartesian  with 
respect  to  the  general  coordinates.  Equation  (186a)  is  suitable  for 
obtaining  the  Laplacian  of  the  general  coordinates.  For,  on  contracting 
the  indices  l  and  m,  viz.,  setting  t.  =  m  and  performing  the  sum  on  m,  we 
get 


-v 


=  -g 


ijrr. 

ij 


(187) 


where 
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Vz  «  - — r —  (sum  on  m)  . 
3x  3x 
m  m 


Thus  in  two  dimensions,  writing  x1  =  x,  x2  =  y,  51  38  £,  C2  =  n  we  have 


v2S  “  (2gi2ri2  '  gllP22  "  g22Vll)/g 


v2n  -  (2812ri2  "  glir22  "  822ril)/g 


(188) 


(189) 


where 


V2  =  -—j  +  ~~~Z 
3x  3y2 


A  second  order  differential  operator  defined  as 


D  '  822  Jj?  "  2g12  3£3n  +  811 


(190) 


and  the  use  of  eqs.  (173)-(18l)  yields  another  form  of  the  Laplacians, 


V2C  =  (x^Dy  -  y^DxJ/g3/2 


V2n  =  (yrDx  -  x  Dy)/g3/2 


(191) 


(192) 


§6.  Gradient,  Divergence,  Curl,  and  Laplacian 
(i)  Scalars: 

There  are  two  types  of  scalar  quantities.  One  is  called  an  absolute 
scalar  or  an  invariant ,  while  the  other  is  called  a  scalar  density. 

Any  function  of  the  coordinates  x1  is  called  an  absolute  scalar  if 
on  coordinate  transformation  from  x^  to  x^  the  value  of  4>  does  not  change. 
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(j>  (x1  ,x2,x3)  =  (Hx^x^x3)  (193) 

There  are  scalars  which  on  coordinate  transformation  do  not  trans¬ 
form  like  (193).  As  an  example  the  function  g,  viz., 

g  =  det  (g„ ) 

is  not  an  absolute  scalar.  On  coordinate  transformation 

g  =  det(g  ) 

On  actual  substitution  of  (161)  in  the  above  determinant  and  by  expansion, 
we  obtain 


g  =  (J)2  g  (194) 

Similarly 

g  =  (J)2g  (195) 

where 

3  i  _  a— i 

J  =  det  (~~r)  ,  J  =  detC^r)  (196) 

3xJ  3xJ 

Thus  g  or  /g  is  not  an  absolute  scalar,  its  value  in  some  other  coordinate 
system  is  given  by  j/g. 

Multiplying  the  absolute  scalars,  vectors,  or  tensors,  by  /g  we  get 
the  corresponding  densities. 

In  §1.6  we  have  already  defined  the  operator  V  or  grad  as 
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.  o  i  3 

grad  =  V  =  a  - r 

~  ~  ~  i 

dX 


If  <j>  is  an  absolute  scalar,  then  grad  4>  is  a  vector  given  as 

grad  $  =  a'*’  -^r 
Sx1 

so  that  the  covariant  components  of  grad  <j>  are 

(grad  $)  =  -^r  =  4>  ,  (198) 

3X1 

Using  the  method  of  raising  an  index  (cf.  §3),  the  contravariant  com¬ 
ponents  of  grad  are 

(grad  (p) 1  =  gik  -^j-  (199) 

8x 


(ii)  Vectors: 

The  divergence  of  a  vector  v  was  defined  in  (44)  as 


div 


On  using  (137),  we  get 


div  v 


-i) 


-  V 


(200) 
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■Ink 


l/Lit  ~i 


From  (136) 
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kj  1 

g  v  ,a.a.  ,  contravariant  components. 
>  K-~  J 


(205d) 


In  the  same  manner  (grad  v)*  can  be  written.  For  example,  one  useful 
representation  is 


(grad  v)*  =  v1  aka . 

9  1 


(206) 


In  mechanics,  we  sometimes  need  the  inner  products  of  grad  v.  Using 
the  definition  (29),  we  easily  obtain 


x  k 

(grad  v)  :  (grad  v)  =  v  v  . 

,  k.  ,  l 


(207) 


(grad  v)  :  (grad  v)*  =  g.  g^v1  vm 
'  '  lm  ,  k.  ,  n 


(208) 


where  both  are  scalars. 

The  curl  of  a  vector  v  is  defined  in  the  usual  way. 


curl  v  =  a"1  x 


3v 


3xJ 


,  j  k. 

=  vkj<?  *  5  ) 


Using  (90b),  we  have 


i  1  ijk 

curl  v  =  —  e  v,  .  a  . 


(209) 


Thus  the  contravariant  components  of  curl  v  are  given  by 


,  3v  3v 

(curl  y)1  =  — ( — r - £) 

Tg'  3xJ  3x 


(210) 
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I 

i 

where  i,  j,  k  are  in  the  cyclic  permutations  of  1,  2,  3,  in  this 
order. 

(iii)  Tensors : 

The  divergence  of  a  tensor  of  second  order  has  been  defined  as 
(cf.  eq.  (52)) 


div  T  = 


k 

a 


Using  the  derivatives  given  in  (147) -(149)  we  obtain 


div  T  =  T  ,  a .  (211a) 

,k~i 

■g^u.k?1  <211W 

-gJVj-k.i  (211c) 

Thus  div  T  is  a  vector  whose  contravar iant  components  are  given  by  (211a, 
c)  and  the  covariant  components  are  given  by  (211b).  The  operation  of 
divergence  thus  reduces  the  order  of  the  tensor  by  one. 

For  a  divergence-free  tensor 


T 


ik 


=  0 


(212a) 


or, 


gjkT..  .  =  0 


(212b) 


If  a  tensor  is  such  that 


T*.  .  =  0  (212d) 

*3  >  1 

then  it  is  called  a  covariant  divergence-free  tensor. 

(iv)  Laplacian  of  a  scalar: 

The  Laplacian  of  an  absolute  scalar  tp  can  now  be  obtained  by  first 
using  the-  formula  for  the  covariant  derivative  of  covariant  components 
of  a  vector,  viz.,  (142),  to  have 


(<P  .) 
>  i 


32<j> 


3x13xJ 


,r 

ij  3xr 


(213) 


r  r 

Since  T..  =  T..,  hence 
iJ  J  i 


(♦.),-  (♦  .)  , 
» J  > 1 


that  is,  the  covariant  differentiation  of  absolute  scalars  is  commutative. 

Having  obtained  the  covariant  derivative,  we  now  obtain  the  Laplacian 
using  (202), 

V2(p  =  div(grad  <j>)  =  g1^  ( — -r-^-r  -  F^  -^r)  (214) 

3x  3xJ  J  3x 

It  is  easy  to  verify  that  if  ^  is  a  curvilinear  coordinate,  xm,  then 


which  was  obtained  earlier  by  another  method  (cf.  (187)).  Similarly, 

if  <p  is  a  Cartesian  coordinate,  x  ,  then 

n 


.  .  32x  3x 

g1J( — : -  rr.  — -)  =  0 

BxV  1J  3xr 


or 


..  32x 

_ n 


dX 


3x13x^ 


-r  V2xr  =  0 


3x 


(215) 


As  an  example,  consider  the  Cartesian  coordinates  x^  =  x,  x^  =  y 
in  a  plane,  and  let  x1  =  £,  x2  =  n  be  the  curvilinear  coordinates.  Then 
introducing  the  operator  D  defined  in  (190),  the  equations  for  x  and  y 
as  dependent  variables  are 


Dx  =  ~g(x^V2C,  +  x^V2n) 

Dy  =  -g(yj.v2c  +  y^v2n) 


(216) 


Equations  (216)  have  been  used  in  Ref.  [29]  to  compute  the  coordinates 
for  arbitrary  shaped  two-dimensional  bodies. 


§7.  Miscellaneous  Derivations. 

In  this  section  we  consider  a  few  derivations  which  are  used  in  the 

study  of  geometry  and  mechanics. 

(i)  Intrinsic  derivative: 
k  k 

Let  x  =  x  (t)  be  the  parametric  equations  of  a  space  curve  with  t 


as  a  parameter.  A  vector  function  u  of  position  will  then  also  be  a 


function  of  t  on  this  curve,  viz.. 


u  -  u(x  (t)) 


The  intrinsic  derivative,  also  called  the  total  or  substantive 

dy 

derivative,  of  u  with  respect  to  t  is  defined  as  —  .  Writing 


u  =  u  a^ 


we  get 


d“  „  du1  i  d~i 

dt  ~dtT  ~i  +  U  “d T 


du1  i  dx£ 

dt  ?i  +  U  3xj  dt 


/du*"  ,  rri  dx^v 
(*dT  +  u  Frj  'dF) -i 


6u 

dt  -i 


The  quantity 


is  called  the  intrinsic  derivative.  We  can  also  write  (217)  as 


Su1  3u*  dx~^  r  i  dx^ 
6t  3xj  dt  +  U  r  J  dt 


i  dx^ 
\j  dt 


(218) 


As  an  example,  in  fluid  mechanics,  the  velocity  vector  u  is  defined 


as 


u  = 


dr 

dt 


where  t  is  the  absolute  time.  Since  dr  =  a^dx  ,  so  that 


dx 

y =  ir 


=  v 


i  dx 

Thus  u  =  are  the  contravariant  components  of  u.  The  components  u* 
can  also  be  explicit  functions  of  t  beside  being  implicit,  but  a^  are  only 
implicit  functions  of  t.  Thus 


61 


dV  _  du1  .  ,  i  d~i 

dt 


dt  dt  “i 


a,  +  a 


j  3u\  irr  j 

U~*  y  *r 


,3u  ,  J  i, 

TtT  +  u  U,j)-i 


(219) 


which  is  the  well  known  substantive  derivative  defining  the  acceleration 
vector. 

The  intrinsic  derivative  of  a  tensor  of  any  structure  can  be  found 
by  using  the  method  followed  in  obtaining  (218).  Thus  for  a  tensor  T, 


.  u,u„*»*u  r,  r 

dT6_12s  1  p 

—  =  -r—  T  a  •••a  a  •  •  *a  r 

dt  fit  r.r  •••r  ~u,  ~u  ~ 

12  pl  s 


Higher  order  intrinsic  derivatives  can  be  obtained  in  a  straight 
forward  manner.  Thus 


fi2u*  6  ,fiu^.  /  i  dx£.  dx^ 

6t2  fit1  fit  '  =  ,j  dt  \k  dt  * 


Intrinsic  differentiation  in  general  is  not  commutative.  Other  uses 
of  intrinsic  derivative  are  in  the  definitions  of  a  parallel  field  of 
vectors  and  of  the  geodesic  curves  in  space. 

A  field  of  vectors  A  along  a  curve  x*  =  x*(t)  and  in  any  space  are 
called  parallel  if  at  all  points  of  this  curve  where 


A  -  A(x  (t)) 
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we  have 


(220a) 


Thus  in  an  Euclidean  space,  the  meaning  of  eq.  (220a)  is  that  the 
components  of  A  referred  to  the  rectangular  Cartesian  coordinates  are 
constants. 

For  arbitrary  coordinates  and  in  any  space,  a  field  of  vectors  is 
called  parallel  when 


dA 

dt 


=  0 


Thus,  the  parallel  field  of  vectors  satisfy  the  equation 


Arr\ 

r  j 


(220b) 


Equation  (220b)  forms  a  system  of  N  equations  in  an  N-dimensional  space 
and  can  be  solved  by  specifying  A*  at  an  initial  point  t  =  t^,  [15]. 

The  geodesics  or  the  geodesic  curves  of  a  space  are  the  curves  along 
which  the  distance  between  two  points  is  minimal.  Let  s  be  the  arc  length 
along  a  curve,  then  we  define  the  unit  tangent  vector  field  t(s)  whose 
contravariant  components  are  given  by 


t(s) 


dx 

ds  ~i 


The  field  t(s)  is  said  to  be  a  tangent  vector  field  on  a  geodesic 


when 


dt 

ds 


0 
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dtiii m 


(221a) 


Using  again  the  definition  (217),  we  obtain  the  equations  for  geodesics 
as 


dxr  dx^ 
ds  ds 


=  0 


(221b) 


Thus  the  geodesics  are  the  solution  of  the  second  order  equations  (221b). 
(ii)  Magnitude  of  a  vector: 

The  magnitude  of  a  vector  u  is  a  scalar.  This  magnitude  is  obtained 
by  taking  the  scalar  product  of  u  with  itself.  Thus 

(u)2  =  u  •  u 
=  (?£  * 

=  gyuV  (222a) 

Also 

(u)2  =  gijuiUj  (222b) 

=  uV  (222c) 

In  the  same  manner,  the  magnitude  of  grad  <f>  denoted  as  |grad  <f>  |  is 
given  by 

|grad  4>  1 2  =  gik  -^r  -^r  (223) 

3x  3xJ 


In  two  dimensions,  writing  x1  =  £,  x2 


n,  we  have 


(224) 


I  grad  $|2  =  [g22  (4>^) 2  -  2gl2d>^4>n  +  gu(*n)2]/g 


(iii)  Angle  between  two  vectors: 

The  angle  0  between  the  two  vectors 


u  =  u  a. 


and 


v  *  vJa. 


is  given  by 


cos  0  =  g 


juW  J (uka 


k)(v  V 


(225) 


The  two  vectors  are  orthogonal  if 


g^uV  =  0 


The  angle  between  any  two  coordinate  curves  at  a  point  is  given  by  the 


base  vectors  corresponding  to  these  curves.  Thus  the  angle  0^  between 


.  i  .  j  . 

the  curves  x  and  x  is  given  by 


cos  Oy  =  (at  •  aJ)//|ai|  |a 


gij/,/gii8jj 


(226) 


where,  since  i  and  j  are  fixed  numbers,  there  is  no  implicit  summation 


on  repeated  indices.  If  x*'  and  x^  are  orthogonal,  then  g^  ■  0  for  i  +  j 


— -  '  Y*  ' *• 
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(iv)  Cross  product  of  vectors: 

For  the  cross  product  of  two  vectors  u  and  v,  the  use  of  eqs.  (90) 
yields  the  result 


j—  i  k  i 

u  X  V  =  /g  e_kuJv  a 


(227a) 


1  ijk 

—  e  u.v  a 

/T  j  k~i 


(227b) 


giving  the  covariant  and  contravariant  components  respectively. 

(v)  Physical  components  of  a  vector: 

In  a  three-dimensional  space  if  all  the  coordinates  are  orthogonal, 
then  as  noted  in  (226) 


g±j  “  0  for  i  +  j 


and  the  non-zero  terms  are  g^,  g22,  g3;j.  It  is  customary  to  use  the 


notation 


hl  =  811  =  ^ 


h2  =  p  =  — — 
2  “22  22 


(228) 


h2  =  g  =  __L_ 
3  °33  33 

8 


(229) 


The  covariant  and  the  contravariant  components  of  a  vector  v 


referred  to  the  orthogonal  coordinates  are  then  related  as 
v  =  h2v3  ,  v2  =  h2v2  ,  v^  =  h3v3 


(230) 


The  physical  components  of  the  vector  v  are  the  orthogonal  projec¬ 
tions  of  the  vector  on  the  coordinate  axes.  Denoting  these  components 
by  V  (subscript  i  is  just  a  label),  we  get 


V1  " 


h,  v 1  =  v 


i/hi 


V 


2 


h^v2  =  v2/h 


2 


(231) 


V 


3 


h^v3  =  v^/h 


3 


The  magnitude  of  the  velocity  vector  v  is  then  simply 

|v | 2  =  Vx2  +  V22  +  V32 

=  V^V1  +  V2V2  +  v^v3 


(232) 


(vi)  Arc  lengths,  elements  of  area  and  volume: 

In  any  coordinate  system  (orthogonal  or  non-orthogonal)  the  arc 
lengths  are  easily  obtained  by  the  metric  equation 

(ds)2  =  gijdxidx^ 

Thus  for  i,  j,  and  k  as  fixed  numbers,  we  have 

i 

! 


(233) 


(ds>xj  =  const>  =  ,/8ii  dx1  (no  summation) 

If 

x  =  const. 

is  an  arc  length  along  the  x1  curve. 

Similarly  denoting  the  element  of  area  on  which  the  curve  x 
const,  as  do^,  we  have 

d°l  =  la2^x2  x  a3dx3| 

=  t822S33  “  (§23^^  1//2dx2<ix3 
dc^  =  |a^dx3  x  a^dx1! 

=  ^811833  ~  (813)2]  1/,2dxldx3 
do^  =  la^dx1  x  ^dx2  | 

=  ^8^^822  “  (812)2]  ^dx^dx2 
The  element  of  volume  is 

dV  =  a^  •  (a^  x  a^)dxJdx2dx3 

=  v'g  dxMx2dx3  . 

§8.  The  Curvature  Tensor  and  Its  Implications. 


(234a) 


(234b) 


(234c) 


(235) 


Questions  regarding  the  nature  of  spaces  have  been  raised  and  dis¬ 
cussed,  mostly  by  philosophers,  at  different  stages  of  human  civilization. 


A  definitive  philosophic  work  on  this  subject  was  published  by  Immanuel 
Kant  in  the  "Critique  of  Pure  Reason"  in  1787.  Despite  a  work  of  such 
brilliance,  the  description  of  space  remained  shrouded  in  mystery  and 
abstract  formalisms.  Scientific  answers  to  the  questions  regarding 
space  started  emerging  after  the  works  of  Gauss  and  Riemann  in  the  first 
half  of  the  nineteenth  century.  In  this  section  we  shall  try  to  define 
a  space  and  its  structure  through  analytic  constructions  as  simply  as 
possible.  The  material  of  this  section  supplements  the  discussions  of 
§2. 

After  gaining  a  working  knowledge  of  basic  tensor  rules  and  par¬ 
ticularly  after  having  the  metric  equation  (74)  at  our  disposal,  we  now 
pose  the  following  simple  problem.  "Is  it  possible  to  devise  a  coor¬ 
dinate  system  x1  =  £,  x2  =  n  in  a  two-dimensional  plane  such  that  the 
element  of  length  between  two  infinitesimally  close  points  be  given  by 
the  metric 


(ds)2  =  (d£)2  +  (cos2£) (dq)2  ?" 


(236) 


In  essence,  the  problem  is  to  find  whether  in  a  two-dimensional  plane 
can  we  have  g^  =  1  and  g^2  =  cos2£? 

The  answer  to  the  above  question  is  that  we  can  never  introduce 
the  above  metric  in  a  plane.  In  fact,  as  we  shall  see  later,  this 
metric  suits  the  surface  of  a  sphere  which  is  a  curved  two-dimensional 
space.  Recall  that  in  a  two-dimensional  plane  we  can  introduce  Cartesian, 
and  plenty  of  orthogonal,  or  non-orthogonal  curvilinear  coordinate  systems. 
Each  chosen  coordinate  system  yields  a  specific  set  of  the  functions  g^ . 


For  example  in  a  plane: 


=  1  ,  g^2  =  0  ,  g22  =  1  ,  for  Cartesian  coordinates, 

g-^  =  1  ,  g12  =  0  ,  g22  =  (£)2  ,  for  polar  coordinates. 

g^  =  1  ,  g^2  =  -2cos  a  ,  g22  =  1  ,  for  oblique  rectilinear 

coordinates  with  a  as  the  included  angle  between  the  coordinates.  We 
can  go  on  adding  to  the  above  list,  but  the  g^j  * s  of  (236)  are  forbidden. 
These  considerations  suggest  that  there  must  be  a  condition  or  a  set  of 
conditions  on  the  g„'s  which  must  be  satisfied  in  each  specific  space. 

To  get  started  on  this  problem,  we  proceed  as  follows. 

Let  A  be  an  arbitrary  vector  and  x1  a  coordinate  system  in  our 
chosen  space  whose  structure  we  wish  to  study.  We  have  the  result  from 
(141)  that  the  partial  derivatives  of  the  entity  A  can  be  expressed  in 
terms  of  the  covariant  derivatives  as 


3A 


3x 


n 


m 


Differentiating  partially  once  more  with  respect  to  x  and  using  (140), 


we  obtain 


where  a  comma,  as  before,  denotes  covariant  differentiation. 


Proceeding  again  from 


we  obtain 


3A 


A.  aJ 


3x 


m  j  ,m~ 


32A 


=  (A.  )  a^  +  rA  A  a^ 

n„  m  j ,m  ,n~  mn  j , l~ 


3x  3x 


Subtracting  (238)  from  (237),  we  get 


32A  32A 


=  .  A„aJ" 


.  m.  n  n,  m  *jnm  i~ 
3x  3x  3x  3x 


(238) 


(239) 


where 


r2  =  _JL_  _  d  p^  +  p^ps  _p2ps 

•jnm  n  jm  m  jn  ns  jm  ms  jn 


It  is  a  direct  algebraic  problem  to  show,  using  eq.  (132),  that 


(240) 


<A. 


n  ,m 


(A. 


m  ,n 


=  R  .  A. 
•jnm  l 


(241) 


H 

The  structure  of  the  quantities  R  .  shows  that  they  are  the 

•jnm 

components  of  a  fourth  order  tensor,  covariant  of  order  three  and  contra- 
variant  of  order  one.  This  tensor  is  known  as  the  Riemann-Chr istof f el 
tensor.  It  is  formed  of  rf,  and  their  first  partial  derivatives.  In 
turn  we  may  state  that  the  Riemann-Chr istof fel  tensor  is  formed  purely 
of  the  metric  coefficients  and  their  first  and  second  partial  derivatives 
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From  (241)  we  conclude  that  the  covariant  differentiation  in  a 


space  is  commutative  provided  that 


(242) 


for  all  values  of  its  indices  and  for  all  coordinate  systems  introduced 
in  the  chosen  space. 

Suppose  in  the  chosen  space  it  is  possible  for  us  to  introduce  a 
set  of  rectangular  Cartesian  coordinates  on  a  global  scale.  The  metric 
tensor  component.'',  are  then  the  Kronecker  deltas  6..  whose  values  are 
either  one  or  zero.  Thus  their  partial  derivatives  and  so  also  all  the 
Christoff el  symbols  are  identically  zero.  The  vanishing  of  all  the 

Christoffel  symbols  makes  (240)  zero  and  so  eq.  (242)  is  satisfied. 

I 

The  vanishing  of  a  tensor  (here  R  .  )  in  one  coordinate  system  means 

•  jnm 

that  all  its  components  should  remain  zero  in  any  other  coordinate 
system  introduced  in  the  same  space.  (Refer  to  the  three  properties  of 
a  tensor  expressions  in  §3.2;  listed  after  eq.  (121)).  It  must  be  noted 
that  when  the  coordinates  are  not  rectangular  Cartesian  then  all  the  g.. 


and  also  the  T.,  are  functions  of  the  coordinates.  Nevertheless,  eq . 
jk 

(242)  will  still  remain  valid.  Spaces  in  which  eq.  (242)  remains  valid 
are  called  Euclidean.  Such  spaces  are  also  called  flat  because  as  will 

j? 

be  seen  shortly,  the  tensor  R  .  determines  the  curvature  of  the  curved 

•  jnm 

space.  Spaces  for  which  eq.  (242)  is  not  satisfied  are  called  Riemannian 
or  non-Euclidean . 

It  is  now  obvious  why  pure  reasoning  fails  to  provide  a  classifica¬ 


tion  of  spaces.  The  idea  of  a  curved  space  is  implicit  not  only  in  the 


(246) 


-  -  ]  U IM*  I  I 


Rrjnp  =^ljP,r]  '^FIjn’r]  +  rjn[rP’"]  "  'jp1™’*1 


,  32g  32g  32g  32g. 

R  =  I( _ 1I£_  +  _ 2IL _ ZIIL  _  _ I1E_) 

rjnp  2  SxV  3xr3xP  3xj3xP  3Xr3xn 


ts 


+g  ([jn,s] [rp,t]  -  [jp,s] [rn,t]) 


(247) 


where  [ij,k]  and  P1  are  the  Christoffel  symbols  of  the  first  and  second 
J  k 

kinds  respectively  as  defined  in  (130)  and  (133). 

From  (247)  it  is  obvious  that 


R  =  — R 
rjnp  jrnp 

R  .  =  -R  . 

rjnp  rjpn 

R  .  =  R 

rjnp  nprj 


(248a) 

(248b) 

(248c) 


and 


R  .  +  R  .  +  R  .  =  0 

rjnp  rnpj  rpjn 


(248d) 


so  that 


R, 

jrnp 


+  R.  +  R. 
jprn  jnpr 


0 


(248e) 


From  (248)  we  also  note  that 


R  4 

rjnp 


0  if  r=j  or,  n=p 
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(249) 


Thus  apart  from  sign,  the  only  non-vanishing  components  are  of  the  form 


R  .  .  ,  R  ,  ,  R  , 

rjrj  rjrp  rjnp 


(250) 


where  r,  j,  n  and  p  are  distinct  from  one  another.  The  total  number  of 
distinct  components  for  a  space  of  dimension  N  are 


-  1)  . 


Thus  the  curvature  tensor  has  only  one  component  in  two-dimensional 
space,  six  in  a  three-dimensional  space,  twenty  in  a  four-dimensional 
space,  and  so  on.  It  can  be  seen  immediately  that  in  a  two-dimensional 
space,  the  component  is  ^212’  a  three-dimensional  space  the 

components  are 


R1212  *  R1313  ’  R2323’  Rl213  *  R1232  *  R1323 


(251) 


(Refer  to  Part  III,  §  2  for  an  expanded  form  of  the  equations  for  a 
three-dimensional  space.) 

It  should  be  noted  that  for  a  flat  space 


R  .  =  0 

rjnp 


(252) 


for  all  values  of  the  indices  r,  j,  n  and  p. 


§8.2.  The  Possibility  of  Local  Cartesian  Coordinates  in  a  Riemannian  Space. 

We  are  now  in  a  position  to  investigate  further  the  curved  nature  of 
Riemannian  spaces.  In  this  section  we  will  show  that  in  a  Riemannian 
space  it  is  possible  to  introduce  a  coordinate  transformation  in  which 
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all  the  metric  coefficients  are  constants  and  all  the  Christoffel 
symbols  are  zero  locally,  (cf.  §2). 

Let  x^  be  a  curvilinear  coordinate  system  in  a  curved  space,  viz., 
not  all  the  components  of  are  zero.  We  now  introduce  a  coordinate 

transformation  from  to  x^  at  a  fixed  point  of  the  space  denoted  by 
subscript  0,  as 


^  =  xj  -  xJ +  •  xo)(x 


s  s. 
-  x0) 


(253) 


5 ■  i + -  # 


6j  =  +  (IJ  )  (xr  _  xr>^ 

m  -to  pr  (T  X(r  — m 
3x  3x 


(254) 


(255) 


Differentiating  (255)  with  respect  to  x  ,  we  get 


a2«J  +  (IJ  ,  iZ  +  (rJ  j  (xr  .  xr)  ?.2»p 

0  a?1  ‘  V0'  0>  aA? 


si's?  pr  0  Si1  3?  'r  0 


<—=>0  ■  -<>0  SK  -  -(rL>0 

3x  3x 


(256) 


Evaluating  (166)  at  x*,  i.e.,  x1  ■+  x*,  and  using  (256),  we  get 


<^>o  ■  (ry>o  -  <>o {?  ■  ° 


(257) 


This  proves  that  the  Christoffel  symbols  at  the  point  x^  in  the  new 
coordinates  are  zero. 

Now  using  (161),  we  find  that 


(gij)0  =  (Sjj  >  o 


(258) 


Further  differentiating  (161)  with  respect  to  x  and  then  using  the 
expression  for  the  derivative  of  g  from  (156b),  we  get 


38ii 

<-£>0  •  0 

dX 


(259) 


The  properties  (257)  and  (259)  are  peculiar  only  to  Cartesian  coordinates 
hence  the  stated  result.  It  must  be  restated  that  the  preceding  results, 
both  for  g  and  rj^,  are  applicable  only  locally  and  not  on  a  global 
scale.  These  results  also  show  that  the  basic  nature  of  a  space  cannot 
be  guessed  simply  by  g_  and  but  by  the  derivatives  of  r^. 

The  coordinate  system  x^  discussed  above  is  also  called  a  geodesic 
polar  coordinate  system. 


§8.3.  Ricci's  Tensor. 

0 

A  contraction  of  l  and  m  in  the  tensor  R  ,  yields  a  tensor  of  the 

•  jnm 

second  order  which  is  called  Ricci's  tensor.  Ricci's  tensor  of  opposite 
sign  will  be  obtained  if  £  and  n  are  contracted.  Thus  Ricci's  tensor  is 


jn 


•jni 


is 

g  R 


sjni 


(260) 


In  expanded  form  it  can  be  represented  in  the  following  two  forms: 


(261) 


32 


Jn  3xJ3x" 


From  (240) 


R.  =  R  , 
jn  *jn£. 


-2-  r* 
n  j*. 


3x 


3x 


L  rl  +  p*1  ps  _  ps 

i.  jn  ns  j£  £s  jn 


(262) 


From  (261)  it  is  obvious  that  the  tensor  R,  is  symmetric, 

jn 

R.  =  R  . 
jn  nj 


since 


(263) 


The  tensor  representation  (262)  is  of  much  importance  in  the  Einstein 

theory  of  relativity  because  it  is  symmetric  and  has  as  many  components 

as  the  metric  tensor  g...  A  scalar  R  can  be  obtained  by  the  inner 

1J 

multiplication  of  g~*n  and  R  ,  viz., 

R  =  gjnRjn  (264) 


and  is  called  the  curvature  invariant. 

§8.4.  Bianchi's  Identity. 

If  we  differentiate  a  second  order  tensor 


T  aU 
Tij-  - 
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and  find  the  partial  derivatives 


3T 


32T 


.  m,  n  ’  .  n.  m 

dX  dX  dX  dX 


then  subtracting  the  two,  we  obtain 


o  0 

/rp  \  „  /'r  \  —  _D  ^  _  D  rjn 

ik,m  , n  ik,n  ,m  ’inm  Hk  ’kran  iH 


(265) 


We  now  take  the  covariant  derivative  of  eq.  (241)  both  sides  and  using 
the  notation 


T  .  =  A.  . 

ij  i.J 


we  write  three  equations  by  cyclic  permutation  as 


(T\  )  ~  (T.  )  =-R*\  A  —  .  T 

jm,n  ,r  jn,ra  ,r  •jnrn.r  l  *jnm  Hr 


(266a) 


(T.  )  -  (T  )  =~R^  A„  -  R*.  T 

jr,m  ,n  jm,r  ,n  •jmr.n  i  ’jmr  Hn 


(266b) 


(T.  )  -  (T,  )  — R*  A  -  Rl .  T 

jn,r  ,m  jr,n  ,in  •jrn,ra  1  *jrn  Hm 


(266c) 


Adding  eqs.  (266),  using  (265)  and  (243b),  we  get 


pee 

R  .  +  R*  +  R*  =0 

•jnp.r  • jpr , n  *jrn,p 


(267) 


This  is  the  first  form  of  the  Bianchi's  identity.  The  second  form  can 
be  obtained  by  using  (244)  and  using  the  fact  that  the  metric  coefficients 
g*^  behave  like  constants  under  covariant  differentiation.  This  form  is 


0 


(268) 


R  ,  +  R  .  +  R  . 

mjnp.r  mjpr.n  mjrn.p 

§8.4.1.  A  Divergence-Free  Tensor. 

The  use  of  Ricci's  tensor  and  the  Bianchi's  identity  produces  an 
important  tensor.  Inner  multiplication  of  (268)  with  gmPg^n  and  use 
of  (248)  and  (260)  yields 


(gjnR.  )  -  (gJnR,  )  -  (gmpR  )  =  0 

jn  ,r  6  jr  ,n  6  mr  ,p 


(269) 


The  first  term  under  covariant  differentiation  in  (269)  is  the  curvature 
invariant  R  defined  in  (264),  so  that 


R 


r 


2(8lV,„ 


o 


If  we  now  introduce  a  mixed  tensor 


(270) 


=  gjnR 


jr 


1  <5nR 

2  r 


then  eq.  (270)  implies  that 


(271) 


Gn  =0  (272) 

•r,n 

That  is,  the  covariant  divergence  of  the  mixed  tensor  is  zero. 

In  place  of  the  mixed  tensor,  we  can  have  a  contravariant  symmetric 


tensor  by  first  writing 


Thus  (270)  becomes 


R 

.r 


Rnq 

,n 


0 


r£ 


Inner  multiplication  by  giA',  and  because  of  symmetry  of  Rnc*,  we  get 


(273) 


where 


Eij  =  Rlj  -  i  giJR  .  (274) 

Equations  (272)  and  (273)  state  that  the  tensor  components  defined  in 
(271)  and  (274)  are  divergence-free.  The  tensor  is  symmetric  and 
is  called  the  "energy-momentum  tensor."  Both  eqs.  (272)  and  (273)  state 
a  conservation  law  of  much  importance  in  physics.  Note  that  the  covariant 
components  of  the  energy-momentum  tensor  can  be  obtained  from  (274) 
by  the  usual  rule  of  lowering  an  index.  Thus 

E«  -  R«  - 1 8ijR  •  (275> 

§9.  The  Geometry  of  the  Event-Space. 

An  event-space  is  a  coordinate  space  in  which  the  time  variable  is 
also  taken  as  one  of  the  coordinates  so  as  to  have  a  space-time  continuum 
in  which  physical  events  occur.  All  the  tensor  theoretic  results  obtained 
so  far  are  obviously  applicable  in  this  space. 

The  geometry  of  the  event-space  has  all  along  been  important  to  the 
theory  of  relativity.  However,  it  is  the  opinion  of  this  author  that  all 


mechanics,  whether  relativistic  or  non-relativistic,  should  be  treated 
at  least  in  the  start  as  a  unified  subject^,  The  rigid  classifications 
of  relativistic  and  non-relativistic  mechanics  deprives  one  from  a 
correct  understanding  of  the  mechanics  and  of  the  geometry  associated 
with  it.  In  this  section  we  first  briefly  state  the  basic  postulates 
of  mechanics  and  go  on  to  explore  some  of  the  consequences  from  a 
geometrical  standpoint. 

An  inertial  frame  of  reference  is  precisely  defined  by  Newton's 
first  law  of  motion.  However,  we  can  rephrase  it  as:  An  inertial 
frame  is  a  coordinate  frame  with  respect  to  which  bodies,  under  the 
absence  of  external  forces,  move  with  zero  acceleration. 

The  two  basic  postulates  of  mechanics  are: 

(I)  All  physical  laws  are  form- invariant  when  transformed  between 
inertial  frames. 

(II)  Light  travels  isotropically  and  with  a  constant  finite 
speed  c(=  2.998  x  10i0  cm/sec)  in  all  inertial  frames. 

In  this  section  we  shall  use  the  Greek  suffixes  for  index  values 
ranging  from  1  to  4  and  Latin  index  values  from  1  to  3. 

Let  there  be  two  inertial  frames  in  which  the  coordinates  are  denoted 
as  xa,  x0.  The  first  three  coordinates  are  the  rectangular  Cartesian 
and  the  fourth  is  the  time.  Thus,  for  example 

x1  =  x  ,  x2  =  y  ,  x3  =  z  ,  x1*  =  t  . 


^The  inspiration  for  the  work  of  this  section  is  due  to  a  paper  by 
G.  C.  McVittie  [30]. 
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Let  the  second  inertial  frame  move  along  x1  with  a  constant  velocity 
V  with  respect  to  the  first.  A  general  coordinate  transformation  between 
xa  and  7*  can  be  written  as  four  equations 


4*  ^  (x1  »x2  »x3  ,  x4) 


(276) 


If  now  the  two  postulates,  and  specifically  of  the  isotropic  propagation 
and  constancy  of  light,  are  used  then  as  shown  by  Tolman  [31]  the  only 

D 

possible  transformation  forms  for  <j>  in  (276)  are 


x1  =  k(x3 

-  Vx4) 

(277a) 

x2  = 

X2 

(277b) 

~3  . 

x3 

(277c) 

x4  =  k(x4 

Vx\ 

c2  } 

(27  7d ) 

where 


V2  -1/2 

k  =  (1  -  -T-)  ' 


The  mapping  (277)  is  called  the  Lorentz  transformation.  It  is  immediately 
seen  from  (27  7),  that  when  V  is  very  small  in  comparison  to  c,  then  k  -*■  1 
and 


x1  =  x1  -  Vx4 


X2  =  x2 


X3  =  X3 


(278a) 

(278b) 


(278c) 


which  is  called  a  Newtonian  or  Galilean  transformation. 

As  a  check  on  the  form-invariancy  of  a  physical  law,  we  can  take 
Newton's  second  law  of  motion  for  a  body  of  constant  mass  m^  and  acted 
upon  by  a  force  system  F1, 

F1  -  ~(m  ~)  ,  i  =  1,2,3  (279) 

dx  0  dx4 

It  can  be  easily  verified  that  under  the  transformation  (278)  the  law 
takes  the  form 

Y1  =  -i-(m  ^)  ,  i  =  1,2,3  (280) 

dx  dx 

so  that  the  form  is  preserved  under  a  Galilean  transformation.  If  we 
repeat  the  same  procedure  using  the  Lorentz  transformation  then  we  find 
that  (279)  cannot  be  transformed  to  the  form  (280) . 

We  now  consider  the  invariant  nature  of  the  element  of  length  for 
an  Euclidean  four-dimensional  space.  Recall  from  eq.  (71)  that  the 
element  of  length  in  E4  will  be 

(ds)2  =  <5  dxadx^ 
up 

This  metric  does  not  remain  invariant  either  for  the  transformation  (277) 
or  (278).  Thus  the  Euclidean  metric  is  completely  unsuitable  for  de¬ 
scribing  a  physical  phenomena  in  an  event  space.  For  the  Lorentz  trans¬ 
formation,  it  can  be  shown  that  the  metric 


(ds)2  =  (cdx4)2  -  (dx1)2  -  (dx2)2  -  (dx3)2  (281) 

transforms  to 

(ds)2  »  (cdx4)2  -  (dx1)2  -  (dx2)2  -  (dx3)2  (282) 

so  that  the  metric  (281)  is  form-invariant  under  Lorentz  transformation. 
We  can  write  (281)  in  the  form  of  eq.  (74), 

(ds)2  =  g  ftdxadx^ 
ap 

so  that 

gn  ~~  — ^  *  ®22  =  *  S33  —  *  ^44  —  ^ 

ga3  =  0  if 

Thus,  for  this  metric 


and  consequently  in  all  the  tensor  formulae  g  must  be  replaced  by  -g. 

To  complete  the  consequences  of  (281),  we  first  write  it  in  the 

form 

(do)2  =  (^)2 

=  (dx4)2 - ^-[(dx1)2  +  (dx2)2  +  (dx3)2]  (284) 

c 

It  must  be  realized  that  the  element  ds  in  (281)  is  not  the  distance 
between  two  closely  spaced  points.  Thus  do  is  an  interval  which  tends  to 
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an  interval  of  time  only  when  c  is  considered  to  be  infinite  rather  than 
finite  as  stated  in  the  second  postulate  of  mechanics. 

Now  writing 


/dx1 . 2 

(tt) 

dx 


+  <^>2  + 

dx  dx 


(u)2 


=  u 


2 


in  (284)  to  have 


dx4 

da 


-1/2 


(285) 


Note  that  in  the  Newtonian  mechanics  c  -*■  °°  and  dx4  =  da,  and  this  is  all 
we  get  out  of  the  metric  (284). 

A  four  dimensional  Minkowski  momentum  vector  is  now  defined  by  the 
components 


m 


dx 
0  da 


,  a  =  1, 2,3,4 


(286) 


where  m^  is  the  mass  of  the  body  when  at  rest,  viz.  u  =  0.  The  fourth 
component  of  (286)  is  defined  as  the  mass  m  of  the  body  in  motion.  That 
is 


m 


0 


dx4 

da 


m 


which  on  using  (285)  gives  the  well  known  relativistic  mass 


The  other  three  components  are 


dx^  _  dx^  dx4 
m0  da  m0  da 


which  are  the  Newtonian  momentum  components. 

The  four-dimensional  Minkowski  force  vector  has  components 


F 


a 


m. 


■  2  a 
d^x 

da2 


0 


d  ,  dx  . 
dx4^ 


Thus 


,,  u2n-1/2  d  ,  dx\ 

(1  ' 


and 


F4 


a  -4>_1/2 

c 


dm 

dx4 


dE 

dx4 


(288) 


(289) 


(290) 


(291) 


where  from  (287) 


(292) 


This  short  discussion  on  the  fundamentals  of  the  special  theory 
establishes  the  connection  between  the  relativistic  and  non-relativistic 
mechanics  and  more  importantly  brings  out  the  structure  of  the  metric 
needed  to  describe  a  space-time  continuum. 

In  the  case  of  general  relativity,  Einstein  proposed  the  principle 
of  covariance  which  states  that  the  physical  laws  under  a  general  trans¬ 
formation  of  coordinates  are  form-invariant.  This  principle  thus 
sweeps  away  the  privileged  position  of  inertial  frames  as  embodied  in 


the  two  principles  of  mechanics.  For  the  description  of  general  rela¬ 
tivity  theory,  the  Riemannian  metric,  eq.  (74),  is  used  in  its  most 


general  form  with  the  metric  coefficients 


as  related  to  the  distribu¬ 


tion  of  matter.  For  details  refer  to  [7  ],  [31],  etc. 


§9.1.  Newtonian  Mechanics  Using  the  Principles  of  Special  Relativity. 

Newtonian  mechanics  with  reference  to  an  inertial  frame  of  reference 
is  described  by  the  three  spatial  Cartesian  coordinates  x1  *=  x,  x2  =  y, 
x3  =  z  at  each  absolute  instant  of  time  x4  =  t.  Thus  time  is  not  affected 
by  the  motion  and  remains  the  same  for  all  coordinate  systems.  In 
essence  time  is  not  a  coordinate  any  more  but  is  a  parameter  which  de¬ 
scribes  the  transformation  of  a  three-dimensional  Euclidean  space  into 
itself  with  the  passage  of  time.  The  geometry  of  this  event-space  is 
then  simply  defined  by  an  Euclidean  E3  metric  in  either  the  Cartesian 
or  any  general  three-dimensional  coordinate  system  with  t  or  x4  as  a 
parameter . 
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If  we  attempt  to  describe  the  motion  of  a  mass  point  as  a  collection 
of  four  numbers  (x1,*2 ,x3,x4)  in  a  four-dimensional  manifold  then  we 
have  to  use  the  metric  of  the  special  relativity  given  in  (284),  with 
the  option  of  performing  a  general  coordinate  transformation  from  the 
Cartesian  to  curvilinear  while  keeping  x4  the  same,  as  described  by 
McVittie  [30]. 

Let  us  introduce  a  transformation  of  spatial  coordinates  to  a 
curvilinear  system  as 


xr  »  c|>raa)  ,  r  =  1,2,3;  a  =  1,2, 3, 4 
x4  =  £4 


(293) 


Thus 


dx 


r 


(294a) 


dx1*  =  d£4 


(294b) 


Writing  (284)  in  the  form 

(do)2  =  (dx4)2 - -  <5  dx'^dx'* 

c^  XJ 

and  using  (294) ,  we  obtain 

(do)2  -  80tedCadC6 

where 

8ag  ^a4^84  c2 
89 


(295) 


(296) 


(297a) 


»** 


and 


\ 


.  3x^  3x^ 

"^44  "  ^il  4  4 

3?  3? 


j.  Sx^  3x^ 

"P4  =  =  «  HP  H4 


-  ($  3x*  3x~^ 

Ypq  "  YqP  '  iJ  HP  3cq  / 


(297b) 


It  must  be  recalled  that  the  Greek,  indices  range  from  1  to  4,  while  the 
Latin  indices  range  from  1  to  3.  Obviously,  as  defined  in  (297b),  the 
metric  (296)  has  symmetric  coefficients,  viz., 

8afS  =  g6a  ' 


Let 


g  =  det(ga6)  ,  A0  =  det(yag) 


and 


(298) 


A  =  Y11Y22Y33  +  2Y12Y13Y23  "  YllYi3  "  Y22Y13  "  Y33Y12 


then 


j^.  ^0 

8  =  c6  c8 


(299a) 


and 
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(299b) 


In  Newtonian  mechanics  we  deal  with  velocities  which  are  much 

IT  Cl 

smaller  than  c.  Thus  x  (£  )  should  be  such  that  y  y  do  not  contain  a 

factor  of  c2.  Hence  g  is  of  the  order  of  one,  while  the  remaining 

g  are  of  the  order  of  -%•  . 
pv  c2 

Let  v  be  a  four  component  velocity  vector,  then  according  to  the 
special  theory  of  relativity  its  ci  . ra variant  components  are 


a  dxa  _  .  (u)2  -l/2  dxa 

U  '  da  1  c2  J  dx4 


For  u  «  c,  neglecting  terms  of  order  — j  ,  we  get 


a  dxa 
U  dx4 


(300) 


Thus,  for  Newtonian  mechanics 


4  ,  i  dx  ,  n  « 

uH  =  l,u  -  —  ,  1  -  1,2,3  . 


In  the  transformed  coordinates  defined  by 


(301) 


i  i,ra. 
x  =  x  (C  ) 


we  have  the  components  as 


v4  -  1  ,  v 


1 


(302) 


1 


The  divergence  of  v  in  special  relativity  is,  (eq.  (201)) 


div  v  =  — - — (/^$T  va) 


/  -  r<* 

/-g  aE, 


(303) 


Using  (299b)  and  neglecting  terms  of  order  -y  ,  we  get 


div  v  =  — - —  (/a  va) 

SE  3  E,a 


(304) 


with  v4  =  1 . 

From  eq.  (211a)  the  divergence  of  a  contravariant  tensor  in  special 
relativity  (four  space)  will  be 


tmv  =  -i-  tmv)  +  rp  Tov 

’v  3CV 


(305) 


Again  using  (299b)  and  neglecting  terms  of  order  —  ,  we  get 


TPV  =  -^^-(/A  TMa)  +  T°a 
,v  ^  oa 


(306) 


§9.1.1.  Application  to  the  Navier-Stokes  Equations. 

As  an  application  of  the  preceding  approximations  of  the  special 
relativity  to  the  Newtonian  mechanics,  we  consider  the  transformation 

of  the  complete r  Navier-Stokes  system  of  equations  to  time  dependent 

+ 

coordinates 1 . 

The  Navier-Stokes  system  of  equations  for  a  viscous  compressible 
fluid  in  the  invariant  vector  form  is 

^McVittie  [30]  has  considered  only  the  inviscid  equations. 

^A  derivation  without  using  the  special  relativity  has  been  obtained  by 
the  present  author,  [32], 


“f 


~  +  div(pv)  =  0 


(307) 


— (pv)  +  div 


=  0 


(308) 


where  p  is  the  density,  t  is  the  stress  tensor,  and  the  div  operator  is 
the  spatial  three-dimensional  divergence. 

We  now  define  a  four-dimensional  energy-momentum  tensor 


such  that 


Tae  =  pv°v6  +  Tae 


v4  =  1 


T44  =  p 


T14  -  t41  =  pv1 


(309) 


(310) 


TiJ  =  TJi 

where  i  and  j  vary  from  1  to  3.  Using  (304)  and  (306)  we  can  write  both 
(307)  and  (308)  as  a  single  equation 

=  0  (311) 

.  p 

As  before,  a  comma  in  (311)  denotes  covariant  differentiation;  a  =  1,  2,  3 
correspond  to  the  three  equations  of  motion  (308)  and  a  =  4  corresponds  to 


the  equation  of  continuity  (307)  . 

Earlier,  in  §8.4.1  the  subject  of  divergence-free  tensors  was  discussed 


and  an  equation  exactly  similar  to  (311)  was  obtained  for  curved  geo¬ 
metries,  eq.  (273).  These  similarities  tell  us  a  lot  about  the  connec¬ 
tion  of  geometry  with  mechanics  and  the  physical  laws. 


t 

h 


Part  II 


The  Geometry  of  Curves  and  Surfaces 
§1.  Theory  of  Curves 

In  this  section  we  shall  consider  only  those  parts  of  the  theory 
of  curves  in  space  which  are  needed  in  the  theory  of  surface.  All  our 
considerations  will  be  confined  to  an  Euclidean  E3  in  which  the  space 
curves  in  parametric  form  are  defined  by 

r  =  r(t)  (1) 

where  t  is  a  parameter  which  takes  values  in  a  certain  interval  a  <  t  <  b. 


c 


Figure  4. 


It  is  assumed  that  the  real  vector  function  r(t)  is  p  >  1  times 
continuously  differentiable  for  every  value  of  t  in  the  specified  interval, 
and  at  least  one  component  of  the  first  derivative 


dr 

dl 


(2) 


is  different  from  zero.  Note  that  the  parameter  t  can  be  expressed  in 
terms  of  any  other  parameter,  say  i,  provided  that  f  0. 
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j L . 


Choose  s,  the  arc  length  along  the  curve  as  a  parameter.  Let 


(5) 


d£ 

k  =  — 
ds 

The  unit  principal  normal  vector  is  then  defined  as 

p  =  k/ |k|  (6) 

The  magnitude 

k(s)  =  |k|  ,  p(s)  =  1/k  (7) 

is  the  curvature  of  the  curve  and  p  is  the  radius  of  curvature.  The 
principal  normal  is  directed  toward  the  center  of  curvature  of  the  curve 
at  that  point. 

(iii)  Normal  plane: 

The  totality  of  all  vectors  which  are  bound  at  a  point  of  the  curve 
and  which  are  orthogonal  to  the  unit  tangent  vector  at  that  point  lie  in 
a  plane.  This  plane  is  called  the  normal  plane. 

(iv)  Osculating  plane: 

Choose  any  three  nearby  points  on  a  space  curve  through  which  a 
plane  can  pass.  Let  the  equation  of  this  plane  be  written  in  the  current 
variable  r  as 

r  •  a  =  c  (8) 

where  a  is  perpendicular  to  the  plane.  We  now  define  a  function  f(u)  of 
the  parameter  u, 

f(u)  =  r  •  a  -  c  (9) 

let  x,  be  the  point  on  the  curve  where  the  parameter  has  the  value 
u..  The  three  points  chosen  on  the  curve  are  denoted  as  u^,  u^,  u^  such 


97 


that  they  satisfy  eq.  (9),  i.e., 

f(u1)  =  0  ,  f(u2)  =  0  ,  f (u3)  =  0 

Hence  according  to  the  Rolle's  theorem 

f ' (Cx)  =  0  »  ui  <  51  <  u2  s 
f ' (C2)  =  0  .  u2  <  ^2  <  u3  ) 


(10) 


(ID 


Because  of  eqs.  (11)  we  can  again  apply  Rolle's  theorem  to  the  function 
f'(u)  in  the  interval  <  u  <  £2>  so  that 

f"U3)  =  o  ,  q  <  ?3  <  s2  (i2) 


As  the  two  points  u2  and  u3  approach  in  the  limit  to  u^,  we  have 

u2>u3>£3>£2»£3  -*■  u^ 

and  equations  (10) ,  (11)  and  (12)  yield 

f(ux)  =  0  ,  f ’ (uL)  =  0  ,  f"(Ul)  =  0  (13a) 


or. 


(13b) 


Combining  equations  (8)  and  (13b),  we  get  the  equation  of  a  plane  at  x^, 

r  =  +  At  +  yp  (14) 

where  A  and  p  are  scalar  parameters.  This  plane  is  called  the  osculating 
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plane,  and  as  shown  by  (14)  it  is  spanned  by  the  unit  tangent  and  the 
unit  principal  normal  vectors. 

(v)  Binormal  vector: 

A  unit  vector  b(s)  which  is  orthogonal  to  both  t  and  p  is  called 
the  binormal  vector.  Its  orientation  is  fixed  by  taking  t,  p,  b  to  form 
a  right-handed  triad  as  shown  in  Fig.  5. 


Figure  5. 


Thus 


b  =  t  x  p  (15) 

Note  that  for  plane  curves  the  binorraal  b  is  the  constant  unit  vector 
normal  to  the  plane,  and  the  principal  normal  is  the  usual  normal  to  the 
curve  directed  toward  the  center  of  curvature  at  that  point. 

The  twisted  curves  in  space  have  their  binormals  as  functions  of  s. 
Because  of  twisting  a  new  quantity  called  torsion  appears,  which  is 
obtained  as  follows. 

Consider  the  obvious  equations 

b*b=l,b*t=0  (16) 


Differentiating  each  equation  with  respect  to  s,  we  obtain 


§1.1.  Serret-Frenet  Equations. 

A  set  of  equations  known  as  the  Serret-Frenet  equations,  which  are 
the  intrinsic  equations  of  a  curve,  are  the  following.  Differentiating 
the  equation 


p  =  b  x  t 


with  respect  to  s,  we  have 


=  xb  -  kt 


Equations  (6),  (18b)  and  (19)  are  the  Serret-Frenet  equations,  and  are 
collected  below 


=  kp  ;  k  =  curvature 


=  -ip  ;  x  =  torsion 


=  xb  -  kt 


For  a  plane  curve,  t  =  0,  so  that 


b  =  constant 


-j—  *  -kt 
ds 
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§2.  Geometry  of  Two-Dimensional  Surfaces  Embedded  in  E3. 

In  the  theory  of  surfaces,  embedded  in  E3,  we  shall  use  u 1 ,  u2 , 
or  ua  as  the  coordinates  in  the  surface  and  x *  (i  =  1,2,3)  as  any 
general  coordinate  system  in  E3.  An  element  of  directed  segment  dr 
is  then  represented  as 

,  ,  i  ,  i 

dr  =  - -  dx  =  a  .dx 


3x 


"i 


(22a) 


If  the  same  element  belongs  to  a  surface  S,  then 


A  J  a  A  01 

dr  = - du  =  r  du 

~  ,,,  a  ~a 

au 


(22b) 


since  in  principle 


r  =  r(xi)  ■  r(xi(u3,u2)) 


for  points  belonging  to  the  surface.  Also  since 


i  i/  i  2^ 
x  =  x  (u  ,uz) 


(23) 


hence 


,  i  9x  ,  a  i,  a 
dx  =  -  du  =  x  du 

-a  a 

rfu 


(24) 


where  we  have  used  the  notation 


i  dx 


(25) 


Now  the  metric  formed  from  (22)  is 


so  that  on  using  (24) ,  we  get 


/ ,  o  i  i,  a  g 

(ds)z  =  XaX&^u  du 


,  a  ,  s 
=  a  .du  du 
aB 


(27) 


The  three  quantities 


a  B 


i  j 
gijXaxB 


(28) 


form  the  components  of  a  symmetric  tensor,  called  the  fundamental  metric 
tensor  of  a  surface.  In  the  old  literature,  the  following  non-tensor ial 
notation  is  also  used. 


11 


12 


=  a 


21 


22 


Since  in  an  Euclidean  space  we  can  always  choose  a  rectangular 
Cartesian  system,  so  that 


8ij 


6 


ij 


and  then 


a 


aB 


,  i  j 

0 . .X 

ij  a  B 


(29) 


From  here  onward  we  shall  return  to  the  symbolism  of  Part  I  and  use 


the  notion 


’ij 


for  a 


ij 


as  there  is  no  chance  for  confusion  if  the  meaning 


is  clear  from  the  context.  Thus,  we  shall  take  the  metric  of  a  surface 


with  coefficients  g  „  rather  than  a  as 

ag  ag 

(ds)2  =  g  dxadx^  (30) 

ap 

Sometimes  when  there  is  no  use  for  an  index  notation  we  have  used  the 
symbols 

x1  =  u  or  5  ,  u2  =  v  or  n 

Therefore 

(ds)2  =  gxi(du)2  +  2g12dudv  +  g22(dv)2  (31) 


The  metric  (30)  for  an  element  of  length  in  the  surface  is  called 
the  "first  fundamental  form"  for  a  surface.  Some  expansions  for  future 
reference  are  listed  below.  A  variable  subscript  in  the  formulae  given 
below  stands  for  a  partial  derivative. 


g  =  x2  +  y2  +  z2  =  r  •  r 
11  u  u  u  ~1  ~1 


(32a) 


g10  =  xx  +  yy  +  zz  =  r  •  r„ 
12  u  v  u  v  u  v  -1  ~2 


(32b) 


=  g 


21 


g,9  =  x2  +  y2  +  z2  =  r 
22  v  v  v  ~2 


-2 


(32c) 


8(v)  =  811822  “  (sl2)‘ 


(33) 


To  distinguish  the  "g"  of  the  general  coordinates  (cf.  Part  I)  from  the 
"g"  formed  by  the  coefficients  of  the  first  fundamental  form  of  a  surface, 


we  denote  the  latter  by  g 
on  the  surface. 


(v)’ 


where  v  stands  for  a  coordinate  held  fixed 


(34a) 


ol&  .a 
ga3S  "  5f3 


gU  =  S22/g(v)  •  S12  =  821  =  -g^/g^A  .  g22  =  S,,/g 


12' 6 (v) 


311/B(v) 


(34b) 


where  v  is  a  parameter  which  remains  fixed  on  a  surface. 


Let  0  be  the  angle  between  r.  =  — —  and  r„  =  .  Then 

~1  3u  ~2  3v 


cos  6  *  (?j_  r2)/'/gii822  8i2^’/gllg22 


I?1  x  =  g11822sin2e 


=  811822(1  "  cos20) 


From  (35)  we  see  that  the  surface  coordinates  are  orthogonal  if 


S12  =  °  * 


The  base  vectors  in  the  surface  defined  in  (22b),  viz., 


-a  .a 

dU 


We  again  emphasize  the  notation  that  a  subscripted  r  such  as  r^  or  r^ 
stands  for  differentiation  with  respect  to  u1  or  u2.  Only  when  expres¬ 
sions  have  been  opened  in  full,  the  notation  u1  =  u,  u2  =  v  has  been 
used . 


define  the  unit  normal  vector  n  at  each  point  of  the  surface  through 
the  equation 


n 


■<v) 


3r 

-( -  X  - ) 

-  3u  3v' 


(37) 


The  Cartesian  components  of  n  will  be  denoted  by  X,  Y, 
from  eq.  (37) 

X  =  Jl//g(v)  ’  Y  =  ’  Z  =  J3/,/8(v) 


and  Z,  so  that 


(38) 


where 


J,  =  y  z  -  y  z 
1  u  v  v  u 


J~  =  x  z  -  x  z 
2  v  u  u  v 


J,  =  x  y  -  x  y 
3  u  v  v  u 


(39) 


§2.1.  Normal  Curvature  of  a  Surface:  Second  Fundamental  Form. 

A  plane  containing  t  and  n  at  a  point  P  of  the  surface  cuts  the 

surface  in  different  curves  when  rotated  about  n  as  an  axis.  Each  curve 

is  known  as  a  normal  section  of  the  surface  at  the  point  P.  Since  these 

curves  belong  both  to  the  surface  and  also  to  the  embedding  space,  a 

study  of  the  curvature  properties  of  these  as  space  curves  also  reveals 

the  curvature  properties  of  the  surfaces  in  which  they  lie. 

We  decompose  the  curvature  vector  k  at  P  of  C,  defined  in  eq.  (5), 

into  a  vector  k  normal  to  the  surface  and  a  vector  k  tangential  to  the 
~n  'g 

surface  as  shown  in  Fig.  7. 
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Figure  7 . 


Thus 

k  =  k  +  k  (40 

~n  ~g 

The  vector  k  is  called  the  normal  curvature  vector  at  the  point  P. 

-n 

is  directed  either  toward  or  against  the  direction  of  n,  so  that 

k  =  n  k  (41 

-n  ~  n 

where  k  is  the  normal  curvature  of  the  normal  section  of  the  surface 
n 

and  is  an  algebraic  number. 

To  find  the  expression  for  k  ,  we  consider  the  equation 

n 

n  •  t  =  0 

and  differentiate  it  with  respect  to  s 


dn 

-r21-  •  t  +  n  •  (n  k  +  k)=0 


Thus  (45)  becomes 


b  duaduB 

k  c  ag 

n  (ds)2 

b  „duaduB 
=  _a3 _ 

,  u ,  v 
g  du  du 
pv 


The  form 

,  ,  ci  ,  3 

b  ndu  du 

ot  3 


(47a) 


(47b) 


(48) 


is  called  the  "second  fundamental  form"  of  the  surface  theory.  The 
expanded  forms  of  b  „  are 


'll 


Xx  +  Yy  +  Zz 


uu 


uu 


uu 


12 


=  Xx  +  Yy  +  Zz  =  b 


uv 


uv 


uv 


'21 


(49) 


22 


Xx 


w 


+  Yy  +  Zz 


w 


vv 


Similar  to  g,  .  ,  we  also  define 
(v) 

b  =  bllb22  '  (b12)2  (50) 

It  is  shown  in  standard  texts  on  differential  geometry,  e.g.,  [6],  [13], 
[17],  etc.,  that  points  on  a  surface  can  be  classified  as  follows. 


b  >  0  ,  elliptic  point 

b  =  0  ,  parabolic  point  (51) 

b  <  0  ,  hyperbolic  point 
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First  notice  that 


We  now  return  to  a  consideration  of  k  . 

n 


n  •  k  =  n  •  (n  k  +  k) 
~  n  ~g 


=  k 

n 


(52) 


Since  p  is  the  unit  principal  normal  to  the  curve,  hence 


k  =  kp 


(53) 


Using  (53)  in  (52)  and  denoting 


n  •  p  =  cos  y 


(54) 


we  obtain 


k  =  k  cos  y 

n 


(55) 


Therefore  if  y  =  0>  then  k  =  k  ;  if  y  =  then  k  =  0  and  the  curve  is 

n  Z  n 

a  plane  curve;  if  y  =  lr»  then  k  =  -k^.  Let 


k  =  -  ,  k  =  — 
p  n  pn 


then 


p  =  p  cos  y  (56) 

n 

Equation  (56)  gives  a  theorem,  called  Meusnier's  theorem:  The  center 
of  curvature  0  of  all  curves  on  S  at  P  having  the  same  tangent  t  lie  on 


Figure  8. 


§2.2.  Principal  Normal  Curvatures. 
From  (47b) ,  writing 


A 


dv 

du 


we  have 


k 

n 


bll  +  2Ab12  +  X2b22 
811  +  2Xg12  +  X2g22 


(57) 


The  coefficients  g  D  and  b  are  constants  at  P,  so  that  k  is  determined 
a  8  a  3  n 

by  the  direction  A.  Thus  all  curves  through  P  having  the  same  tangent 

and  the  same  sense  of  n  have  the  same  normal  curvature  k  . 

n 

To  find  the  extreme  values  of  k  ,  we  differentiate  k  with  respect 

n  n  r 

to  A  and  set  it  equal  to  zero.  Thus 


dk 

TT  =  0  (58) 


The  roots  of  the  above  equation  determine  those  directions  for  which 

the  normal  curvature  k  assumes  extreme  values.  These  directions  are 

n 


called  the  principal  normal  curvatures  at  P  of  the  surface.  To  find 


these  values,  we  first  write  (57)  as 


where 


=  A  +  AB 
n  C  +  AD 


A  =  bH  +  Ab12 


B  "  b12  +  Ab22 


c  "  gll  +  ^g12 


D  =  812  +  Xg22 


On  using  (58) ,  we  get 


B  =  A  _  A  +  AB  = 

D  C  ~  C  +  AD  n 


Elimination  of  A  between  the  two  equations 


B  -  D  k  =0 
n 


A  -  C  k  =0 
n 


gives 


kn  "  g,  ,(gllb22  "  2g12b12  +  g22bll)kn  +  a,  ,  0 


'  ,f¥r 


or 


k2 

n 


,  a$ 

b  .g  k 
a0  n 


=  0 


(59b) 


The  roots  of  the  above  equation  denoted  as  k^  and  k2  are  the  principal 
normal  curvatures.  Obviously 


kl  +  k2  =  baBS 


ctB 


k.  k„  = 


12  g 


(v) 


(60) 

(61) 


Some  definitions  based  on  the  above  derivations  are  given  below. 

(i)  Asymptotic  directions: 

Points  on  a  surface  where  k^  =  0  give  two  directions  through  the 
equation 

f>22<V2  +  2b^2^  +  b^  =  0 
or 


d v 
du 


'b12  1 


*/(b12)z  -  b 


llb22 


22 


(62) 


These  directions  are  called  the  asymptotic  directions.  If  a  straight 
line  can  be  drawn  on  a  surface  then  it  is  obviously  an  asymptotic  curve, 
(ii)  Lines  of  curvature: 

The  line  of  curvature  is  a  curve  whose  direction  is  a  principal 

direction  at  any  of  its  points.  That  is,  at  every  point  of  a  line  of 

curvature  the  normal  curvature  is  either  k^  or  k2-  Thus  the  lines  of 

curvature  are  the  solutions  of  the  equation  (58),  viz., 
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Note  that  the  equation  (63)  is  equivalent  to  two  equations  of 
the  first  degree.  Thus  eq.  (63)  defines  two  families  of  curves  on  a 
surface.  Moreover,  the  two  curves  are  also  orthogonal.  (Refer  to 
eqs.  (123)  for  proof.) 

(iii)  Coordinate  curves  as  lines  of  curvature: 

If  the  curves  u  and  v  on  a  surface  are  lines  of  curvature,  then 
from  (63)  we  have  for 

u  =  const. 

812b22  "  822b12  =  °  ’ 


and  for 


v  =  const. 


8llb12  ~  S12bll  0 


Thus  in  these  coordinates 


g12  =  °  ’  b12  =  ° 


(64) 


Because  of  (64),  (57)  becomes 


for  v  =  const.  ,  (u-curve)  :  k,  =  - 

1  811 

(65) 

^22 

for  u  =  const.  ,  (v-curve)  :  k„  =  - 

2  822 


(iv)  Gaussian  and  mean  curvatures: 

The  product  of  the  principal  normal  curvatures  as  defined  by  (61) 
is  the  Gaussian  curvature  K. 


K  =  k^  =  b/g(y)  (66) 

Similarly  the  mean  curvature  K  is  defined  by  (60)  as 

m 

V  -  X  +  kj)  «7, 


Surfaces  for  which  K  =  0  are  called  minimal  surfaces. 

m 

The  structure  of  the  formula  for  the  Gaussian  curvature  K  given  in 
(66)  shows  that  it  is  an  extrinsic  property.  In  fact  K  is  an  intrinsic 
property  of  the  surface,  viz.,  it  depends  only  on  the  coefficients  g 
of  the  first  fundamental  form  and  their  derivatives.  Refer  to  eq.  (91) 
for  this  aspect  of  K. 

Note  that  if  K  >  0  then  both  the  principal  normal  curvatures  have 
the  same  signs,  while  if  K  <  0  then  they  differ  in  sign.  For  example 
K  >  0  for  ellipsoids,  elliptic  paraboloids  and  spheres,  etc.,  while  K  <  0 
for  hyperbolic  paraboloids,  hyperboloids,  etc. 

§2.3.  Equations  for  the  Derivatives  of  Surface  Normal  (Weingarten 
Equations) . 


Since 


n  •  n  =  1 

hence 


n  •  n  =  0  ,  a  =  1,2 

-a 


(68) 


which  shows  that  lies  in  the  tangent  plane.  Consequently  n^  and  n2 
must  be  linear  in  and  r2- 

n  =  P^1  +  Q~2  (69a) 

n2  =  Rrx  +  Sr2  (69b) 


To  find  P,  ••*,  S,  we  first  note  that 

n  •  r  =  0  ,  n  •  r2  =  0  (70) 

Differentiating  the  first  equation  in  (70)  by  u^  and  the  second  by  u1, 
we  get 


?1 


=  "2 


*1 


(71) 


Thus  from  (46b) 


-bll  =  9i  *  Hi  V 


"bi2  ?i  ' 


?2  • 


b22  -2  ’  ~2 


(72) 


hLr'Tj'i.  -r. 
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I 


|  Since  n  is  orthogonal  to  b  ^  and  r hence  from  (46a)  we  find  that 

I 

!  BaS  =  bag  * 

f 


Next  taking  the  dot  product  of  (75)  with  r^  we  get 


e  =  r  •  r 

a6sy<5  '«6  '6 


(76) 


We  now  write 

?a<3  *  =  [a3’6] 


(77) 


where  as  before  (refer  to  eq.  (129)  of  Part  I),  the  quantities  [a3,6] 
are  called  the  Christoffel  symbols  of  the  first  kind.  If  we  now  take 
the  inner  multiplication  of  both  sides  of  eq.  (76)  with  g°3,  we  get 


=  g6a(ae,6] 


(78) 


The  quantities  defined  in  (78)  are  the  Christoffel  symbols  of  the 
second  kind.  (Refer  to  eq.  (133)  of  Part  I.) 

The  equations 


~a3 


T  „r  +  b  Qn 
aS'Y  aB~ 


(79) 


are  called  the  formulae  of  Gauss  for  the  second  derivatives  r  a. 

~aB 

The  Christoffel  symbols  defined  in  (77)  and  (78)  have  got  exactly 
the  same  structure  as  in  the  general  case  discussed  in  Part  I,  §4. 


Note:  The  symbol  T  is  the  capital  upsilon  of  the  Greek  alphabets. 


However,  because  of  the  two-dimensional  manifold  under  consideration, 
the  Greek  indices  range  only  from  1  to  2.  We  shall  keep  the  notation 
T  in  place  of  r  for  the  Christoffel  symbols  of  the  second  kind  so  as 
not  to  cause  confusion  in  their  use  to  be  discussed  in  Part  III. 


§2.4.1.  Christoffel  Symbols. 

For  future  references,  we  now  list  the  expanded  forms  of  the 
Christoffel  symbols  for  a  surface. 


[aS,  6]  =  + 

3ua 


3g 


a6 


3g 


3u 


cxg. 


3u 


(80) 


Tag  " 


(81) 


Writing  =  u,  u2  =  v,  we  have 


^8-1  i  ^8-1  i  ^8-»  9 

T11  =  tg22  3u  +  812(  3v  ~  2  3u  ^/2g(v) 


(82a) 


tg 


22  l6ll  3v 


9g22  3g22  d&12 

+  g12  3u  2  3v  ^/2g(v) 


(32b) 


T1 

22 


3g12  3822 


[g22(2  3v  "  3u  J  ~  8 


3g22 

12  3v  ^2g(v) 


(82c) 


3g,n  3gn 

T2  =  fe  (2  — — - — ) 

11  L811V  3u  ’ 


9g 


11 


3v 


*12  3u  J/"6(v) 


]/2g, 


(82d) 


T  ^  =  T  1 

12  21 


(g 


3g 


11 


22  3v 


3g22 

‘  8i2  _Tir)/28(v) 


(82e) 
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(82  f ) 


3g 


T2  -  ^2  -  /„ 

12  *21  lgll  3u 


22 


^11 

"  g12  3v  )/2g(v) 


,  Tl 

2g(v)  3u  11 


T2 

12 


(82g) 


1 


2g 


<v) 


d&(»)  „  Ti 

3v  12 


T2 

22 


(82h) 


Note  that  the  g  in  eqs.  (82)  are  those  which  have  been  defined  in  (32) 
ap 

and  (33) . 

From  eq,  (134)  of  Part  I,  we  have  the  result 


(83) 


where 


32r  3?i 

'lj  3x13xj  3xj 


It  is  worthwhile  to  compare  (79)  and  (83)  in  the  same  coordinates  in  the 
sense  that  at  the  surface  both  should  coincide.  This  idea  will  be 
explored  fully  in  Part  III. 

§2.5.  Intrinsic  Nature  of  the  Gaussian  Curvature  (Equations  of  Codazzi 
and  Mainardi). 

The  position  vector  r  in  an  Euclidean  space  can  always  be  represented 
in  terms  of  the  constant  unit  vecotrs.  Thus  it  is  clear  that 
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(84) 


9  ,  x  3  ,  . 

Su’  5“6  '  3u6  ^ 


for  any  choice  of  a,  0  and  y.  We  now  use  eq.  (84)  to  obtain  some 
important  results  of  the  surface  theory. 

On  differentiating  eq.  (79)  and  using  (74c),  while  properly  taking 
care  of  the  dummy  indices,  eq.  (84)  yields 


9T 


[{- 


a8 


3T 


9u ' 


+  Ta  T6 


9u 


-  TA  T6  } 
aB  ya  ay  BX 


~g0(S(b  b  -  b  b  )]r^ 
aB  yo  ay  8a  ~6 


9b  _  9b  . 

+  (——■ - ^  +  T6  b  -  TX  b  )n  =  0 

9uY  9u6  aB  y&  ay  BX  ' 


(85) 


Since  r^  (6  =  1,2)  and  n  are  independent  vectors,  hence  the  coefficients 
of  r^  and  n  must  vanish  separately. 

The  term  in  curly  brackets  in  (85)  is  the  two-dimensional  version 
of  the  Riemann-Christof fel  tensor  defined  in  Part  I,  eq.  (240).  For 
the  sake  of  clarity,  we  use  R*  in  place  of  R  for  the  two-dimensional  case. 
Thus  (85)  yields  the  equations 


R 


*6 

•ay  8 


06 

g 


(b  Qb 
aB  yo 


b  b  )  =  0 
ay  Bo 


(86) 


9b 


aB 


9b 


9u 


.HI 


9u 


T*  b  6 
aB  yo 


-  T>e» 


=  0 


(87) 
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A  two-dimensional  Riemann  curvature  tensor  (similar  to  eq.  (244),  Part 
1)  is  now  introduced  as 


*  *8 
R  .  =  g  -R 
payB  By<S  *ayB 


(88) 


Thus  (86)  becomes 


* 

yayB 


b  b  -  b  b 

aB  yy  ay  By 


(89) 


As  discussed  in  Part  I,  §8.1,  the  covariant  Riemann  tensor  has  in 
all  16  components  in  a  two-dimensional  space  out  of  which,  apart  from 
sign,  only  four  are  non-zero.  Thus  the  four  components  are 


*  *  *  * 
R1212  ’  R2121  ’  R2112  ’  R1221 


where 


*  * 

R1212  =  R2121 


and 


* 

R2112 


* 

R1221 


Therefore  from  (89) 


*  * 
R1212  =  R2121 


(90a) 


*  * 
R2112  =  R1221 


(90b) 
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Using  eq.  (66),  we  get  from  (90a) 


K  R1212/g(v) 


(91) 


Equation  (91)  shows  that  the  Gaussian  curvature  is  an  intrinsic 

k 

property  of  the  surface,  since  R^2i2  f°rmec*  purely  of  the  coef¬ 
ficients  g^g  and  their  derivatives.  This  is  the  "Theorems  egregium" 
of  Gauss. 

Equations  (87)  are  known  as  the  Codazzi-Mainardi  equations.  The 
only  two  possible  equations  from  (87)  are 


3b 


11 


3b 


12 


3v 


3u 


"  6-. ,  T 


11  12 


(TU  -  T12)b12  + 


Tllb22  0 


(92a) 


3b 


22 


3b 


12 


3u 


3v 


+  b„Tl0  +  (T|2  -  T| 0)b 


11  22 


12  12 


T12b22  =  0 


(92b) 


§2,5.1.  A  Particular  Form  of  Codazzi  Equations. 

Consider  the  case  when  u  and  v  are  the  lines  of  curvature,  so  that 

g12  =  °  ,  b12  =  0  •  (93) 

From  (82) : 

Ti  T2  = 

12  2gn  3v  *  11  2g22  3v 

(94) 

2  1  3822  j  1  3g22 

12  =  ?g22  3u  ’  22  2gn  3u 
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The  set  of  eqs.  (92)  take  the  form  (by  the  use  of  eq.  (65)), 


3bll  1  3811 

=  ±(k.  +  k0)—ii 

3v  2  1  2  3v 


(95a) 


3b22  1  3®22 

-nr  -  2<ki +  k2)if 


(95b) 


3  3 

Opening  the  derivatives  '^■(k1811)  and  g^(k2g22^  and  using  (95),  we  8et 


3kl  1  3gll 

nr  -  2^<k2  -  Vir 


(96a) 


3k2  1  3g22 
-(k,  -  k  )-  22 


3u  2g 


22 


2'  3u 


(96b) 


§2.5.2.  The  Third  Fundamental  Form. 

Let  all  the  unit  normal  vectors  n  to  a  surface  have  been  translated 
parallel  to  themselves  such  that  their  initial  points  are  tied  at  the 
origin  of  coordinates.  The  terminal  points  will  then  lie  on  the  surface 
of  a  unit  sphere,  for  ordinary  surfaces.  The  first  fundamental  form  for 
this  sphere  will  then  be 


.  a,  6 
c  „du  du 
ab 


where 


c 


aS 


n 

~a 


(97) 


As  before,  denoting  the  components  of  n  with  respect  to  the  Cartesian 
coordinates  as  X,  Y,  Z,  we  have  the  expansions 
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Cll  =  X2  +  Y2  +  Z2 
11  u  u  u 


c.  „  =  X  X  +YY  +  ZZ 
12  u  v  u  v  u  v 


(98) 


'22 


X2  +  Y2  +  Z2 

V  V  V 


Using  the  Weingarten  equations  (74)  or  (74c),  we  can  also  write  (98)  as 


C11  (g22bll  2g12b12bll  +  gllb12)/g(v) 


C12  [g22bllb12  g12(bllb22  +  b12)  +  gllb12b22]/g 


(v) 


(99) 


c22  (8o')^^2  ~  “®i2b12b22  gllb22^^g 


(v) 


Also 


C  =  C11C22  -  (C12)2  =  b2/g(v) 


(100) 


§2.6.  The  Geodesic  Curvature. 

In  §2,  eq.  (40)  we  wrote  the  curvature  vector  k  as  the  sum  of  the 

normal  curvature  k  and  a  tangential  curvature  k  .  The  vector  k  is 

called  the  geodesic  curvature  vector,  and  its  magnitude  as  the  geodesic 

curvature  k  . 

g 

Since  the  vector  k  lies  in  the  tangent  plane  to  a  surface,  we 
8 

define  a  unit  vector  e  as 


e  =  n  x  t  (Fig.  7) 


(101) 


and  write 
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k  =  e  k 
'8  '  g 


'(102) 


Now 


k  =  e 
g 


A 

•  k 


=  e 


(n  x  t) 


dt 

ds 


=  (5 


n 


(103) 


Further 


t 


dr 

ds 


r 

a 


dua 

ds 


(104a) 


and 


dt  d2r 
ds  ds2 


~ot3 


dua  duB 
ds  ds  ~a 


a2  01 

d^u 

ds2 


(104b) 


Using  (79)  in  (104b)  and  putting  in  (103) ,  we  get  after  some  simplifica¬ 
tions 


“g  ■  ^‘V£>3  -  ^2<S>3  +  «12  -  Til)C^>  = 


dv 

ds 


-  (2T 


1 

12 


t2  \  /dV\  2  d_u  du  d2v  dv  d2u, 
22Mds'  ds  +  ds  ds2  “  ds  d s2 


(105) 


Thus,  on  the  curves  u  =  const.,  and  v  =  const.,  we  have 


<v 


-•s 


u  =  const. 


(v) 


L22 


/g3/2 

°22 


(106a) 
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v  =  const. 


(106b) 


<V 


11 


/g 


3/2 

11 


If  the  coordinates  are  orthogonal,  then 

(k)  =  (107a) 

g  .  / -  du  LL 

u  =  const. 

(k  )  -  -i-  ^(Hn^TT)  (107b) 

g  _  / -  dv  11 

v  =  const.  /g22 

A  curve  C  on  a  surface  S  is  called  a  geodesic  curve  or  simply 
geodesic  if  its  geodesic  curvature  vanishes.  Therefore  for  u  =  const, 
to  be  a  geodesic 

Tj2  =  0  (108a) 

similarly  for  v  =  const,  to  be  a  geodesic 

Tj  -  0  (108b) 

52.6.1.  Geodesics  and  Parallelism  on  a  Surface. 

Having  defined  the  geodesics  as  curves  on  a  surface  whose  geodesic 
curvature  at  each  point  is  zero,  we  must  now  find  the  differential 
equations  of  the  geodesics. 

A  vector  in  the  tangent  plane  of  a  surface  is  known  as  a  surface 
vector,  or  a  vector  in  the  surface.  If  u1  and  u2  are  the  surface  coor¬ 
dinates,  then  according  to  eq.  (22b)  the  vectors  r^  are  the  surface  base 
vectors,  and  they  are  related  with  the  space  base  vectors  a^  as 


r 

~a 


i 


a.x 
~i  a 


i  _  3x 

a  .a 
3u 


Let  A  be  a  surface  vector  field.  Then 


A  =  r  B 

-a 


Since  A  can  also  be  regarded  as  a  vector  field  in  E3,  hence 


A  = 


Thus 


A 


i 


DCt  1 
B  X 

a 


Let  us  consider  a  curve  on  the  surface  whose  parametric  equations 


are 


U  -  U  (t) 


Then 


A  =  A(ua(t)) 


Consequently 


dA  _  d]£  ^  na  du^ 

dt  dt  ~a  5ag  dt 


Using  (79) ,  we  get 


dA 

dt 


[Ba  r  +  Bab  n]~~ 
,  g-a  a&~  dt 


3 


(109a) 
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(109b) 


where  BaQ  is  the  covariant  derivative1  of  Ba  defined  as 
»  p 


> 


0 


3B 

3u£ 


+  BYTa 


Y& 


The  quantity 


._a  ,  B 

6B  _  _a  du 

6t  ,B  dt 


(109c) 


is  called  the  intrinsic  derivative 


of 


A  on  the  curve  u 


a 


a 


u 


(t). 


It  is  interesting  to  note  that 


A 


dA 
*  dt 


B 

a 


6Ba 

fit 


where  B^  are  the  covariant  surface  components  of  A.  Consequently 

£<|A|2>  =>>/)  (109d) 

In  place  of  an  arbitrary  vector  A,  let  us  consider  the  tangent  vector 

Ct  Cl 

field  t(s)  on  a  surface  curve  u  =  u  (s) .  Then 


and 


dt 

ds 


r .  u  .  ^.a,  ,du 
[t  cr  +  t  b  0n]~ 3 — 
,C~a  aB-  ds 


(110) 


The  definition  of  the  covariant  derivative  in  a  surface  is  the  same  as  in 
any  other  space.  Refer  to  Part  I,  §4.  The  only  care  one  should  take  is  to 
replace  T  by  T  for  the  Christoffel  symbols. 
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k 


k  +  k 
~n  ~g 


Since  on  a  geodesic  curve  kj  =  0,  hence  from  (110)  we  find  that  for 


the  geodesics 


f"  ‘ 0 

,S>  ds 


dll) 


Thus  the  geodesics  in  a  surface  are  the  solution  curves  of  the  equations 


d2ua  +  du^  du^ 
ds2  By  ds  ds 


a 


1,2  . 


(112) 


The  definition  of  a  parallel  field  of  vectors  in  a  space  of  any 
dimension  was  given  in  Part  I,  §7,  eq.  (220b).  The  same  definition  is 
applicable  in  the  surface,  viz.,  for  a  parallel  field  of  vectors  A,  the 
intrinsic  derivative  is  zero. 


dB 

dt 


Y  a 
+  BtT 


By 


du 

dt 


6 


(113) 


Equation  (113)  is  also  called  the  condition  of  parallel  displacement  in 
the  sense  of  Levi-Civita.  This  means  that  in  the  covariant  differentia¬ 
tion  the  Christoffel  symbols  T®  are  used.  It  must  be  noted  that  a 
covariant  differentiation  can  also  be  defined  in  which  another  three 
index  symbol,  say  Gg^,  is  introduced,  Weyl  [33]. 

§2.7.  Differential  Parameters  of  Beltrami. 

E.  Beltrami  in  1864  introduced  four  differential  parameters  which 
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greatly  simplify  the  representation  of  some  formulae  in  the  surface 
theory.  If  $  is  a  function  of  the  surface  coordinates  u  and  v,  then  the 
differential  parameters  of  the  first  and  second  orders  are  as  follows, 
(i)  First  order: 


V  =  (gll^  '  2812*U*V  +  822*u)/g(v) 


(114a) 


V*’*)  =  tSnVv  '  812(Vv  +  KV  +  822Vu]/8(v)  (114b) 


(ii)  Second  order: 


3  ,822^u  ~  812<t>vN  .  3  Al^v  ~  812^ 


(v) 


-)  + 

3v  / 


(v) 


(v) 


(114c) 


S22*  -  -  Vu  -  TnV(»w  -  T22»u  "  TL*v> 


<+„v  -  T12*„  -  T?2*v)21/8(v) 


(114d) 


The  parameter  A^<(>  is  the  surface  gradient  of  <j>,  viz.. 


.  .  a8  3<f>  34 

V  =  8  .a  -  8 


3u  3u 


The  parameter  ( 4> , i);)  ifi  related  with  the  angle  0  between  two  curves 
const.,  ij/  =  const.,  as 


(4> , ^)  =  -/ (A1<(»)  (A^)  cos  6 


(115) 


Thus  the  curves  <p  *  const.,  and  ip  =  const.,  are  orthogonal  if 


A1  (<J>  ,if/)  =  0  (116) 

The  condition  equation  (116)  can  be  linked  with  a  second  degree 
equation  in  dv/du  whose  solution  curves  are  orthogonal.  We  recover  this 
theorem  in  a  way  different  from  that  in  Ref.  [  6  ],  p.  80,  as  follows. 

For  p  =  const.,  and  \ p  =  const.,  we  have  the  auxiliary  equations 

A  du  +  d>  dv  =  0 
u  v 


Thus 


ilj  du  +  p  dv  =  0 
u  v 


dv 

du 


(117) 


so  that 


p  lb  -  p  p  =0 
u  V  V  u 


Using  (118)  in  (116) ,  we  have 


hlKK  ~  2g12Vu  +  822^^0  =  ° 


We  now  introduce  the  symbols 


p  p  =  T  ,  p  p  =  S  ,  p  p  =  R 
VTV  V  u  u  u 


so  that  the  condition  of  orthogonality  is 


811T  2§12S  +  g22R  =  ° 


(118) 


(119) 


(120) 


(121) 
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Using  (117)  in  (120),  we  also  have 


T  =  <f>  (~)2 

u  u  dv 


S  =  -<p  ill  ^ 
u  u  dv 


(122) 


R  =  <p  \p 
u  u 


From  (122),  we  find  that 


R(du)2  +  2S  dudv  +  T(dv)2 


is  identically  zero.  We  therefore  state  the  main  result  as  follows. 


"The  ordinary  differential  equation 


R(du)2  +  2S  dudv  +  T(dv)2  =  0 


(123a) 


for  arbitrary  R,  S  and  T  yields  orthogonal  solution  curves  if  and  only  if 


811T  "  2g12S  +  g22R  =  0  *" 


(123b) 


The  second  order  differential  parameter  given  in  (114c)  also  allows 
us  to  define  a  second  order  differential  operator. 


a  _  1  r  3  r  1  /  9  9  ■.  ■, 

A2  "  j - 4u  {  , - Cg22  2u  "  812  3v 


_9  r  1  /  9  9  .  , . 

3v  , - (811  9v  ~  g12  3u)}J 


(124) 


133 


Thus 


A2u 


1  r  3  ,  g22 


^7T8u  Vi 


3  ,  S12 


(v) 


(v) 


)  -  — (• 

-  3v  / 


■)] 


(125) 


(v) 


a,v  -  -i-(f  -  f-<- 

l  1 -  3v  ) -  9u 


3  ,  812 


■)] 


(126) 


’(v) 


In  Part  I,  eq.  (190)  we  introduced  a  second  order  differential 
operator  in  a  plane.  For  a  two-dimensional  surface  it  assumes  the  same 
form 


D  =  g„„3  -  2g  .3  +  g  3  (127) 

22  uu  °12  uv  611  vv 

Using  (125)-(127)  with 

F  =  g(y)A2u  ,  G  =  g(v)A2v  (128) 

we  can  write  (114c)  as 

A24>  =  (D<j>  +  F<J>u  +  G<t»v)  /s(v)  •  (129) 

§2.7.1.  First  Differential  Parameters. 

Let  x,  y,  z  be  the  Cartesian  coordinates.  It  is  easy  to  show  based 
on  the  expansions  (32)  and  the  components  of  the  unit  normal  to  a  surface, 
eq.  (38),  that 

AjX  =>  1  -  X2 
Ajy  =  1  -  Y2 
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A  z  =  1  -  Z2 


A1(x,y)  =  -XY 


A^y.z)  =  -yz 


A1(x,z)  =  -XZ 


If  u  and  v  are  the  surface  coordinates,  then 


Alu  =  8U  =  822/g(v) 


U30) 


V  ‘  b22  '  8ll/8<») 


(131) 


A^u.v)  »  g12  -  -gl2/g 


(V) 


§3.  Mapping  of  Surfaces. 

Let  there  be  two  surfaces  S  and  S  in  which  the  parametric  coordinates 
are  denoted  as  (£,n)  and  (£,n)  respectively.  The  mapping  of  a  portion 
of  S  onto  a  portion  of  S  is  called  a  one-to-one  mapping  when  it  is 
possible  to  establish  the  functional  relations 


£  ■  f1(£,n)  ,  n  =  f2U,n) 


(132) 


where  and  f2  are  differentiable  functions  of  the  desired  orders,  and 
the  Jacobian  of  the  transformation  is  not  zero,  viz.. 


Vn  -  {A  *  °  • 
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In  a  three-dimensional  space  E3,  the  two  points  on  S  and  S,  which 
transform  into  one  another  under  the  mapping  (132) ,  are  respectively 

r  =  (133a) 

r  =  J(?,n)  .  (133b) 

Using  (132)  in  (133b) ,  we  get 

r  =  4<(f1(C,u) ,  f2(C.n)) 


=  x(c» n)  • 

Thus  the  two  points  which  are  the  images  of  one  another  are  representable 
through  the  same  parametric  coordinates  as 


r  =  $U,n) 

r  =  X(?,n) 


(134) 


Equation  (134)  expresses  the  meaning  of  the  sentence,  "the  coordinate 
systems  on  S  and  S  are  the  same. " 

Below  we  discuss  various  mappings  from  one  surface  to  another.  Some 
definitions  have  been  taken  directly  from  Ref.  [18]. 

(i)  Isometric  mapping: 

A  mapping  of  a  portion  S  of  a  surface  onto  a  portion  S  of  a  surface 
is  isometric  if  and  only  if  at  corresponding  points  of  S  and  S,  when 
referred  to  the  same  coordinate  systems  on  S  and  S,  the  values  of  g 
on  S  and  S  are  the  same. 
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Thus  for  isometric  mapping 


8ag^,n)  ~  8ag^,n^ 


(135) 


(ii)  Equiareal  mapping: 

A  mapping  of  a  portion  S  of  a  surface  onto  a  portion  S  of  a  surface 
is  equiareal  if  and  only  if  at  corresponding  points  of  S  and  S,  when 
referred  to  the  same  coordinate  systems  on  S  and  S,  the  values  of  g 


’(v) 


and  g(v)  °f  the  first  fundamental  form  are  equal. 
Thus  for  equiareal  mapping 


8(v)(5’n)  =  8(v)U*n) 


(136) 


(iii)  Geodesic  mapping: 

A  mapping  of  a  portion  S  of  a  surface  onto  a  portion  S  of  a  surface 
is  geodesic  if  and  only  if  at  corresponding  points,  when  referred  to 
the  same  coordinate  systems  on  S  and  S,  the  following  relation  holds. 


ty  _ 


ag 


ag 


+  5 


Y  3f 


a  —  g 

3u 


+  6 


y  3f 


g  -.—a 

9u 


(137) 


where 

f  =  "  tn(g/g)  (138) 

and 

u1  =  ?  =  C 
u'-  =  n  =  n 
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The  proof  of  (137)  follows  directly  from  (105)  equated  to  zero. 

Since  the  proof  of  all  the  statements  in  (i),  (ii) ,  and  (iii)  above  are 
already  available  in  Ref.  [18],  the  reader  is  referred  to  that  work. 

(iv)  Conformal  mapping: 

A  mapping  of  a  portion  of  a  surface  S  onto  a  portion  of  a  surface 
S  is  conformal  if  and  only  if,  when  referred  to  the  same  coordinate 
systems  on  S  and  S,  the  coefficients  g  and  g  are  proportional  at 
each  point ,  viz .  , 

f„6  '  x({-n)8„e  '  (139) 

As  the  name  "conformal"  suggests,  the  angle  between  the  two  inter¬ 
secting  arcs  in  S  is  preserved  in  mapping  to  S. 

(v)  Conformal  mapping  of  surfaces  in  a  plane: 

A  theorem  on  conformal  mapping  of  surfaces  in  a  plane  states  that: 
"Every  portion  of  a  surface  S,  which  is  at  least  three  times  continuously 
differentiable,  can  be  conformally  mapped  into  a  plane." 

In  a  plane  it  is  always  possible  to  introduce  Cartesian  coordinates. 
If  we  denote  these  coordinates  by  u  and  v,  then  we  will  first  show  that 
from  a  general  coordinates  £,n  in  a  surface  in  which  the  metric  is 
given  as 

(ds)2  =  g^^  (d£) 2  +  2g12d£dq  +  g^dn)2  (140) 

we  can  devise  a  transformation  such  that  the  same  ds  is  given  by 

(ds)2  =  A(u,v)[(du)2  +  (dv)2]  (141) 


V-  ,1^,^  - 
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The  coordinates  u  and  v  are  called  the  isothermic  coordinates. 

First  note  that  (140)  can  be  factored  as 

(ds)2  =  [/g^  +  — ~(g12  +  iv'gTJdnH/i^  d£  +  — —  (g12  -  i/gTTdri)] 

where  i  =  /-I .  For  each  term  in  the  brackets  there  exists  an  integrating 
factor.  Let  f^(£,n)  and  f2(£,n)  be  real  functions,  then  we  can  form 
perfect  differentials 

d£ '  =  (fx  +  if2)[>/g^  d£  +  — —  (g12  +  i/g.  Jdn] 

’^n 

dn'  =  (f x  -  if2)[>/gil  dc  +  — zz(g12  -  i4^)dn] 

/gll 


Thus 


(ds)2 


dc'dn* 


£12 


+  f22 


(142) 


The  curves  £'  =  const.,  and  n'  =  const.,  are  called  isotropic  curves. 
Since  £'  and  n'  are  complex  conjugates,  hence 

£'  =  u  +  iv  ,  n'  =  u  -  iv  (143) 


Using  (143)  in  (142) ,  we  get 

(ds)2  =  A(u,v)[(du)2  +  (dv)2] 


where 


We  thus  find  that  a  coordinate  transformation  from  (C,h)  to  (u,v) 


exists  in  which 


822  =  S]_i  ’  gi2  =  0  (144) 

The  above  analysis  proves  the  theorem  of  conformal  mapping  of 
portions  of  S  into  a  plane,  and  also  introduces  the  concept  of  the  iso¬ 
thermic  coordinates.  In  essence,  the  isothermic  coordinates  in  a  surface 
are  those  coordinates  which  are  orthogonal  and  in  which  g^2  =  g^  =  X, 
so  that  the  metric  in  the  surface  is  given  by  (141), 

§3.1.  Isothermic  and  Equiareal  Coordinates  on  a  Sphere. 

We  take  the  parametric  equation  of  a  sphere  of  unit  radius  as, 

(refer  to  Fig.  9  and  eq.  (153)), 

x  =  sin  0  cos  <j>  ,  y  =-sin  0  sin  <p  ,  z  =  cos  0 

where  ^  and  Q  are  measured  clockwise  from  the  x-  and  z-axes  respectively, 

and  0  =  0,  0  =  it  represent  the  north  and  south  poles  respectively. 

0  <  0  <  TT 
0  <  <p  <  2tt 


Figure  9. 

Thus  g^  =  sin20,  =  1,  g^  =  sin20,  so  that  the  metric  on  the  surface  is 


(ds)2  =  sin20(diji)2  +  (d6)2  (144) 

(i)  Isothermic  coordinates  on  a  sphere: 

We  follow  the  technique  shown  previously  from  eqs.  (141)-(143). 

First 

d£’  =  (f^  +  if sinfi  •  d(<£  +  i  in  tan  y) 
dq'  =  (f^  -  if2)sin0  •  d  (<p  -  i  in  tan 


Second 


f  =  — = —  f 

1  sine  »  2 


0 


then 


+  i  in  tan 


=  u  +  iv 


,  e 

n  =  -  i  in  tan  — 

=  u  -  iv 

Equating  the  real  and  imaginary  parts,  we  get  the  mapping 

(145a) 

(145b) 


=  u  ,  0  =  2tan  1 (eV) 


or 


u  =  <\>  ,  v  =  in  tan  y 


n 


1 1 

i! 

'! 


N 


141 
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The  equations  (145)  define  the  isothermic  coordinates  on  a  unit  sphere. 
It  is  an  easy  matter  to  verify  that  using  (145) ,  the  Cartesian  coordi¬ 
nates  are 


i  v 

2e  cos  u 


1  +  e 


2v 


y  = 


-2e  sin  u 


1  +  e 


2v 


1  -  e 
1  +  e 


2V 

2V 


in  which  the  metric  coefficients  are 


(146) 


S22  811 


4e 


2v 


(1  +  e2v)2 


and  the  metric  has  the  required  form, 


(ds)  =  - -S-  ■  -  -[  (du) 2  +  (dv)2] 

(1  +  e2V)Z 

(ii)  Equiareal  (mapping)  coordinates: 

The  mapping 


u 


v  =  2sin  | 


(147) 


(148) 


(149) 


yields  the  metric 

(ds)2  =  +  v2(l  -  ^)  (du)2  (150) 

1  -T 

The  mapping  (149)  is  equiareal,  for  if  we  take  an  auxiliary  Cartesian 
*  * 

plane  x  y  in  which  the  polar  coordinates  are  v  and  u  so  that 

* 

x  =  v  cos  u 


.Li 
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k 

y  =- v  sin  u 

then  the  metric  in  this  plane  will  be  given  by 

(ds  )2  =  (dv)2  +  v2(du)2  (151) 

Thus  the  value  of  g^  in  both  spaces  are  the  same,  which  is  the 
condition  of  equiareal  mapping. 

The  transformation  (149)  establishes  a  one-to-one  correspondence 
between  the  points  of  a  unit  sphere  and  the  points  of  a  plane.  As  9 
varies  from  0  to  it,  v  varies  from  0  to  2.  The  north  pole  is  the  center 
of  the  concentric  circles.  The  limiting  circle  on  the  outside  is  the 


§4.  Some  Standard  Parametric  Representations. 

For  reference  purposes,  we  list  some  parametric  representations  for 
known  surfaces.  In  the  following  we  have  used  ua(a  =  1,2)  to  represent 
the  surface  coordinates.  (Taken  from  Ref.  [18 J). 

(i)  Sphere  of  radius  a: 


r(uu)  =  (a  cos  u2cos  u1,  a  cos  u2sin  u1,  a  sin  u2) 


(152) 


Part  III 


Basic  Differential  Models  for  Coordinate  Generation 

§1.  Problem  Formulation 

The  problem  of  generating  spatial  coordinates,  either  by  analytic 
or  numerical  methods,  is  a  problem  of  much  interest  in  practically  all 
branches  of  engineering  mechanics  and  physics.  A  look  at  the  older 
literature  shows  that  most  problems  in  fluid  mechanics,  electrostatics, 
potential  theory,  space  mechanics,  even  relativity,  etc.,  which  have 
been  classified  as  solutions  of  permanent  value,  are  for  discs,  flat 
plates,  circles,  spheres,  spheoriods,  cones,  ellipsoids,  and  paraboloids, 
etc.  The  main  reason  for  interest  in  these  shapes  is  because  of  the 
availability  of  exact  analytic  coordinates  which  are  body  conforming, 
so  that  the  physical  conditions  at  their  respective  surfaces  can  be 
exactly  imposed.  In  some  cases  the  governing  equations  in  these  coordi¬ 
nates  are  much  simpler  than  in  any  other  coordinates. 

The  coordinates  for  the  above  mentioned  shapes  and  a  score  of 
others  are  obtained  by  the  use  of  elementary  geometrical  and  algebraic 
methods,  which  are  introduced  at  a  very  early  stage  of  one's  mathematical 
and  engineering  education.  Later,  at  a  slightly  higher  level,  in  courses 
on  differential  geometry,  these  coordinates  are  repeatedly  used  in 
exercises  to  investigate  the  geometric  properties  of  surfaces  and  of  the 
curves  which  are  formed  in  them.  These  geometric  properties  are  obtained 
by  using  the  differential  relations  which  have  been  developed  in  Part 
II  of  this  report.  A  question  which  naturally  arises  at  this  stage  is 


.  ■*  .  _ _ 
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this:  Is  it  possible  to  develop  a  set  of  consistent  differential 
relations  and  equations  from  the  available  body  of  differential- 
geometric  results  so  as  to  generate  coordinates  for  arbitrary  shaped 
given  bodies?  In  fact  this  question  has  been  addressed  by  various 
researchers  after  Gauss,  not  from  the  point  of  view  of  arbitrary  shaped 
bodies,  but  for  specific  characteristics  of  a  body.  As  an  example,  the 
most  widely  studied  problem  has  been  of  generating  a  surface,  and  so 
its  coordinates,  when  the  mean  curvature  is  zero  everywhere  in  the 
surface.  Such  surfaces  are  called  the  "minimal"  surfaces.  Weingarten 
surfaces  provide  another  example.  Eisenhart  in  1923  published  a  book 

[34]  on  conjugate  and  other  forms  of  coordinate  net  in  surfaces. 

The  material  of  this  chapter  should  not  be  taken  as  a  review  of 

the  existing  methods  of  coordinate  generation  but  rather  as  an  attempt 
to  bring  in  the  ideas  of  tensors  and  differential  geometry  in  formulating 
problems  of  coordinate  generation.  (A  comprehensive  review  of  the 
existing  methods  in  coordinate  generation  is  to  be  published  shortly 

[35] .)  The  following  two  basic  criteria  have  been  used  in  the  selection 
of  material  for  this  chapter. 

(i)  Derive  only  those  differential  relations  and  equations  which 
have  a  direct  bearing  on  the  geometry  of  the  generated  surfaces,  and 
which  are  of  a  nature  of  lasting  interest  for  future  research. 

(ii)  Methods  to  obtain  solutions  of  the  developed  equations,  if 
possible. 

Two  methods,  which  satisfy  the  above  criteria  are  discussed  below. 
However,  It  is  important  first  to  list  a  few  expansions  for  the  ensuing 
material  from  Parts  I  and  II. 


§1.1.  Collection  of  Some  Useful  Expansions  and  Notation. 

In  what  follows,  the  general  curvilinear  coordinates  are  again 
denoted  as  x'*'.  However,  when  an  expression  has  been  expanded  out  in 
full  and  there  is  no  use  for  an  index  notation  then  we  shall  use  the 
symbols  £,  n,  £,  where 

x1  =  £  ,  x2  =  n  ,  x3  =  5  (1) 


Rectangular  Cartesian  coordinates  are  the  components  of  the  position 
vector  r,  i.e.. 


r  =  (x,y,z) 


(2) 


From  Part  I,  eq.  (39),  the  covariant  base  vectors  in  space  are 

3r 


Thus 


?3 


(3) 


where  a  variable  subscript  will  always  denote  a  partial  derivative. 

The  metric  tensor  g^  in  three  dimensions  has  six  distinct  components. 
The  determinant  g  is  then 


8  *  811822833  +  2g12813823  “  ^823^2gll  ~  ^813^2g22  “  ^812^2g33 


(4) 


Writing 
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(5a) 


G2  =  811833  ~  ^s13^ 2 
G3  =  811822  ~  (g12)2 
G4  =  813823  ~  812g33 
G5  =  812823  ~  813822 
G6  =  812813  ~  811823 


(5b) 

(5c) 

(5d) 

(5e) 

(5f  ) 


we  have 


g11  =  G^g  ,  g22  =  G2/g  ,  g33  =  G3/g  (6a) 

g12  =»  G^/g  ,  g13  -  G5/g  ,  g23  =  G&/g  (6b) 


The  space  Christoffel  symbols  [ i j , k ]  and  T  ,  defined  in  Part  I, 

J  R 

eqs.  (130)  and  (133)  have  been  expanded  for  a  three-dimensional  space 
and  listed  in  Appendix  1.  The  surface  Christoffel  symbols  of  the  second 
kind  for  various  coordinates  held  fixed  are  listed  in  Appendix  2. 
(Refer  also  to  Part  II,  eqs.  (82).) 

As  stated  earlier,  the  x1  or  £,  n,  C  are  the  coordinates  in  a 
three-dimensional  space.  In  place  of  using  different  symbols  for  a  two- 
dimensional  surface  imbedded  in  the  three-dimensional  space,  we  have  used 
£,  n  as  coordinates  on  a  surface  on  which  5  is  held  constant.  This  and 
two  other  possibilities  are  listed  below. 
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const . 


(i) 

Coordinates 

(x3,x2) 

or 

U,n) 

on  a  surface 

on  which 

X3 

=  G 

=  const. 

(ii) 

II 

(x3,xJ) 

or 

(c,U 

II  M 

tl 

II 

(t 

X2 

*  n 

=  const. 

(iii) 

It 

(x2,x3) 

or 

(n.O 

It  It 

II 

II 

It 

X1 

=  5 

=  const. 

(7) 

Note  that  the  right-handed  convention  is  implicit  in  the  ordering  of 
the  coordinates. 


The  index  symbol  v  is  used  in  parentheses  to  denote  which  index  or 
coordinate  has  been  held  fixed,  with  the  exception  of  G  (without 
parentheses)  which  stands  for  the  value  of  G  at  v  =  const,  as  defined 
in  eqs.  (5) . 

For  variations  from  1  to  2,  or  3  to  1,  or  2  to  3,  we  use  Greek 
indices.  Thus,  according  to  Part  II,  eq.  (37),  the  unit  surface  normal 
on  a  surface  v  »  const,  will  be 


s(v)  -  <E„  -  V'ls.  -  r6l 


(8) 


where 


vs*l:o*2,3  =  3  (surface  x1  *  £  =  const.) 


v  *  2  :  a  *  3  ,  6  =  1  (surface  x2  ■  n  =  const.) 


(9) 


v*3:a-l,S=2  (surface  x3  *  £  =  const.) 


The  rectangular  components  of  n^  are 


In  Part  II,  the  G  appearing  here  was  denoted  as  g 


(V)’ 
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(10) 


nM  -  tt(v>,  Y<u>,  z(v))  . 


The  coefficients  of  the  second  fundamental  form  for  a  surface,  b 

ag 

have  been  defined  in  Part  II,  eq.  (46a).  We  now  adopt  the  following 
notation  in  place  of  b 

ap 


=  n^  •  r  (no  sum  on  a) 

-  -act  ' 


T(V>  =  n<V>  •  r 

ag 


(11) 


=  n^  •  roa  (no  sum  on  6) 
'  ~pp 


where  (v,a,g)  are  in  the  permutational  sequences  of  (1,2,3)  as  shown  in 
(9). 

The  Gauss  equations,  defined  in  Part  II,  eq.  (79)  are  now  written  as 


r  =  TY  r  +  S^n^ 
-ota  aa-y 


_  _  „y  _  .  _(v)  (v) 

r  0  -  i  r  -r  i  n 
-ap  ap~y 


(12) 


be  ~  Teeb  +v  - 


where  the  summation  is  to  be  performed  only  on  y,  and  (v,a,g)  are  in  the 
permutational  sequences  of  (1,2,3)  as  shown  in  (9). 

The  sum  of  principal  curvatures  of  the  surface  v  *  const.,  (defined 
in  Part  II,  eq.  (60))  is  now  written  as 
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(13) 


k  (v)  +  k  <v) 
1  2 


<S«/V>  -  2*a6I<V)  +  ^8S<V))/Gv 


where,  in  writing  eq.  (13)  for  a  particular  value  of  v,  use  must  be 
made  of  (9) . 

The  two  second  order  differential  operators  introduced  in  Part  II, 
eqs.  (124)  and  (127),  are  now  written  as 


42<V>  S  iZtS“{i=<lWla  ' 

V  V 


+  3  (g  3  -  2  3  ) }  1 

8 /q-  *cu  e  8a0  or  J 

v 


(14) 


D  ^  =  8003aa  "  2gae9a0  +  gact30B 


(15) 


§2.  Differential  Equations  for  Coordinate  Generation  Based  on  the 
Riemann  Tensor. 


Earlier  in  Part  I,  58,  we  discussed  the  curvature  of  a  general  space 
in  terms  of  the  Riemann  tensor.  It  was  shown  in  Part  I,  eq.  (251)  that 


the  six  distinct  components  of  the  Riemann  tensor  R  ^  for  a  three- 


dimension  space  are 


R1212  ’  R1313  ’  R2323  ’  R1213  ’  R1232  ’  R1323  ’ 

If  the  space  is  Euclidean,  i.e.,  E3,  then  the  above  components  are 
identically  zero  no  matter  which  coordinate  system  is  introduced  in  this 
space.  Thus 


R1212  ”  °  ’  **1313  “  0  *  R2323  “  ° 


(16a) 
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R1213  =  °  ’  **1232  =  °  ’  *1323  =  ° 


(16b) 


Equations  (16)  are  those  consistent  set  of  partial  differential 
equations  which  must  always  be  satisfied  by  the  metric  coefficients  g^ . 

It  should  be  noted  that  there  are  six  distinct  coefficients  to  be  obtained 
from  the  six  equations  (16),  so  that  we  have  a  closed  system  of  equations. 
In  contrast,  a  two  dimensional  space  has  only  one  curvature  equation  and 
three  metric  coefficients,  the  four-dimensional  space  in  general  rela¬ 
tivity  has  twenty  curvature  equations  for  the  ten  metric  coefficients. 

In  these  cases  the  system  is  either  under  determined  or  over  determined, 
respectively. 

Using  eq.  (247)  of  Part  I,  we  now  write  the  six  second  order  partial 
differential  equations  as  dictated  by  (16). 


1212 


3Zgl1  _  2  dZg12  +  &2g22 


3^ 


3£3n 


3?" 


+  2gtS( [22 , s] [11 , t ]  -  [12, s] [12, t])  =  0 


(17) 


1313 


>2gn  3g, 

-  2  - —  + 

3£3S  35" 


gll  „  8  g13  “  633 


+  2gts([33,s] [11, t]  -  [13, s] [13, t])  =  0 


(18) 


82g22  a2g23  92g33 

^2323  "  3;2  “  +  3r)2 


ts 


+  2g  ([33,s] [22, t]  -  [23,s] [23,t])  =  0 


(19) 
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1 


r  _  9  SU  3  812  3  g13  f  9  g23 

1213  ana£  3^3 ?  3£3n  3£2 


+  2gCs([23,s][llft]  -  [12 ,s] {13 , t ] )  =  0 


32g22  32§i2  92g23  +  32g13 


1232  3q3;  a^an  3n2 


ts 


+  2g  ([22,s] [13, t]  -  [23, s] [12, t])  =  0 


R  =  32gl33  _  a2g13  __  32g23  +  a2g12 
1323  3£3n  3n35  3£3?  3?2 


+  2gtS( [33 , s] [12 , t ]  -  [23, s] [13, t])  =  0 


For  a  triply  orthogonal  system  of  coordinates 


812  =  g13  =  g23  “  0 


(20) 


(21) 


(22) 


(23a) 


[12,3]  -  [13,2]  =  [23,1]  =  0 


r\3,  =  r?„  =  rL  -  o 


12  ‘13  23 


(23b) 

(23c) 


=  811§22g33 


(23d) 


Under  the  constraint  of  orthogonality,  eqs.  (17)— (22)  reduce  somewhat. 


Using  eqs.  (23)  and  then  multiplying  the  first  equation  by  1/i/g  g^* 


second  by  1/vg^g^,  and  the  third  by  l/Zg^g^,  we  can  put  the  equations 
in  the  following  form. 


±{—l 

9«V~ 


gll822 


+  !fu.)  * 

34  '  3n V -  ' 


/gllS22  28J3'/gng22 


9gll  dg22 

K  3; 


0  (24) 
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1 


1 


1 


JL(. 

^  r  ' 


3g 


a«Vi 


33)  +  -bi¬ 


ts 


11 


3n 


811833 
1  3§ 


3£  '  3? 


|/811833 


3? 


)  + 


2g22*/giig33 


3gll  3g33 

3n  3n 


=  o  (25) 


A 


”>+4*  1 


3g 


22,  .  1  3g22  3g33 


822833 


3n  '  3c 


A 


822833 


3C 


-)  + 


2g  JZ  g  ^  3? 

811  ®22°33 


-  0  (26) 


9  811  _  1  3gll,  1  3gll  JL.  5g22. 

3n3C  2  3n  K&u  3C  g22  3C  ’  + 


3  g22  x  3g22,_l__  3g22  1  3g33, 

'Nr  "**  _  *\  r  ' 


2  ac  vg22  ac 


g33  3? 


3  g33  ^  3  3g33  i_  38n  3g33, 

'  ~  ■x _  '  «.  x-  /  ' 


3?3n  2  3C  'g13  an 


g33  3n 


1 

3gll 

3g33 

2g33 

3C 

3n 

1 

3gH 

3g22 

2gll 

3C 

H 

1 

3g22 

3g33 

2g22 

H 

3n 

(27) 


(28) 


(29) 


Equations  (24) — (29)  are  the  celebrated  Lame's  equations,  which  he 
obtained  in  1859  by  following  a  different  approach. 


§2.1.  Laplacians  of  £,  q,  and  4  and  Their  Inversions. 

In  Part  I,  eq.  (214),  we  obtained  the  Laplacian  V2<j>  of  a  scalar  <p, 

where 


V2  =  3  +3  +3 

xx  yy  z  z 

From  the  equation  for  V^<p  we  obtained  the  Laplacian  of  any  curvilinear 
coordinate  xm,  the  coordinates  being  assumed  to  be  functions  of  the 
Cartesian  coordinates  x,  y,  z.  Thus 

_  o  m  i i _m  , _  _  . 

V2x  =  -g  JTij  (30a) 
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I<GirU  +  G2r22  +  G3r33  +  2V?2 


+  2G5rl3  +  2G6r23) 


(30b) 


where  x1  =  5,  x2  =  n,  x3  =  5,  and  G^  have  been  defined  in  (5). 

The  inversion  of  these  equations  can  be  written  down  by  using  eq. 
(215)  of  Part  I  as 


,ij  a2* 

’  3xi3xj 


3x 


V23 


ij  l2y  =  _  Jz.  V2xm 


3x^3x^ 


3x 


m 


M  -  2  =  _  J2.  V2xm 


3x13xj 


3x 


m 


Introducing  the  operator 


L  =  G-  3  +  G03  +  Ga 3  +  2G,3r  +  2G,3rr  +  2G,3 

1  2  nn  3  55  4  5n  5  5c  6  nc 


(31a) 


(31b) 


(31c) 


(32) 


we  can  write  eqs.  (31)  as 

Lx  =  -g(x^V25  +  x^V2n  +  x^V2^)  (33a) 

Ly  =  -g(y^25  +  ynv2n  +  y^2t>)  (33b) 

Lz  -  -g(z?V25  +  z^2r\  +  z^V2t)  (33c) 

The  operator  L  reduces  to  the  operator  D  (defined  in  Part  I,  eq.  (190)) 


155 


for  the  two-dimensional  case.  The  corresponding  equations  are  then 
eqs.  (216)  of  Part  I  for  the  surface  £  =  constant. 


§2.1.1.  Laplacians  in  Orthogonal  Coordinates. 

In  the  case  of  orthogonal  coordinates  the  equation  (30b)  can  be 
simplified  to  have  the  following  forms  for  £,  n,  £. 


2  _  _  1  3  ,  /  822833> 


85 


/ 


’ll 


811833^ 

g22 


V2C 


811822 

g33 


) 


(34a) 


(34b) 


(34c) 


§2.2.  Riemann  Curvature  Tensor  for  Specific  Surfaces. 

It  is  worthwhile  for  us  at  this  stage  to  list  the  Riemann  curvature 

tensor  for  specific  surfaces,  £  =  const.,  n  =  const.,  and  £  =  const. 

We  refer  to  Eqs.  (89)  and  (91)  of  Part  II  where  the  single  Riemann 
* 

tensor  for  a  surface  R  „  and  the  Gaussian  curvature  K  were  defined. 

payS 

* 

In  expanded  form  the  expression  for  R^  ^ 

CyB  =  ~"VaP’Ml  "  -VaY.M]  +  T®  luB.o] 

3u  3u 


-  T°g[PY,o] 


(35) 


Thus,  for  £  =  const. 
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R1212  =  Ji[22'1]  ~  ^t21»1]  +  T21[12,1]  +  T21 112,2^ 

-  T^lll.l]  -  T|2[ll,2]  .  (36) 

for  n  *  const. 

R1313  =  "  ^[31,1]  +  T 3! f  13 , 3 ]  +  ^[13,1] 

-  t|3[H»3]  -  ?33 [11,1]  •  ,  (37) 


for  C  =  const. 


* 

R2323 


-  33 » 2 ]  -  ^[32,2]  +  t|2[23,2]  +  t|2  [23,3J 


-  T23[22,2]  -  T|3[22,3] 


(38) 


For  the  expressions  of  T's  refer  to  Appendix  2. 

Each  one  of  (36)-(38)  can  be  reduced  to  different  forms.  For 
example,  eq.  (36)  can  also  be  expressed  as  follows. 


* 

R1212 


-  Tu 


3  ^ 

)  -  TF<r. 


H'g 


11 


T12>1 


(39a) 


3 

=  3 


3  3C  g 


22 


T22) 


3  ^7 

J-C-2  T{ ,)) 

3n  g22  12 


(39b) 


where  is  defined  in  (5c).  The  forms  (39a, b)  are  due  to  J.  Liouville. 
The  Gaussian  curvature  for  each  surface  is 


157 


-  *m2/G3 


k  (n)k  (n) 
*1  k2 


k  (°k  (C) 
1  2 


R1313/G2 


R2323/Gl 


(40a) 


(40b) 


(40c) 


It  must  be  noted  that  in  any  one  of  the  formulae,  from  (36)-(40),  all 
quantities  have  to  be  evaluated  for  the  coordinate  held  fixed.  There 
is  no  difficulty  in  this  process,  since  any  of  these  quantities  have  no 
derivatives  with  respect  to  the  variable  held  fixed. 

As  obtained  in  Part  II,  eq.  (90a),  we  can  also  write  the  equations 
in  (40)  as 


K(S)  =  [S(C)U(C)  _  (T<C))2j/g 


K(n)  =  [s(n)u(n)  _  (T<n))2]/G2 

r(0  =  ^(OyCC)  _  (T(£))2]/G^ 


(41a) 


(41b) 


(41c) 


In  the  representation  (41),  the  quantities  S,  T,  U  can  also  be  determined 
through  other  quantities  which  are  dependent  on  the  derivatives  with 
respect  to  c.  These  representations  thus  establish  a  connection  of  the 
given  surface  with  the  neighboring  surfaces.  This  idea  has  later  been 
used  (see  §3,  eqs.  (74))  to  develop  a  method  of  coordinate  generation 
from  the  data  of  the  given  surfaces. 


§2.2.1.  Coordinates  in  a  Plane. 

If  the  surface  t,  =  const,  (say),  on  which  £  and  n  are  the  parametric 
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coordinates,  reduces  to  a  plane,  then  for  this  surface  =  0.^ 

Consequently,  we  have  a  single  equation  (selecting  either  (39a)  or 
(38b)), 


a  ^  a 

— ( — —  t2  ) — ( — —  T2  )  =  0 

3n'gu  ll;  35'gu  *12; 


In  contrast  to  the  six  equations  (24) -(29)  for  a  three-dimensional 
space,  we  have  only  a  simple  equation  for  two-dimensional  space.  All 
the  three  coefficients  g  cannot  be  determined  from  this  single  equation 
and  additional  relations,  either  algebraic  or  differential,  have  to  be 
imposed  to  solve  eq.  (42).  We  shall  return  to  these  problems  in  §2.4. 


§2.3.  Determination  of  the  Cartesian  Coordinates. 

The  solutions  of  eqs.  (17)-(22)  under  the  prescribed  boundary 

conditions  should  provide  all  the  metric  coefficients  as  a  field 

distribution,  so  that  by  differentiation  one  can  calculate  also  all  the 

Christoffel  symbols  F.,  .  Now  in  any  physical  problem,  e.g.,  the  Navier- 

J  k 

Stokes  equations,  only  the  metric  coefficients  and  the  Christoffel 

symbols  rl,  appear  in  the  transformed  equations,  so  that  the  solutions 
J  k 

of  the  equations  (17)-(22)  provide  all  the  essential  coefficients  to 
solve  the  physical  problem.  Nevertheless,  one  sometimes  also  needs  the 


values  of  the  Cartesian  coordinates  x,  y,  z  as  functions  of  £,  n, 

Our  purpose  is  now  to  describe  a  technique  for  the  determination  of  x, 
y,  z  based  on  the  availability  of  the  metric  coefficients. 

We  define  the  unit  base  vectors 


Refer  also  to  Appendix  3  for  the  Beltrami  equations  in  a  plane. 
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=  =  aj^r± 


(43) 


where  there  is  no  summation  on  i.  The  components  of  along  the 
rectangular  Cartesian  axes  are  denoted  as  u^,  v^,  w^  respectively,  i.e.. 


*i  =  ^i’VV 


(44) 


In  total  there  will  be  nine  values  of  u.,  v.,  w.. 

1  l  l 

Now 


dr  =  d£,  +  ^2*^22  dn  +  ~3'^33  di° 


so  that  the  values  of  x,  y,  and  z  are  given  by  the  line  integrals 


/  ul’/%^  d£  +  u2'/^22  dn  +  U3,/®33  dC 


(45a) 


y  *  /  v1*/i^L  d^  +  v2*^22  dn  +  V3*^33  dC 


(45b) 


z  *  /  w^^/g^"  d^  +  v2'/^22  dn  +  W3V^®33^  d& 


(45c) 


The  determination  of  u^  v  ,  w^  (i  =  1,3)  which  is  needed  in  eqs. 
(45)  poses  another  problem.  Their  derivatives  can,  however,  be  expressed 
in  terms  of  u^,  v^,  w^  by  substituting  (43)  in  eq.  (134)  of  Part  I  as 


!~l=x  far  1  +A  / £22  p2  +A  /!l3  f3 
J  >1 J  gu  U  ~2 J  ij  -3  J  gu  ij 


3x 


’ii 


*i  98ii 


2g 


ii  3xJ 


(46) 
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where  there  is  no  summation  on  i. 

On  changing  i  and  j  from  1  to  3,  we  find  twenty  seven  values  of  the 
derivatives  from  (46).  Thus,  they  form  a  system  of  twenty  seven  first 
order  partial  differential  equations  to  be  solved  under  a  prescribed 
Cauchy  data. 


§2.4.  Coordinate  Generation  Capabilities  of  the  Developed  Equations. 

The  derivation  of  eqs.  (17 ) — ( 22)  has  demonstrated  quite  clearly 
that  these  equations  are  neither  arbitrary  nor  randomly  selected  to 
generate  some  coordinate  system.  They  are  actually  the  very  basic 
equations  which  every  coordinate  system  in  E3  must  satisfy.  Any  six 
symmetric  functions  of  a  coordinate  system  £,  n,  £  which  satisfy 
eqs.  (17)-(22)  are  qualified  to  be  called  the  metric  coefficients  of 
the  introduced  coordinate  system  in  E3. 

Despite  the  versatility  and  power  of  these  equations,  the  solution 
of  these  equations  is  difficult  to  obtain.  The  set  of  equations  (17)- 
(22)  form  a  highly  nonlinear  system  of  coupled  partial  differential 
equations.  Even  if  they  have  been  solved,  the  determination  of  x,  y, 
z  requires  a  solution  of  twenty  seven  first  order  partial  differential 
equations  as  shown  in  §2.3.  Nevertheless,  these  equations  must  form  a 
basis  for  future  developments.  An  indepth  study  into  the  nature  of 
these  equations,  e.g.,  the  comparability  conditions  and  the  type  of 
data  to  be  prescribed  as  the  boundary  conditions,  etc.,  has  to  be 
investigated.  In  the  following  sub-sections  we  consider  two  particular 
cases  of  these  equations. 


§2.4.1.  Two-Dimensional  Orthogonal  Coordinates  in  a  Plane. 

For  the  case  of  orthogonal  coordinates  in  a  plane,  the  basic 


equation  to  start  with  is  eq.  (42).  When  the  constraint  of  orthog¬ 


onality,  viz.. 


812 


=  0 


is  imposed,  we  get  the  equation 


-0 


3g 


/gllg22 


3£ 


22)  +  — c 


3g 


11, 


8ll822 


3n 


(47) 


In  Ref.  [36]  a  method  has  been  developed  to  compute  orthogonal 
curvilinear  coordinates  about  arbitrarily  given  inner  and  outer  boundaries' 
Equation  (47)  is  first  simplified  for  the  case  of  isothermic  coordinates 
(refer  to  53  of  Part  II,  particularly  eq.  (144)  and  the  definition 
that  follows) ,  in  which 


822  “  gll 

and  eq.  (47)  takes  the  much  simpler  form 


92P  3^P 

H2  3n 


=  0 


(48) 


(49) 


where 


P  =  ing 


11  ’ 


From  eqs.  (34a, b)  we  also  have  the  additional  conditions,  that  for 
isothermic  coordinates 
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v2C  =  o  ,  v2n  =  0 


(50) 


where  now  V2  =  3  +3  and  -  =  1. 
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Equation  (49)  can  be  exactly  solved  in  a  Fourier  series  form  by 
prescribing  the  values  of  or  P  at  the  inner  boundary,  denoted  at 
n  =  n^,  and  the  outer  boundary  denoted  as  n  =  as  shown  in  Ref.  [36]. 
The  equations  (50)  are  then  used  to  pick  out  those  £  distributions 
which  establish  an  orthogonal  correspondence  between  the  points  of  the 
inner  and  outer  boundaries. 

§2.4.2.  Three-Dimensional  Orthogonal  Coordinates. 

The  governing  equations  for  the  three-dimensional  orthogonal 
coordinates  are  the  Lame's  equations  and  have  been  stated  earlier  in 
eqs.  (24)~(29).  These  equations  are  as  complicated  as  their  non- 
orthogonal  counterparts.  In  this  section  we  shall  study  two  particular 
forms  of  these  equations  which  are  amenable  to  analysis  and  computation. 

In  this  connection  we  need  the  following  definitions. 

(i)  Gaussian  curvature  in  orthogonal  coordinates: 

For  a  surface  4  =  const,  in  which  f;  and  n  are  the  orthogonal 
coordinates,  and  c  is  the  coordinate  normal  to  the  surface,  we  have 
(cf.  Appendix  2) 


1  3gll  p  1  Sg22 

2g 2 2  ^  ’  12  *  2S22  H 


(51) 


where  G^  *  8^822'  Thus  from  (39a),  we  have 


1  / - ,  ')  ,  i 

2  ’/gllg22  tr/  —  -- 

/gl lg22 


1 


1212 


>R. 


H 


11 


81 ls22 


■)] 


(52) 
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Using  eq.  (24)  in  (52),  we  obtain 


=  R*  /G 
K  R1212/G3 


3gU  3g22 


48n822833 


34  34 


at  4  =  const. 


(53) 


(ii)  Surfaces  of  constant  Gaussian  curvature: 

Surfaces  for  which  the  Gaussian  curvature  has  a  constant  value 
at  every  point  in  the  surface  are  known  as  surfaces  of  constant  curvature. 
The  Gaussian  curvature  can  be  either  zero,  positive,  or  a  negative 
constant  throughout  the  surface.  A  surface  £  =  const,  for  which  =  0 

is  a  developable  surface.  A  developable  surface  can  be  mapped  isomet- 
rically  onto  a  plane.  Recall  from  Part  II,  eq.  (135)  that  the  isometric 
correspondence  between  two  surfaces,  when  the  coordinates  on  the  two 
surfaces  are  the  same,  is  such  that  the  length  element  ds  between  two 
corresponding  points  remains  the  same. 

The  simplest  example  of  a  surface  for  which  >  0  is  a  sphere 

of  radius  R  for  which  =  1/R2.  If  a  sphere  or  a  spherical  cap  can 

be  deformed  in  any  other  shape  whatsoever  without  stretching,  then  its 
Gaussian  curvature  will  not  be  altered.  In  the  case  of  >  0,  every 

surface  of  constant  curvature  can  be  mapped  isometrically  on  a  sphere 
of  radius  (K^)  All  these  results  are  explained  in  books  on  dif¬ 

ferential  geometry,  e.g.,  [17]. 

We  now  consider  the  following  two  cases  of  orthogonal  coordinates. 
Case  I: 


Since  the  coordinates  in  the  surface  are  orthogonal,  the  length 


element  for  C  *  const,  is 


(ds^)2  =  8-qW?)2  +  S22  2 


(54 


We  now  select  C  =  const,  as  any  arbitrary  curve.  Through  every  point 
of  this  curve  a  geodesic  can  be  drawn.  We  call  these  curves  as  n  = 
const.  Obviously  £,  n  are  orthogonal.  From  Part  II,  eq.  (107b)  we 
have  the  result  that  for  n  =  const,  to  be  a  geodesic  we  must  have 


3g 


11 


3n 


=  0  . 


Thus  the  metric  which  we  select  for  Case  I  is 


^11  (^*^  >  *  833  i 


(55 


where  4  remains  fixed  on  each  selected  surface.  Since  is  not  a 
function  of  n,  we  can  define  an  arc  legnth  as  a  perfect  differential 
du  for  each  5  =  const,  as 


du 


Thus 

(ds^)2  =  (du)2  +  g22(dn)2  (56 

In  the  literature,  the  coordinates  u,  n  are  known  as  geodesic  polar 
coordinates,  since  a  point  0  on  the  surface  can  be  chosen  where 


u  =  0  ,  g  ,(0,n)  =  0 


so  that  ds^  =*  0. 


We  now  substitute  (55)  in  (52) ,  and  have 


*  _  _1  / -  a  ,  1  3g22. 

K1212  2  811  822  3uk  , -  3u  ' 

*®22 


.(O  _  l  a  ,  1  °22} 

2,^22  ^  S“ 


Solving  eq.  (58) 


822 


(u  _/k^)  , 


(57a) 


(57b) 


the  parameter  u  being  the  arc  length  along  the  geodesic  coordinates 
n  =  const. 

A  study  of  eq.  (53)  shows  that  for  to  remain  constant,  the 

forms  of  g^^  and  g^2  should  be 


so  that 


gu  =  *U)  f(5) 


&22  =  'Kc)  F(C»n) 


Z  =  const. 


K(5)  =  ( 


V»\ 

4<fi(/ 


The  form  of  the  function  F(£,n)  is  fixed  by  the  solution  (59). 

The  preceding  method  can  be  made  a  basis  of  numerical  coordinate 
generation  for  those  surfaces  which  can  be  isometrically  mapped  on 
spheres  of  varying  radii. 

Case  II: 

In  this  case  we  select  the  metric  such  that  in  the  surface  £  = 
const.,  the  coordinates  are  isothermic.  Thus  we  take 

822  =  811  *  g33  =  1  (60) 


Under  the  constraint  (60),  eqs.  (24)-(29)  simplify  to  the  following 
equations: 


Ifilj  +  _L(J_  +  2  . 0 

H  hi  H  3n  811  3n  2gn  3c 


(61) 


±(~i-  m  o 

8n 


(62) 


!!ii) .  o 

9rT  _  H* 


g 


11 


(63) 


-0 


3C  8  ^ 
gll 


(64) 


while  eq.  (53)  becomes 
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2 


(65) 


K<0  =  1 


A  study  of  eqs.  (62)-(64)  suggests  that  the  only  form  g^  can  have 


Is 


gu  -  (A  +  BU2f(C,n) 


(66) 


where  A  and  B  are  arbitrary  constants,  and  f  >  0. 

Substituting  (66)  in  (61) ,  we  obtain  the  equation  for  f  as 


_LA  lix  +  j_/l  M\  +  2B2f  -  0 

dCf  +  3n(f  3n}  +  2B  f  "  0 


(67) 


Writing 


Q  =  Jtnf 


we  get 


(68) 


which  is  an  equation  similar  to  eq.  (49)  except  for  the  last  term. 

We  now  substitute  (66)  in  (65) ,  and  have 

=  B2/(A  +  BO2  (69) 

Thus,  for  each  5  *  const,  the  surfaces  generated  will  be  of  constant 
Gaussian  curvature.  Numerical  techniques  can  be  used  to  solve  eq.  (67). 
In  the  context  of  isothermic  coordinates  in  the  surface  and  g^j  ■  1, 
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L. 


•  •• 


•'■■■r-r.  ,T,n . 


we  have  the  following  additional  equations  from  eqs.  (34). 


V2£  =  0 

(70a) 

v2n  =  o 

(70b) 

i  3gn 
v24  -  1 

8ll  3? 

(70c) 

With  eqs.  (70)  available,  it  is  possible  to  develop  a  complete 
algorithm  for  numerical  coordinate  generation. 


§3.  Differential  Equations  for  Coordinate  Generation  Based  on  the 

Formulae  of  Gauss. 

In  this  section  we  shall  discuss  another  method  of  coordinate 
generation  suitable  for  three-dimensional  situations  and  which  has  the 
added  property  that  the  method  reduces  to  the  method  of  Ref.  [29]  for 
two-dimensional  plane  regions.  Some  details  are  available  in  a  previous 
publication,  [37]. 

Before  developing  the  proposed  method  it  is  important  to  have  the 
following  formulae. 

From  Part  I,  eq.  (134),  we  have 


r  =  r1  r  +  r  +  rJ  r 
11=5  11-n  ll£c 


■  _  =  r  j  „r  +  r2  r  +  r.3„r_ 

-£n  12-s  I2~n  12 


r  =  r*  r,  +  r2  r  +  r|9r 
~nn  22'C  22~n  22~c 


(71a) 

(71b) 

(71c) 


where  the  3-space  Christoff el  symbols  are  given  in  Appendix  1. 
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We  now  consider  a  surface  designated  as  ;  »  const,  on  which  £  and 


D  are  the  parametric  coordinate.  Then  from  eq.  (12),  we  have 


?££  ll~£  +  il~n  D 


(72a) 


(C)_(C) 


r-  *  Tj.r.  +  T?0r  +  T"*'n 
~£n  12~£  12~n 


(72b) 


r  *  Ti~r  +  T?,r  +  U^n^ 
~nn  22-5  22-n 


(72c) 


where  the  2-space  Christoff  el  symbols  for  £  -  const,  are  given  in  Appendix 

2. 

Taking  now  the  dot  product  of  every  term  with  n^  in  both  eqs. 

(71)  and  (72),  we  obtain 


S(C)  =  (n(C)  ♦ 


(73a) 


T(°  «  (n(0  •  ra)f32 


(73b) 


U(c)  »  (nU)  •  r^)r|2 


(73c) 


All  £ -derivatives  in  eqs.  (73)  are  assumed  to  be  evaluated  at  £  =  const. 

The  above  procedure  can  be  repeated  for  constant  5  and  n  surfaces. 
However,  in  what  follows  we  shall  be  obtaining  formulae  only  for  £  » 
const,  surfaces,  and  for  brevity  of  notation  drop  the  superscript  (£) 
from  the  formulae.  Thus 

S  -  n  •  r?5  -  (n  •  r^rj,  (74a) 
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T  ■  n  •  r 


in 


(n  •  55)r|2 


0  ‘  9  •  Enn  '  <5  •  Es>rJ2 


(74b) 

(74c) 


where 


n  •  r  -  Xx  +  Yy  +  Zz  . 

S  S  S  S 


(74d) 


S3.1.  Formulation  of  the  Problem. 

We  multiply  eqs.  (72a) -(72c)  respectively  by  ~2812‘  81X > 

adding  and  using  eqs.  (13) -(15)  to  have 


Dr  +  G3(r^A2C  +  r^n)  -  +  k2) 


(75) 


where 


D  ”  8223cc  "  2g123cn  +  8ll9nn 


A2C  =  G3(2812T12  "  822T11  '  811T22) 


A2n  "  G3(2812T12  '  822T11  “  gllT22) 


G3  =  811822  “  (g12^2 


(76) 


To  obtain  an  expression  for  +  k2>  we  take  the  scalar  product  of 
(75)  by  n  and  use  eqs.  (74) 


G3(kx  +  k2)  -  (n  •  (811r22  "  2gl2ri2  +  g22ril)  (77) 
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We  now  propose  to  put  constraints  on  the  coordinates  £  and  n  such  that 


A2£  -  0 

(78) 

A2n  *  0 

(79) 

With  these  constraints,  the  differential  equations  for 

of  the  Cartesian  coordinates  are  given  by 

the  determination 

Dr  *  G3n(k1  +  kj) 

(80) 

In  expanded  form,  eqs.  (80)  are 

822x££  ■  2gi2x£n  +  8llxnn  ”  X  ‘  R 

(81) 

822y££  -  2gi2y£n  +  8iiynn  "  Y  *  R 

(82) 

822z££  "  2g12z£n  +  8Hznn  "  Z  ’  R 

(83) 

where 

R  •  +  +  Z-i><Slirl2  -  2*12ri2  +  *22ril> 

(84) 

and 

x  -  (y{*„  -  y„x{)/-^7 

Y  ‘  <Vt  ‘  *?n),/S3 

(85) 

z  -  <Vn  -  ,nyt)/'®3 


The  proposed  constraining  equations  form  the  core  of  the  method. 


Firstly  it  must  be  noted  that  &2  is  neither  a  Laplace  operator  in  the 

Cartesian  plane  (x,y),  nor  in  the  Cartesian  space  (x,y,z),  though  it 

reduces  to  a  two-dimensional  Laplace  operator  when  the  surface  reduces 

to  a  plane  surface,  viz.,  no  dependence  on  z.  Secondly,  the  eqs.  (78) 

and  (79)  express  an  attempt  in  providing  a  set  of  basic  constraints  on 

the  distribution  of  g  .  in  a  surface,  which  is  perfectly  legitimate. 

ap 

Another  important  observation  in  favor  of  these  equations  is  the 
following.  Using  the  expressions  for  the  Christoff el  symbols  appearing 
in  eqs.  (76),  we  can  also  write 


-A+i/fy.)  - 


■S 


(86) 


h2T) 


3n 


7=)  - 


(87) 


In  the  case  of  Isothermic  coordinates,  viz.,  when  g22  ■  g^,  eqs.  (86) 
and  (87)  are  identically  satisfied.  There  is  a  parallel  situation  in 
the  case  of  conformal  coordinates  in  a  plane  where  Laplace  equations  are 
satisfied  identically. 


S3. 1.1.  Particular  Case  of  Eqs.  (81)-(83).  (Minimal  Surfaces). 

For  surfaces  in  which  Isothermic  coordinates  have  been  introduced 
and  at  each  of  its  points  the  mean  curvature  is  zero,  we  have  from  eqs. 
(81)-(83), 
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where  A  and  B  are  arbitrary  constants.  Thus  when  (92)  is  substituted 
for  f  and  g  in  (89)  we  obtain  a  minimal  surface  of  revolution. 


§3.2.  Coordinate  Generation  Between  Two  Prescribed  Surfaces. 

We  now  consider  the  problem  of  coordinate  generation  between  two 
surfaces  denoted  as  n  a  %  and  n  =  n  in  FUg.  11,  where  £  and  £  are  the 

Xj  00 

parametric  coordinates  in  these  surfaces. 

Equations  (81)-(83)  form  a  quasilinear  system  of  partial  differential 
equations  in  which  the  components  of  the  vector  r^  are  assumed  to  be 
prescribed  or  available  through  some  interpolation/extrapolation  numerical 
scheme.  Since  the  values  of  x,  y,  z  are  known  on  the  basic  inner  and 
outer  boundaries  (cf.  Fig.  11),  the  values  of  (r  )  and  (r  ) 

-C  n  =  nB  n  =  noo 

are  known.  Thus  a  suitable  way  of  prescribing  r^  in  space  can  be 


He  ‘  fl<n)(Ee)n  -  „  +f2(n)(r  )n 

iJ  0 


(93) 


where  f^(n)  and  f2(l)  are  suitable  weights  having  the  properties 


fl(V  =  1  ’  f2(V  =  0 

fl<n00)  "  0  ,  f^nj  =  l 


(94) 


Referring  to  Fig.  11,  we  now  solve  eqs.  (81)-(83)  for  each  £  * 
const.,  on  a  rectangular  plane  by  prescribing  the  values  of  x,  y,  and  z 
on  the  lower  side  (C^)  and  upper  side  (C^)  which  represents  the  curves 
on  B  and  °°  respectively.  The  side  (C^)  and  (C^)  are  the  cut  lines  on 
which  periodic  boundary  conditions  are  to  be  Imposed.  The  preceding 
analysis  thus  completes  the  formulation  of  the  problem. 
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§3.3.  Coordinate  Redistribution. 

For  the  purpose  of  generating  coordinates  between  the  space  of  the 
inner  and  outer  boundary,  which  can  be  distributed  in  a  desired  manner, 
we  consider  a  coordinate  transformation  from  £  to  X  and  n  to  o.  Let 


5  =  £(X)  +  £0 

n  =  n(o)  +  nfi 


(95) 


then 


£  =  £0  at  X  =  XQ  ,  £<Xq)  =  0 
n  =  nB  at  a  -  oB  ,  n(a  )  =  0 


(96) 


Writing 


»(X)  ,  8(0) 


(97a) 


and  denoting  the  transformed  metric  tensor  as  ,  we  have 


8h/A2  ’  hi  =  Xx  +  yX  +  Zx 

(97b) 

!/9X  *  *12  =  XXXa  +  Va  +  Va 

(97c) 

gn/e2  ,  i22  =  +  yj  +  Z2 

(97d) 

9  —  8^^822  "" 

(97e) 

X  =  X,Y«Y,Z  =  Z 

(97f ) 

177 


<97g) 


h  +  k2  =  kl  +  k2 


R  =  R/02X2 


(97h) 


Further  noting  that 


r  A. 

(r  -  ^)/X2 
~XX  X 


(98a) 


rc  =  r  /8X 

~5n  -Xo 


(98b) 


r  9 

/  ~q  °\q2 

r  =  (r  -  o 
~nn  ~oo  0 


(98c) 


Substituting  e,s.  (97)  .«d  (98)  In  e,s.  (81)-(83)  »e  no»  have  the 


following  set  of  equations 


x  -  Px„  +  Qx  +  X  •  R 


g22*XX  "  2g12XX0  +  gHXqq  “  «X  +  "Xq 


S22yXX  "  2gl2yXO  +  gllyoo  =  ryx  +  ^q 


=  Py  +  Qy  +  Y  •  R 


(100) 


=  Pz^,  +  Qz„  +  Z  •  R 


g22ZXX  '  2gl2Zxo  +  gHZqo  ZX  0 


(101) 


where 


s22  , 

P  m  -  X 

XX 


611 

0  *  — ^  0 
^00 


(102) 


Thus,  hy  chousing  l  end  0  erbitr.rlly  ve  cen  redistribute  the  coordinates 


-K- 


in  the  desired  manner.  An  example  of  this  choice  is  given  in  the 
next  section. 

§3.4.  An  Analytical  Example  of  Coordinate  Generation. 

In  this  section  we  shall  consider  the  problem  of  coordinate 
generation  between  a  prolate  ellipsoid  (considered  as  an  inner  body) 
and  a  sphere  (considered  as  an  outer  boundary),  with  coordinate  contraction 
near  the  inner  surface.  This  problem  yields  an  exact  solution  of  the 
equations  (99)-(101). 

Let  n  =  n„  and  n  =  n  be  the  inner  prolate  ellipsoid  and  the  outer 
sphere  respectively.  The  coordinates  which  vary  on  these  two  surfaces 
are  £  and  £;.  We  now  establish  a  net  of  lines  made  of  5  =  const,  and  £  = 
const,  on  both  surfaces.  A  curve  on  the  inner  surface  designated  as 

4  =  4q  is 

x  =  cosh  nD  cos 

D  U 

y  =  sinh  rig  sin  Sq  cos  £  (103) 

z  =  sinh  ri  sin  £  sin  £ 

D  U 

Similarly,  the  curve  C ^  corresponding  to  ?  =  ^  on  the  outer  surface  is 


I 


z  =  e  sin  sin  E, 


(104) 


Based  on  the  forms  of  the  functions  x,  y,  z  in  (103)  and  (104),  we 
assume  the  following  forms  of  x,  y,  z  for  the  surface  c,  = 


x  =  f  (o)cos 


y  =  <p  (a)  sin  cos  ? 


z  =  4>(a)sin  sin  £ 


(105) 


The  boundary  conditions  for  f  and  <(>  are 


f(ofi)  =  cosh  ng 


f(o  )  =  e 


<J)  (ct  )  =  sinh  n, 


$(o  )  =  e 


(106) 


Calculating  the  various  derivatives,  metric  coefficients,  and  all 
other  data  needed  in  eqs.  (99)-(101),  we  get  on  substitution  an  equation 
which  has  sin2?0  and  cos2Cq-  Equating  to  zero  the  coefficients  of  sin2^^ 
and  cos24q,  we  obtain 


11  1Lj.1L 

f  0  * 


(107) 


1 


$1  =  11  +  41 

<j>'  e  4> 


(108) 


where  a  prime  denotes  differentiation  with  respect  to  a.  Equations 
(107)  and  (108)  can  be  directly  integrated.  The  solution  under  the 
boundary  conditions  (105)  and  (106)  is 


f(o)  =  A  exp (Bn (o))  +  C 

(109) 

<p(a)  =  D  exp(Bn(o)) 

(HO) 

where 

A  =  (exp(nco)  -  cosh  nfi)  sinh  rig/ (exp (n^)  -  sinh  nB)  (111a) 

l/ (n,  -  n  ) 

B  =  £n[(exp(noo)  -  sinh  ng)  ]  (111b) 

C  =  (cosh  n  -  sinh  n  )exp(n  )/(exp(n  )  -  sinh  n  )  (111c) 

D  iJ  03  oo  g 

D  =  sinh  n  (Hid) 

D 

As  an  application  we  take  the  functions  £(X)  and  n(o)  from  Ref.  [38], 

e(x)  =  ax 

n(o)  =  b(o  -  0g)K° 

where  a  and  b  are  constants.  Since  at  n  , 

CO 

n(o  )  =  n  -  n„ 

oo  00  D 

hence 
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(112) 


(n„  "  O(o  -  a  )  ,  . 

n(o)  -  — B- •_ - K  0  “  V 


a  -  a. 


B 


where  K  >  1  is  an  arbitrarily  chosen  constant.  A  value  of  K  s  1.1 
gives  sufficient  contraction  of  coordinates  near  the  inner  surface. 

For  the  chosen  problem,  since  the  dependence  on  £  is  quite  simple, 
we  find  that  the  coordinates  between  a  prolate  ellipsoid  and  a  sphere 
with  contraction  are  given  by 

x  =  [C  +  A  exp(Bn(o)) ]cos  £ 

y  =  D  exp (Bn (o)) sin  5  cos  C  (113) 

2  =  D  exp (Bn (a) ) sin  5  sin  £ 

where  A,  B,  C  and  D  are  given  in  eq.  (111). 

A  computer  program  based  on  eqs.  (99)-(101)  has  been  developed  by 
* 

Zlebarth  by  using  the  method  of  finite  difference  approximation.  The 
differenced  equations  are  solved  by  using  the  point-successive  over 
relaxation  method.  Complete  duplication  of  the  exact  solution  obtained 
above  has  been  achieved. 


k 

John  Ziebarth,  private  communication. 
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Fjk  “  +  G4[jk,2]  +  G5£ jk,3] } 


Fjk  =  ^{G4[jk’1]  +  G2[jk,2]  +  G6Cjk,3J} 


r?k  =  ^{G5[jk,l]  +  G6[jk,2]  +  G3[jk,3]} 


ril  -  2i{Gl  —  +  G4(2  IT  -  — >  +  G5(2  IT 


rl  =J_{G  !!H+G  ^21  (^fl3  !^23  _  !£l2 

ri2  2giGi  3n  +  G4  ac  G5(  an  +  ac  ac  ; 


■1  i  ^8i  o  ^890  ^S-io 

rl3  -  ilGi  if  +  Vlf  +  if  -  lf>  +  g5  if1 


1  ^8?2  ^®23  ^^22 

r22  -  f  V2  if  -  if) +  g4  if  +  V2  -if-ir" 


rl  _  i{G  (Ifl2  +  !!l3  _  !!23  ^22  ^33} 

r23  2g{Gl(  34  an  ac  '  G4  34  g5  an 


,1  .  J_fG  (2  !f±i  !!llx  +  G  (2  ^23  ^33,  ZM] 
33  "  2gtGl(2  34  a$)+CV2  84  3n  >  *  S  84  J 


3g23  3g33N 


1  1  ^8io  9  ^^11 

ril  =  2gt(j4  ~w  +  G2(2  ~~di  3rT)  +  G6^2  ~H  84”^ 


3g13  3gll, 


f2  ,  J_{G  !!n  !^22  (^13  !^23  _  !il2 

ri2  2g  G4  an  +  g2  ac  +  G6(  an  +  ac  a; 


r  2  „  j_fG  !hi  +  G  (!hi  +  !hi  !hi)  +  G 

r!3  2g  °4  84  +  G2  34  3C  3n  )  S  H 
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3g,  0  Sg™  3®22  ,»  3g23  3g22.. 

r22  *  3(G4<2  Tn - 5T’  +  °2  ~  +  6  3,1  "  Sc 


3g12  .  3g13  9823 


r23  '  2gtG4(  35  +  3n 


.  ^  „  3g22  .  r  1^33 

3^  }  +  G2  "TT  6  3n 


1  3gl3  3g33 

r 2  =  -i-{G,  (2  — rr~  - 

1 33  2g  4  35 


3T>  +  V2  35 


3g23  3833 


3g 


33, 


■)  +  G,  TZ  ) 
3n  6  35 


i  3gU 

r  3  „  J_{GC  ~T  +  G S1 

F11  2gl  5  3£  6 


3gl2  3gU 


U 


3n 


)  +  G_ (2 


3gl3  3gH 


3£ 


35 


)} 


3g 


11 


3g 


r3  =  -i-fG  -  +  G,  ~~T7 

ri2  2g1  5  3n  6  3{> 


22 


3gl3  3823  3gl2 


+  G3(  3n~  +  3^ 


3g,,  3gi 2  3g23  3g13 


11  +  G,(— + 


35 

3g 


)} 


Vh  =  2^{G5  ~W  t  ”6'  35 
,  Sg12  3g22 


35 

3g 


3„  >  +  S  TT1 


r22  "  2g{G5(2  3n 


38, 


3£ 

3g. 


)  +  G 


22 


6  3n  +  G3(2  3t1 


3g23  _  3g22. 
35 


i.„  +  u.^)+o  58 


22 


3g 


33 , 


Is  ‘ 

38, 


1  °13 

r  3  =J_{GC(2  — 

r33  2g  5V 


3n 

3g 


3C 


6  35  +  G3  3h 


33> 


35 


3C-)+G6(2  35  _  3n 


3g23  ^  !f33} 

'  +  G3  35 


For  triply  orthogonal  systems:  r*3  =  r^3  =  ri2 


0. 
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Appendix  2 

Christoffel  Symbols  Based  on  Surface  Coefficients 


Appendix  3 


The  Beltrami  Equations 

For  a  study  of  the  curvilinear  coordinates  in  plane  two-dimensional 
regions  the  technique  of  quasiconformal  mapping  is  frequently  used. 
Quasiconformal  mappings  are  more  general  and  flexible  than  the  usual 
conformal  mappings.  For  details  on  the  mathematical  aspects  of  the 
quasiconformal  mappings,  refer  to  [39]  and  [40]. 

* 

A  quasiconformal  mapping  of  a  region  D  onto  a  region  D  is  given  by 
a  one-to-one  continuous  mapping  whose  inverse  is  also  continuous 
(homeomorphism) .  The  mapping  function  w  =  f(z,z)  for  this  case  is  taken 
as  a  solution  of  the  complex  equation 

f--H(z,I)f  =0  (1) 

z  z 


where 


z  =  x  +  iy,z  =  x-iy,  i  =  /-L  . 


The  complex  equation  (1)  is  called  the  Beltrami  equation,  which  is 
equal  to  the  two  real  equations 


-n  =  f?C  +  VC  (2) 

x  x  y 

n  =  <*C  +  6C  (3) 

y  *  y 


where 


Y  =  [(1  +  u)2  +  v2]/A 


A  =  1  -  (p2  +  v2) 


(4e) 

(4f) 


Note  that 


ay 


S32  =  1 


a  +  y  =  2(2  -  A) /A 


(6) 


(7) 


A  quasiconformal  mapping  becomes  conformal  when  H  =  0,  or,  equivalently 
a=y=l,3=0.  In  this  case  eqs.  (2)  and  (3)  reduce  to  the  Cauchy- 
Riemann  equations 


=  \  s  = 


and  f(z)  is  then  a  holomorphic  or  analytic  function  in  D. 
Now,  from  eq.  (1) 


|  H | 2  =  p2  +  v2  =  |f_/f  | 2 

11  1  Z  2  1 


so  that  on  using  eqs.  (175)— (177)  of  Part  I,  we  obtain 


A  =  4/g/[2/g  +  (g^  +  g22)] 


(8) 


(9) 


1 


a 


Substituting  (9)  in  (7),  we  get 


.  . 
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Equations  (2)  and  (3)  can  also  be  written  by  using  the  inversion 
relations  given  in  Part  I,  eq.  (177),  as 

x_  =  ay  -  gx  (11) 

£  n  n 

yc  =  %  -  yxn  (12) 

Solving  eqs.  (11)  and  (12)  for  and  y^,  we  have 

x  =  8x  -  ay  (13) 

n  £  £ 

yn  =  yxc  -  ey5  (14) 

The  Beltrami  equations  (13),  (14)  form  a  system  of  first  order 
partial  differential  equations  for  numerical  coordinate  generation. 

The  coefficients  a,  B  and  y  are  related,  as  can  be  seen  by  solving  eqs. 
(6)  and  (10), 

a.Y  =  [gu  +  g22  +  (fe11  +  §22^  2  ~  +  ^  I//2]/2v/g  (15) 

The  choice  of  8  can  be  based  on  the  minimization  of  a  certain  functional 
to  ensure  uniqueness.  This  algorithm  has  been  followed  in  Ref.  [41]. 

If  orthogonal  coordinates  are  desired,  then  using  eqs.  (13)  and  (14) 


in  the  orthogonality  condition 


we  obtain  3  through  the  algebraic  equation 
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